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THE FREEZING POINTS OF HIGH PURITY METALS AS 
PRECISION TEMPERATURE STANDARDS 


VI. THERMAL ANALYSES ON FIVE SAMPLES OF LEAD WITH PURITIES 
GREATER THAN 99.999+-©;! 


E. H. McLAREN AND E. G. Murpbock 


ABSTRACT 


An investigation has been made of the freezing and melting temperatures of 
five samples of high purity lead (supplier’s analyzed impurity contents < 0.7 to 
< 4 p.p.m.) including zone refined metal. Using the induced freezing technique, 
plateaux of essentially constant (< +0.0001° C) temperature with durations of 
over 1 hour are readily obtained on the cooling curves of these samples. A standard 
deviation in plateau temperature (liquidus point) of +0.0001° C was obtained 
from a series of 30 induced freezes on a particular sample. The pressure effect 
on the freezing temperature of lead was found to be 0.0080° C for 1 atmosphere. A 
value of 327.426° C (Int. 1948) was determined for the standard liquidus point 
of pure lead. 

The liquidus points of the samples were intercompared with a precision of 
about 0.0002° C, and alloy melting ranges were examined following different 
types of freezing with and without overnight anneals near the solidus tempera- 
ture. Alloy melting range parameters were found to be useful in the selection 
of the samples of highest purity and at the same time showed that an uncertainty 
of 0.002° C in the above value of the liquidus point of pure lead may exist because 
of residual impurity contents in the purest samples that were examined. 


INTRODUCTION 

An investigation of the freezing and melting temperatures of five samples 
of high purity lead has been carried out to determine the suitability of this 
transition as a thermometric fixed point and to examine the usefulness of 
precision resistance thermometry to discriminate among samples of high purity 
lead. 

MATERIAL 

Table I lists the origin, grade, and form of the lead metal as received, and 
Table II shows the chemical analyses supplied by the sources. Since the lead- 
rich portions of the respective binary phase diagrams for these analyzed 


impurity solutes appear to be reasonably well established, estimates, listed 
in Table III, were made of the expected lowerings of the liquidus and solidus 
from that of pure element and the associated alloy melting ranges for the 


samples. 


‘Manuscript received January 27, 1960. } 
Contribution from the Division of Applied Physics, National Research Council, Ottawa, 
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TABLE I 
Identification of the lead samples 











Sample Origin Grade Form 

Pb.3 The Consolidated Mining High Purity No. 50 Bar, 2.5X2.5X25 cm 
and Smelting Co., 

Pb.4 Trail B.C. High Purity Bar, 2.5X2.5X25 cm 
Canada (ZR-2) No. 186 

(Zone Refined) 

Pb.5 American Smelting Special High Purity Bar, 1.2 triangular sec- 
and Refining Co., Lot HPL No. 7 tion X25 cm 
South Plainfield, N.J., U.S.A. 

Pb.6 Philipp Bros. Zone Refined Bar, 1.2X1.0X30 cm 


(Canada) Ltd., Montreal 25, (Japanese) 

Que., Canada 
Pb.7 Associated Lead Special High Purity Small bars 
Manufacturers Ltd., No. HPEL/2 
Middlesex, England. 


TABLE II 
Impurity contents of the lead samples supplied by the sources 


Weight, p.p.m. 








Sample Bi Cu Tl Ag Sb Fe Mg Mo 

Pb.3 1 0.1 <0.3 0.1 <0.1 

Pb.4 0.1 0.1 <0.3 0.1 <0.1 

Pb.5 Su | al <1 

Pb.6* tr. V.w. tr tr. tr 
<0.5 


Pb.7 <1 <1 <1 


NOTE: tr., trace; v.w., very weak. 
*Qualitative spectrochemical analysis only. 


TABLE III 








Estimated change from pure element (°C) 











Slope,* ———— - ~ - — ———— 
Impurity °C/p.p.m. Pb.4 Pb.3 Pb.6 Pb.7 Pb.5 
Of liquidus In liquidus 

Cu —0.0023 —0.00023 —0.00023 tr. —<0.0023 — <0.0023 

‘Bi —0.00044 —0.00004 —0.00004 — —<0.00044 —<0.00044 
Ag —0.00092 —0.00009 —0.00009 v.w. —<0.00046 — <0.00092 
Sb —0.00067 — <0.00007  —<0.00007 — — _ 

Tl +0.000088 +<0.000026 +<0.000026 — + <0.00009 =~ 


Estimated liquidus 











depression <0.0003 <0.0003 — <0.0031 <0.0037 
Of solidus* In solidus 
Cu —0.0200 —0.0020 —0.0020 tr. — <0.0200 — <0.0200 
Bi —0.00066 —0.00007 —0.00007 —_— — <0.0007 — <0.0007 
Ag —(0.0230 —0.0023 —0.0023 vw. — <0.0115 — <0.0230 
Sb —0.0021 — <0.00021 — <0.0021 — — — 
Estimated solidus 
depression <0.0046 <0.0046 — <0.0322 <0.0437 


Estimated melting range <0.0043 


<0.0043 = <0.03 <0.04 


*Phase diagrams, Metals Handbook (1948); Pb-Tl, Hansen (1958). These estimates assume that the liquidus 
and solidus are linear with composition: the solidus for Pb-Cu and Pb-Ag are more likely to be convex (possibly 
retrograde) with respect to the base of the diagram. 
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EXPERIMENTAL 


The apparatus and techniques used in these thermal analyses have been 
described previously* (McLaren, Parts I-V, 1957a, 6; 1958a, 6; McLaren and 
Murdock 1960). Samples of about 1750 g contained in high purity graphite 
crucibles under dry nitrogen atmospheres were used. This investigation 
includes an examination of the alloy melting ranges of the samples following 
different types of freezing with and without overnight anneals near the solidus 
temperature and a precision intercomparison of their liquidus points. 


A. Thermal Survey 
Figure 1 shows typical cooling curves that were obtained when a sample of 
high purity lead (Pb.3) solidifies under different rates and types of freezing. 
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Fic. 1. Typical freezing curves on high purity lead. 


The very fast normal freezes are made by extracting the melts from the furnace 
and allowing solidification to take place in room air; on slow induced and slow 
normal freezes the melts remain in the slowly cooling furnace block and on 
nucleation are, respectively, either assisted or not (by the induced freezing 
technique, see McLaren, Part II, 19576) to overcome a slow recalescence. A 
variable supercooling of up to 0.14° C was found on solidifying the high purity 
lead in this experimental arrangement. The plateaux of nearly constant 
temperature found on the cooling curves of slow induced freezes provide, as 


*Reference to preceding papers in this series is made by Part number and date of publication. 
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did similar plateaux on the induced freezing curves of high purity zinc, 

cadmium, and indium, an excellent index for a precision temperature standard. 

In order to assess the thermal effects of dissolved chemical impurities in the 

five samples of high purity lead, alloy melting ranges were gauged following 

different types of freezing and varying amounts of solid anneal near the 

melting temperatures. Table IV gives the approximate durations and trans- 
TABLE IV 


Transformation conditions for the thermal survey 


Melting 


Prior Solid 
Designa- freeze anneal, duration, 
Transformation Type tion Figure duration hr hr 
Melt After very fast normal freeze C1 2 4 min Z z 
Melt After very fast normal freeze C2 3 4 min 17 2 
with overnight anneal 
Melt After slow induced freeze £1 4 2 hr 2 2 
Melt After slow induced freeze ge 5 2hr 17 2 


with overnight anneal 


formation conditions for several series of heating curves from which estimates 
of the alloy melting ranges were made: samples were annealed at 320-+5° C. 

Two types of derived curves (see McLaren and Murdock, Part V, 1960) 
were determined from these melting curves for each melting: type a, inverse- 
“> duration against a melting temperature interval and 
The type } 


rate histograms of 
type }, incremented % duration against minimum melting range. 
derived curves are shown in Figs. 2-5 inclusive. 
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the total melting time spent in given melting ranges following slow induced freezes with 
2-hr solid anneals: derived from ‘g;’ curves. Fic. 5. Fraction of the total melting time spent 
in give n melting ranges following slow induced freezes with 17-hr solid anneals: derived from 
‘go' curves. 
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From the type a histograms the temperature widths at the 1% level of the 
total melting time were determined both for the over-all sample and, where 
found, the width at this level of any prominent core revealed on the melting 
curve; from the type 6 plots melting ranges for both 50% and 90% durations 
were determined. Table V summarizes the values of these melting range 
parameters that were made on all the samples. Consideration of these alloy 
melting range measurements will be deferred until after presentation of the 


TABLE V 


Alloy melting range comparisons 





Transformation 








type Ingot Pb.4 Pb.3 Pb.6 Pb.7 Pb.5 
Melting width (°C) at 1% level from inverse curves 
C1 w* 0.0019 0.0019 0.0029 0.0067 0.0069 
W 0.0019 0.0019 0.0040 0.0067 0.0175 
Ce w 0.0012 0.0029 0.0042 0.0064 0.0072 
W 0.0043 0.0029 0.0042 0.0064 0.0140 
£1 w 0.0020 0.0014 0.0008 0.0036 0.0083 
W 0.0028 0.0050 0.0047 0.0108 0.0083 
£2 w 0.0028 0.0019 0.0006 0.0061 0.0089 
W 0.0028 0.0046 0.0039 0.0061 0.0089 
Melting range (°C) for 50% duration 
C1 0.0004 0.0004 0.0007 0.0020 0.0072 
C2 0.0004 0.0004 0.0009 0.0025 0.0064 
£1 0.0020 0.0004 0.0004 0.0022 0.0039 
ge 0.0014 0.0004 0.0004 0.0050 0.0014 
Melting range (°C) for 90% duration 
Cy 0.0036 0.0036 0.0060 0.012 0.030 
Cs 0.0050 0.0033 0.0028 0.010 0.024 
£1 0.0028 0.0053 0.0047 0.013 0.022 
£2 0.0026 0.0097 0.0033 C011 0.018 
Table III (estimate) <0.0043 <0.0043 a <0.03 <0.04 





Note: w., grain cores; W., over-all sample. 


results of the liquidus point intercomparison because it is not possible to decide 
which of the samples of high purity lead is the purest solely on the basis of 
narrow alloy melting range. 


B. Intercomparison of Liquidus Points 

(a) Pressure Effects on the Freezing Temperature 

The effects of the external nitrogen pressure over the melts and the immer- 
sion pressure in the liquid metal on the plateau freezing temperature were 
experimentally determined in order that accurate corrections could be applied 
in temperature calculations at the lead point: these results are shown in Fig. 6. 
The increase in the freezing temperature of lead for 1 atm gas pressure over 
the melt was 0.0080° C. An immersion of about 10 cm (of the mid-point of the 
sensing element of a Meyers thermometer) was required before the measured 
temperatures agreed within 0.0001°C of the variation deduced from the 
Clausius—Clapeyron relation using the above measured gross pressure effect and 
a liquid density at the freezing point of 10.6 g/cc (Liquid Metals Handbook 
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Fic. 6. Pressure effects on the plateau freezing temperature of lead. 


1952). This expected variation of 0.0000819° C per centimeter is shown as the 
solid line in curve I of Fig. 6. 

(b) Freezing Temperature Comparison 

A careful intercomparison of the plateau freezing temperatures obtained on 
slow (~3 hours) induced freezes on the five samples of lead was carried out: 
the temperature of four samples were compared to that of sample Pb.4 using 
the measuring scheme that is regularly employed for precision temperature 
intercomparisons (see McLaren, Part III, 1958a). The fully corrected zero 
current resistances at the liquidus points of a single Meyers thermometer, 
S155, at 15 cm mid-point immersion corrected to 1 atm dry nitrogen pressure 








TABLE VI 
Intercomparison of plateau temperatures 
Plateau Liquidus 
Ratio Rpp/Rrp Reproducibility temp. depression 
~ - — rel. to from 
Mean No. of 10°X (equivalent 10-*°C) Pb.4, Table III, 
“G 


date Sample freezes Mean _ std. dev. Std. dev. Spread "€ 





14.9.59 Pb.3 10 2.24182546 +36 +1.2 3.9 +0.00012 <0.0003 
Pb.4 10 2504 +26 +0.9 2.7 0.0 <0.0003 

20.9.59 Pb.6 10 2.24181749 +30 +0.8 2.9 —0.0021 _ 
Pb.4 10 2503 +47 +1.3 5.0 

25.9.59 Pb. 7 5 2.24181402 +68 +1.9 4.6  —0.0031 <0.0031 
Pb.4 5 2521 +18 +0.5 1.1 

30.9.59 Pb.5 5 2.24181436 +56 +1.6 3.6 —0.0030 <0.0037 
Pb.4 5 2531 +24 +0.7 1.6 

(All freezes) Pb. 30 =.2.24182511 +36 +1.0 5.0 
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Fic. 7. Intercomparison of plateau freezing temperatures of Pb.3 and Pb.6 with Pb.4: 
Rrp and Rp», zero current fully corrected resistances of the thermometer at the triple point 
of water and the lead point respectively. 


over the melts were determined. Measurements at the triple point of water 
(triple cell 29) were made before and after the measurements at the freezing 
points of the lead samples and the average Ryp was used to determine Rp,/Rrp 
ratios. Figure 7 gives the results of comparing samples Pb.3 and Pb.6 against 
Pb.4, 10 determinations being made on each sample. Table VI summarizes 
the measurements that were made on all the samples. 

(c) Freezing Temperature of Lead 

A value of 327.4258° C (Int. 1948) was determined for the standard liquidus 
point of sample Pb.4. 

DISCUSSION 

(a) Purest Sample 

An inspection of Figs. 2-5 together with a consideration of the values of the 
melting range parameters in Table V afford the means to appraise differences 
in the effects of dissolved impurity solutes among the samples. These results 
show that samples Pb.4, Pb.3, and Pb.6 are clearly purer than the remaining 
two samples because they melt over distinctly narrower ranges following 
different types of freezing with and without the overnight solid anneals. On 
the other hand it is difficult to distinguish among the three purest samples on 
the basis of melting range alone. Evidence of strong coring can be seen on 
melting curves of all these samples which is variously dependent on prior 
thermal history, e.g. Pb.4, c2; Pb.3, gi, go; Pb.6, g1, ge. This coring is the result 
of considerable impurity segregation on freezing and shows that all the samples 
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contain residual impurity solutes in sufficient concentration to have associated 
alloy melting ranges that can be readily detected by the resistance 
thermometer. 

Further evidence of significant impurity segregation in all these samples may 
be deduced from Fig. 8, which shows a selection of the front ends of melting 
curves determined on the samples. Similar structures were found on the 
melting curves of samples of high purity tin (McLaren and Murdock, Part V, 
1960) and were thought to be associated with gross impurity segregations in 
the ingots because they were not observed on the melting curves of the purest 
samples. 

In the absence of a purer sample of lead that could be shown to have a 
narrower transition range that would be invariant with rate of freezing and 
amount of solid anneal, it is only possible to estimate which of these three 
samples is most likely to be the purest one. The liquidus points of Pb.4 and 
Pb.3 are nearly identical and were found to be about 0.002° C higher than that 
of Pb.6. The problem thus reduces to deciding whether the liquidus points of 
Pb.4 and Pb.3 lie closer to that of pure element than does the liquidus point 
of Pb.6. 

From Tables II and III it may be inferred that the three samples would melt 
in a temperature range bounded by the Pb-TI liquidus and the Pb—Ag solidus. 
The following factors weigh in favor of samples Pb.4 and Pb.3 having liquidus 
points closer than that of Pb.6 to the liquidus point of pure element. 

(i) From Table V it appears that sample Pb.4 generally had narrower over- 
all melting ranges than Pb.6. 

(ii) Reasonable agreement was found between the measured liquidus depres- 
sions (Table VI) and alloy melting ranges (Table V) and the values of these 
quantities that were estimated from known phase diagrams and the supplied 
chemical analyses. 

(iii) Sample Pb.6 showed a greater amount of coring (W/w..... £1, £2, 
Table V) following slow freezing: this must be attributed to a segregation of 
slightly solid-soluble solutes such as Cu and Ag. This observation is sub- 
stantiated by the qualitative spectrochemical analysis which indicates the 
presence of a trace of copper and silver in this sample. It is also significant that 
only 1 p.p.m. Cu would depress the liquidus about 0.002° C. 

(iv) For the liquidus points of Pb.4 and Pb.3 to be raised by 0.002° C from 
that of pure element would require the presence of 23 p.p.m. Tl, an amount 
nearly 100 times greater than their analyzed thallium contents. 


(b) Thermometric Standard 

The reproducibilities of the liquidus points of these samples of high purity 
lead are comparable with that of any fixed temperature point on the Inter- 
national Temperature Scale. Alloy melting range parameters may be used to 
select samples of lead that would have liquidus points within 0.001° C of 
pure element: samples with an over-all width at the 1% level of the total 
melting time on inverse melting curves of < 0.002° C following fast and slow 
freezing would meet this requirement. 
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Fic. 8. Portions of melting curves showing anomalous front end structures: types (see 
Table IV) of melting curves are indicated. 

The cadmium point (321.032° C Int. 1948) has three slight advantages over 
the lead point as a fixed temperature standard: samples are already available 
that have liquidus points within 0.001° C of pure element (McLaren, Part IV, 
19585); a smaller supercooling (0.03° C for cadmium compared with up to 
0.14° C for lead) decreases the amount of initial dendritic freezing which 
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possibly could interfere with the smooth interface coring that is responsible 
for the plateau temperatures; and cadmium has a larger (1.7 times) latent heat 
of fusion per cubic centimeter, which assists in the ready attainment of steady 
plateaux of long durations on induced freezes in simple block furnaces. Inter- 
polation errors arising in the use of the lead point as a fixed point in the 
quadratic calibration of resistance thermometers would be similar to those 
found for the cadmium point and discussed by McLaren (1959). 


CONCLUSIONS 


This investigation has shown that the liquidus point of high purity lead 
(327.426° C Int. 1948) may be used as a reliable precision fixed-temperature 
point and has an experimental realization that is comparable in simplicity and 
ease of attainment to the freezing temperatures of the pure metals previously 
studied. The resistance thermometer has been shown to have the sensitivity 
required to discriminate among samples of nominal 99.9999% purity lead on 


the basis of narrow alloy melting range. 
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THE WHITE TIN — GREY TIN TRANSITION 
IN TIN-MERCURY ALLOYS! 


R. W. SmitH? 


ABSTRACT 


The effect of the addition of mercury on the linear rate of white tin (@) — grey 
tin (a) transformation in a number of dilute tin-binary alloys has been examined. 
Rates of transformation at —78° and —25°C are shown for binary alloys of 
zone-refined tin and aluminum, arsenic, bismuth, germanium, indium, lead, 
antimony, and zinc, and also for a number of commercially pure tins. The resultant 
rates of all of these upon further alloying with 0.6 at.9 Hg are likewise given. At 
—78° C, the addition of mercury always produces a very low transformation 
rate, approximately 410-4 mm/hr. However, at —25° C, an increase in rate 
is observed in certain cases. 

The results are discussed and a theory is developed to account for them and 
also to predict the effect of mercury on the 8 — a transition in any tin-binary 


system. 


1. INTRODUCTION 


In a recent study of the a = 8 transition in tin-mercury alloys (Smith 
1959), it was shown that the introduction of mercury into pure tin yields 
transformation data valuable in understanding the general characteristics of 
this classic phase change. The rate of conversion of white tin into the grey 
allotrope was considerably reduced and by an amount which was virtually 
independent of the mercury content of the alloy. Also, the rates of trans- 
formation in alloys of zone-refined tin and mercury were found to be similar 
to that in one made from a less pure tin. A mechanism was proposed to 
for these observations based on the fact that mercury is mobile well 


below the transition temperature. The following experiments were devised 


to verify and extend this postulate. 
2. EXPERIMENTAL PROCEDURE AND RESULTS 
\ number of non-mercurial alloys were first prepared from zone-refined 
tin and certain metallic additions.* Then 0.6 at.% mercurial alloys were 


nade from one half of each. All were cold-rolled to a thickness of 1 mm, cut 
nto 0.1 K0.15%3.0 cm strips and annealed in argon at 200°C for 1 hour. 
Fach specimen was then fixed at one end toa metric scale and, after inoculation 
vith yrey tin at the other end, sealed in a glass tube to be maintained at the 
Test temperature 

When a grey tin wart had established itself across a specimen’s entire cross 
section, the further progress of the phase boundary was determined by accu- 
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rately checking its position against the scale at regular intervals. It was impor- 
tant that the position of a particular interface be determined in a few seconds 
to prevent change in the specimen’s temperature. A graphical plot of interface 
position against time gave an accurate estimate of the rate of transformation. 
This method of rate determination had been found to be the most convenient 
of a number tried (Smith, to be published). 

Two temperatures were employed, — 25° C, using a refrigerator, and —78° C, 
obtained by surrounding the specimen tubes with solid carbon dioxide. Four 
specimens of each material were tested at both temperatures, except in the 
case of tin-germanium, where the confusion incident upon spontaneous 
nucleation required more. 

The transformation rates for zone-refined tin, three commercially pure tins, 
and a number of alloys are shown in Table I. No values are given for the 
tin-lead alloy at —25°C since its linear transformation was immeasurably 
small over the 60-day period of the test. 

3. DISCUSSION 

One of the most significant points to note from the results is that the intro- 
duction of mercury into tin alloys will always produce a reduction in the rate 
of 8—a transformation when maintained at —78° C, whilst exposure at 
— 25° C may in some cases result in a rate increase. Also, the reductions in rate 
at —78° C are usually drastic and all to approximately the same final value, 
4X10-! mm/hr. Such a 8— a transformation rate is only brought about in 
pure tin by the addition of a relatively large amount of a potent transformation 
inhibitor, e.g. lead. Explanation of these observations must be associated 
with the melting point of mercury, —38.8° C. 

The 8— a transformation of a tin-mercury alloy requires the decomposition 
of any existing precipitates of HgSn,.* and the rejection of mercury at the 
a~—B phase boundary.t At —78° C, diffusion is slow and rejected mercury is 
solid. This results in considerable hindrance to the movement of an a-3 
interface and so impedes transformation, predominantly intergranular at this 
temperature (Bornemann and Gela 1953; Smith, to be published). 

It is not so directly obvious why alloying with mercury resulted in widely 
different transformation rates at —25° C, only producing an increase over the 
original value in certain cases. EFlowever, from an examination of the tollowing 
Cases, It is proposed to supply the reason. 

é 25° 


(a) Transformation of White Ttn Containing Only Mercury, at —25° C 

As a small region of grey tin forms in the white tin matrix, mercury is 
rejected at the interface (Ewald and Tufte L958, 1959). Lf possible, this 
mercury will combine with adjacent white tin, presumably to form ElgSnjy 
Further transformation will require this compound to decompose, again 
liberating mercury. Thus in the steady state, a cross section through a 3— a 


interface would reveal the white tin matrix, possibly with tree regions of ElgSnye, 


*HySnie is a hexagonal compound formed in white tin when mercury is present in excess 
of 0.2 at.{¢ (Groen 1954) 
tMercury is insoluble in grey tin (Ewald and Putte 1958, 1959) 
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a layer of HgSni2, then a layer of mercury, and finally the grey tin matrix. 
Since the mercury layer can only have a limited width, pockets of free mercury 
will also appear in the grey tin. Transformation can be pictured as the trans- 
ference of tin atoms from the 6-matrix, through the layers of HgSnj. and 
mercury, to the grey tin lattice for depositions. Such a transformation path is 
energetically more difficult than direct diffusion from 8 — @ and soa reduction 
in transformation rate is therefore to be expected when mercury is alloyed with 
pure tin. 


(b) Transformation of White Tin When Soluble Additions Are Present, Together 
with Mercury, at —25° C 

During the very early stages of transformation, the conditions described in 
(a) would exist at an a—8 interface. Continued transformation would result in 
the mercury layer becoming saturated with the other solutes originally present 
in the white tin. (Most elements have only a very limited solubility in mercury 
at —25°C and are insoluble in grey tin.) Thus a solute-rich layer would 
gradually form between the HgSnj. and the mercury. Accordingly, diffusion 
of tin atoms across the boundary region is further impeded and a lower trans- 
formation rate is observed than in (a). 


(c) Transformation of White Tin Containing Mercury and Also Other Solutes 
Present as Local Precipitates, at —25° C 

Initial transformation will follow (6) until the mercury layer becomes 
saturated and the localized precipitates in the white tin are no longer attacked. 
These particles will present local barriers to diffusion but, between them, the 
B-matrix will gradually dissolve away to be deposited as grey tin. This process 
will eventually result in the particles becoming completely detached from the 
B-surface. Hence only a slightly smaller transformation rate is to be expected 
when compared with that of pure tin and mercury alone. 


(d) Transformation of White Tin Containing Mercury and Grey-Tin-Soluble 
Additions, at —25° C 

Grey tin would only be expected to dissolve an appreciable amount of those 
elements which possess its own diamond cubic structure, e.g. germanium. 
However, a diamond cubic element would have only low solubility in tetragonal 
tin and thus it would be energetically favorable for any such solute appearing 
in the HgSnj,.—mercury interface to be taken up by the mercury and deposited 
on the grey tin. Such practice would tend to an increase in the rate of trans- 
formation over (c). 

These cases may be summarized thus: The solute properties tending to make 
the 8 — a transformation rate of a ternary tin—mercury alloy approach that of 
zone-refined tin and mercury are a localized dispersion in white tin, high 
solubility in mercury at —25° C, and a similar structure to that of grey tin. 

As yet only a prediction has been made of the transformation rate at 
—25°C of a ternary tin-mercury alloy. How this value differs from that 
observed in the binary alloy in the absence of mercury will depend on the 
particular solute and also on its concentration. Any decrease results from the 
reduced ease with which tin atoms diffuse through the mercury-rich phase 
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boundary region. An increase in transformation rate is effected through 
mercurial rearrangement of the solute which originally caused so much 
hindrance to the advance of an a—§ interface. 

The data shown in Table I agrees with the predictions. The rates of trans- 
formation of the ‘pure’ tins reflect the purity of each. Of the metallic aaditions 
to the zone-refined tin, only antimony, bismuth, and indium have any appre- 
ciable solid solubility in white tin at room temperature (Hansen 1958). That 
of zinc is only 2 at.% at the eutectic temperature and is probably very small 
at room temperature. However, the expected tendency for a very low transfor- 
mation rate at —25°C in the tin-indium-—mercury alloy, is more than offset 
by indium’s high solubility in mercury, 68.9 at.©@ at 10° C. The other metallic 
additions have very small or negligible solubilities in mercury below room 
temperature. The rather high rate of transformation in the alloy containing 
mercury and germanium is presumably due to the latter’s solubility in grey 
tin and that it may also act as a deposition center. 

Thus it is seen that the results of this study are in close agreement with the 


theory developed herein. 
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ON THE GEOMETRY OF RADIO REFLECTIONS FROM AURORA! 
P. A. ForsytTH 


ABSTRACT 


By assuming that auroral radio reflections are produced by volume scattering 
in clouds of ionization having the same spatial configuration as the visible auroral 
structures, and by taking into account the radar pulse duration and antenna 
beam width, it is possible to predict the probability of echo occurrence as a func- 
tion of range and azimuth. This echo distribution is quite similar to that observed 
experimentally even when ‘aspect sensitivity’? of the individual scatterer is 
neglected. Unfortunately, the optical evidence is not sufficiently extensive to 
permit precise calculations to be made, nor the radio evidence to permit detailed 
comparisons, but previous estimates of the shape of the scattering structures 
that have been based on the azimuthal echo distribution without regard to the 
factors discussed here are likely to be seriously in error. 


INTRODUCTION 


Although the reflection of radio waves from auroral ionization has been 
observed regularly for more than 10 years, there is no general agreement 
concerning the nature of the reflection mechanism. Soon after the first observa- 
tions it was noted that for stations in north middle latitudes radar echoes 
could be obtained only from auroral displays that occurred several hundred 
kilometers to the north of the radar station, that is, from auroral displays that 
were located at low elevation angles as seen visually from the radar station. 
This characteristic of the radar reflections still awaits a completely satisfactory 
explanation. The first attempt to provide such an explanation was made by 
Herlofson (1947), who suggested that the reflection process was specular and 
that reflections would occur only when the direction of incidence of the radio 
waves was accurately normal to the auroral structures, which, in turn, are 
known to lie along the earth’s magnetic field lines. At first, this explanation 
was accepted enthusiastically by most researchers, but as more experimental 
evidence was accumulated it became necessary to reject the specular reflection 
process in favor of a less rigid geometric control that has come to be known as 
‘faspect sensitivity’. Booker (1956) suggested that the observed aspect 
sensitivity could result from the diffraction pattern produced by long thin 
ionic inhomogeneities aligned with the magnetic field. Kaiser (1956) in England 
and Forsyth and Vogan (1957) in Canada concluded independently that the 
observed degree of aspect sensitivity could result from simple scattering by 
prolate spheroids having a ratio of length to diameter of about three or four 
to one. In this work the usual experimental evidence is the distribution in 
range and azimuth of the density of occurrence of radar echoes as seen from a 
single radar station. For stations such as Saskatoon (52.1° north, 106.6° west) 
the echo occurrence shows a strong tendency to peak in the approximate 
direction of magnetic north and at distances of 600 to 800 kilometers. This 

1Manuscript received December 7, 1959. 
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paper grew out of an attempt to arrive at a more accurate estimate of the 
degree of aspect sensitivity involved in auroral radio reflections by taking into 
consideration some of the known characteristics of visible aurora. As the work 
progressed it became apparent that echo distributions approximating closely 
to those observed experimentally could be derived without recourse to the 
“‘aspect sensitive’ scattering structures and that the statistical characteristics 
of visible aurora are not known at the present time with sufficient accuracy to 
arrive at a meaningful estimate of the shape of the scattering structures by 
this method. Even though these echo distributions can be calculated only 
crudely, it is felt that it is worth while to present them at this time because 
they depend critically on factors that have been almost completely ignored in 
previous attempts at quantitative analysis of the observed distributions. 


CALCULATION OF ECHO PROBABILITY DISTRIBUTIONS 

There is still some uncertainty concerning the extent to which the radio 
reflections can be associated with visible auroral displays, but there seems to 
be no doubt that at least a significant fraction of the radar echoes observed at, 
say, Saskatoon are produced by ionization associated with auroral structures 
that could, under the proper conditions, be observed visually. What evidence 
there is on the subject suggests that the auroral echoes tend to come from the 
more well-defined auroral forms. One feature of these forms is a characteristic 
configuration in space. They frequently extend for many hundreds and perhaps 
even thousands of kilometers in the east-west direction (St6érmer 1936), some 
30 or 40 kilometers in the vertical direction (Currie and Weaver 1955), and 
only about a kilometer or perhaps less in the north-south direction. The first 
basic assumption used in the computations is that the echo-producing auroral 
ionization has the same configuration as the auroral luminosity. The second 
assumption is that the reflection process is a non-coherent, volume effect. 
Another way of stating this assumption is to say that many scatterers are 
involved, and at least some of the scatterers are separated by several wave- 
lengths. No specific assumption is made about the size or shape of the 
scatterers. 

A wide variety of radar parameters have been used to investigate aurora. 
-Pulse durations varying from about ten to several hundred microseconds have 
been used. Aside from the usual expectation that the radar sensitivity in- 
creases with increasing pulse duration, relatively little attention has been paid 
to this important parameter in making comparisons between data obtained 
with different radars. If, however, the auroral ionization has the spatial 
configuration outlined above, then the distribution of echoes in range and 
azimuth should depend markedly on both the pulse duration and antenna 
radiation beam width of the radar. Two situations are illustrated in Fig. 1, 
which shows a plan view of the radar location in relation to two auroral forms. 
In Fig. 1(@) the auroral structure is located well to the north of the radar and 
the direction of travel of the incident radio waves is nearly normal to the length 
of the auroral form. If the range increment dR corresponds to the pulse dura- 
tion 7, of the radar (i.e. dR = cT/2) then the volume of auroral ionization 
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Fic. 1. Plan view of area surrounding a radar station, illustrating the geometric factors 
that enter the calculation. 


that can simultaneously contribute to the echo is proportional to the length 
AB. In Fig. 1(6) the radio wave front makes an appreciable angle ¢ with the 
length of the auroral structure and the corresponding volume is proportional 
to the much shorter length EF. Actually, in order to determine the distribution 
in azimuth of the echoes it is necessary to use a directive antenna (i.e. an 
antenna of limited beam width). In this case, a further restriction on the 
simultaneous scattering volume may arise. In Fig. 1(@) is sketched an antenna 
beam width of 20° for which the scattering volume would be proportional to 
CD rather than to 1B. This simple approach has been followed in the calcula- 
tions on which this paper is based. 

The computation of the scattering volume as a function of range and azimuth 
was carried out for Saskatoon, initially for the assumption that all auroral 
forms were parallel and roughly aligned with the lines of geomagnetic latitude. 
The volume parameter x (corresponding to EF in Fig. 1(0)), was calculated for 
about sixty field points in a region representing one (the eastern) half of the 
radar coverage. For this calculation the following approximate expression was 
used, 


x = R (sin? ¢+2dR/R)'—R sin ¢, 


where the symbols have the meaning indicated in Fig. 1. The quantity x, 
according to the initial assumption, should be proportional to the power 
reflected by the auroral form. There are some conditions under which the 
probability of observing an echo may be proportional to this quantity raised 
to some fractional power, but for the present purpose, it seemed wisest to 
assume the more usual simple proportionality between echo power and echo 
detectability. In general, there are two segments, x, which contribute to the 
scattering volume. Whenever the centers of these two segments lay within an 
azimuth sector of 20° the larger figure, 2x, was used to compute the echo 
probability. This procedure is equivalent to assuming an antenna beam width 
of 20°. Although the use of a specific value detracts from the generality of the 
treatment, some antenna beam width must be assumed and 20° seems fairly 
representative of those that have been used. This value was used also through- 
out the computation to restrict the volume in the manner indicated in Fig. l(a). 

Similarly, some range dependence must be assumed. The particular range 
dependence chosen corresponds to that frequently found for v.h.f. radars. 
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Fic. 2. Contour maps for the neighborhood of Saskatoon showing the calculated distri- 
butions of echo occurrence for a number of different radar pulse durations (7) and constant 
antenna beam width (6). The maps (a), (0), and (c) assume that all the auroral arcs are parallel 


while (d), (e), and (f) are smoothed for variability of arc orientation. 
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The radar sensitivity was assumed to be proportional to 0.04, 1.0, 1.0, and 0.2 
at ranges of 400, 600, 800, and 1000 kilometers respectively, and zero at ranges 
appreciably greater or less than these values. 

Finally, the distribution of density of occurrence of visual aurora must be 
taken into account. This distribution must, of course, be specific to the location 
of the radar, in this case Saskatoon. The values used were those given by 
Vestine (1944). Once again, the use of these particular values might be 
questioned as being inaccurate, but no more reliable information is at hand. 
The probability of echo occurrence for each field point was calculated by 
multiplying together the factors representing the scattering volume, the radar 
sensitivity, and the occurrence of visible aurora. In Figs. 2(a), 2(0), and 2(c) are 
plotted contours of the probability of echo occurrence for the region surround- 
ing Saskatoon for three different pulse durations, 20, 100, and 500 micro- 
seconds corresponding to range increments of 3, 15, and 75 kilometers 
respectively. 

The assumption that all auroral structures are parallel provides an instruc- 
tive example, but is not in accord with the observations. Jensen and Currie 
(1953) found that at Saskatoon auroral arcs tended to have directions that 
varied through about +10° about the mean direction. In order to take account 
of this variability the value of @ for each field point was varied through +10° 
about the value used for the first computation and the average resulting value 
for x was used. The resulting echo probability distributions, smoothed in this 
way for variability in direction of the auroral forms, are plotted in Figs. 2(d), 
2(e), and 2(f). The same three values of pulse duration were used for these 


calculations. 
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DISCUSSION 

The most striking feature of the echo occurrence graphs of Fig. 2 is the very 
strong concentration in the northward direction. The concentration is even 
more marked than is usually encountered in experimental investigations. It 
should be noted that while the north direction in these graphs was referred to 
earlier as being coincident with the geomagnetic meridian, it should be taken 
more precisely as the normal to the mean local orientation of the auroral forms. 
For Saskatoon, Jensen and Currie (1953) found that this mean orientation 
changed seasonally through more than 25°, a factor that could contribute to 
the width of the observed distribution. 


Variation with Latitude 
One of the strongest arguments in favor of the existence of aspect sensitive 


scattering structures is based on radar observations made at College, Alaska. 
There, even though the location is close to the maximum of the auroral zone, 
radar echoes are obtained most frequently from aurora well to the north of the 
station. An approximate calculation of the kind carried out for Saskatoon was 
also carried out for College. The curvature of the lines of geomagnetic latitude 
in the neighborhood of College differs slightly from the curvature in the neigh- 
borhood of Saskatoon but for the calculation this slight difference was ignored 
and the values of x calculated for Saskatoon were used. College is not only 
farther north geomagnetically than Saskatoon, but in that area the auroral 
zone is farther south. The geomagnetic latitude of College was taken as 64° and 
that of the auroral zone as 67°. (The corresponding figures for Saskatoon were 
60.5° and 69°.) The calculation was carried out only for the 100-microsecond 
pulse duration with smoothing for variability of direction of auroral forms. The 
results are shown in the form of a contour diagram in Fig. 3, where it may be 
seen that there is still a strong concentration of predicted echo occurrence in 


the northern direction. 


Azimuthal Distribution 

Another way of studying the relative occurrence of echoes is to examine the 
variation of echo occurrence as a function of azimuth angle for constant range. 
Several curves derived in this way for a range of 700 kilometers are shown in 
Fig. 4. The results of the three calculations carried out for 7 = 100 micro- 
seconds; unsmoothed and smoothed for Saskatoon and smoothed for College 
are shown. Also shown is a histogram of the experimental distribution found 
at Saskatoon in an earlier investigation by Currie, Forsyth, and Vawter (1953). 
For that investigation, using a frequency of 56 Mc/s, the pulse duration was 
only 16 microseconds and the antenna beam width was about 45°, much wider 
than that used for the calculation. Both these factors should tend to sharpen 
the distribution in the northern direction, but it is seen that the observed 
distribution is much wider than the calculated distribution. 

The simplest way in which the computed distributions could be modified to 
agree more closely with the observed distribution is to increase the degree of 
smoothing, that is, to assume a greater variability in orientation of auroral 
structures. While the evidence of Currie and Weaver is probably valid for most 
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Fic. 3. Contour map of echo distribution in the neighborhood of College, Alaska, calcu- 
lated for the radar parameters of Fig. 2(e). 


extended auroral formations, large local deviations in orientation frequently 
occur, particularly for the more active forms of aurora. It may be that this is 
the reason that auroral radar echoes seem to be preferentially associated with 
the more active auroral forms. A recent study by Lyon (1960) has shown that 
for a radar beam directed at an angle of 45° from the geomagnetic meridian the 
echoes occur most frequently during the period of breakup of the visible aurora, 
the period in which the quiet arcs break up into bands and then into the rayed 
forms. In particular, there seems to be a tendency for the echoes to be produced 
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Fic. 4. Comparison of various distributions of echo occurrence with azimuth for constant 
range (700 km). 
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most frequently by the so-called ‘‘curl” form first discussed by Webster (1957) 
and later by Gartlein (1959). It might be noted that if the scattering volume 
effects discussed in this paper are important then the results of Lyon’s investi- 
gation are entirely reasonable. In Fig. 5, which again shows a plan view of the 
area around the radar station, this time including two auroral curl forms, it may 
be seen that for certain combinations of pulse duration and antenna beam 
width the azimuthal distribution of echo occurrence could be modified by the 
appearance of curl forms. The probability of echo occurrence for a form at the 
azimuth of (6) would be increased relative to that of the simple form in Fig. 1, 
resulting in an increase in the width of the distribution. At the same time it 
may be noticed by comparing Figs. 5(@) and 5(0) that this effect will only 


(a) % be / ~ 
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Fic. 5. Plan view of area surrounding a radar station, illustrating the effect of the occurrence 
of ‘‘curl” forms on the echo distribution. 


modify the width of the peak in the distribution. As long as the extent of the 
auroral form is substantially greater in the east-west direction than it is in the 
north-south direction, then the northern peak in echo occurrence will persist, 
even for stations located under the auroral zone. It would be necessary to 
locate a radar well to the north of the auroral zone in order to cause this 


northern peak to disappear. 


Frequency Dependence 

The experimental evidence concerning the frequency dependence of the 
azimuthal distribution of echoes is confusing. In the earlier work at Saskatoon 
using 56 Mc/s and 106 Mc/s the higher frequency radar had the longer pulse 
duration and the narrower beam width. The resulting azimuthal distribution 
was slightly more peaked for the higher frequency. A recent investigation by 
Presnell et al. (1959) at College used frequencies of 216 Mc/s and 398 Mc/s. In 
this case various pulse durations were used at both frequencies and the 
antenna beam width was inversely proportional to frequency. The peak in the 
resulting azimuthal distribution is slightly sharper for the lower frequency. It 
is unfortunate that there seems to have been no investigation of the azimuthal 
distribution of echoes carried out at two or more frequencies with radars having 
equal pulse durations and accurately scaled antennas. On the other hand, there 
is ample evidence to show conclusively that the occurrence of radar echoes falls 
off rapidly with increasing frequency. It is reasonable to infer from this that 
as the radar frequency is increased, the number of effective scattering centers is 
reduced so that the scatterers become more widely separated in space as well 
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as in time. If this be the case, a larger volume of the auroral structure must 
contribute in order to produce a threshold echo at the higher frequencies. Thus 
as the frequency is increased the importance of local variations in orientation 
would decrease and the resulting azimuthal distribution of echoes would 
approximate more closely to the unsmoothed distributions of Figs. 2 and 4. 

Since the treatment has made use of parameters typically associated with 
v.h.f. radars, some modifications would be necessary in a consideration of the 
distribution to be expected for a u.h.f. radar. In addition, the most recent 
results obtained with u.h.f. radars operating at frequencies between 400 and 
500 Mc/s suggest that the assumptions of this paper may not apply to all 
echoes observed with these radars. Many echoes display the characteristics to 
be expected for reflection from ionization having the same spatial configura- 
tion as the visible aurora, but others appear to come from patches of ionization 
which occur within a relatively thin layer. The importance of this second type 
of echo relative to the first is not yet established but the present treatment 
could not apply directly to reflection from a horizontal layer. It may be that 
both of these kinds of ionic volume distribution contribute also to v.h.f. 
radar reflections. Some of the earlier observations could be interpreted in this 
way. The fact that the two distributions have not been clearly distinguished 
before might be due simply to the poor angular resolution that is characteristic 
of radars operating in the lower frequency range. 


Vertical Extent of urora 

The calculations and discussion up to this point have left out of consideration 
the effect of the vertical extent of the auroral forms. It seems likely that the 
vertical extent of the echo-producing volume is somewhat less than that of the 
corresponding visible form. Nevertheless, it may be instructive to consider the 
effect of the vertical extent of the radio-aurora, assuming it to be identical with 
that of the visible aurora. Clearly, the two-dimensional computations can be 
extended to three dimensions. For short pulse durations the echo amplitude 
(and hence echo occurrence) will fall off for the higher elevation angles. The 
elevation angle dependence that arises in this way will be much reduced for 
pulse durations corresponding to range increments of more than the vertical 
extent of the forms, say for pulse durations greater than about 300 micro- 


seconds. 


Variation with Elevation Angle 

There is another variation with elevation angle to be expected that is 
independent of the vertical extent of the auroral forms. This is due to the fact 
that the lower edge of an auroral form is very nearly horizontal (located at a 
constant height near 100 kilometers). As the elevation angle increases, the 
range to the 100 kilometer level decreases and the curvature of the wave front 
at that level increases. As a result of this increasing curvature the scattering 
volume of an auroral form, as seen by a radar of fixed pulse duration looking 
toward the north, decreases with increasing elevation angle. Other factors such 
as antenna beam width and auroral absorption probably affect, and may even 
dominate, the distribution of auroral echoes in elevation angle. Here again, few 
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experimental data exist, largely because of the difficulty in producing an 
antenna beam of small vertical dimension that is variable in elevation angle. 


Continuous-Wave Experiments 

Quantitative experiments using continuous-wave, as opposed to pulsed, 
radar systems are comparatively rare. If the approach followed in the fore- 
going discussion is pertinent to the natural situation, then there should be 
quite striking differences between the two types of system. It would be expected 
that the continuous-wave systems would experience a large advantage in 
sensitivity over the pulsed system, in addition to the usual band width 
advantage. The few experimental results now available, such as those of 
Forsyth and Vogan (1957), Collins (1958), and Collins and Forsyth (1959) 
tend to confirm this expectation. Some of these results were obtained with 
bistatic continuous-wave systems in which the transmitters were separated 
from the receivers by a distance of about 1000 kilometers. For bistatic systems 
the geometric considerations become more complex, but it is worth noting 
here that for these systems also there should be major differences between 
pulsed and continuous-wave systems. The pulsed system will have a sharp 
maximum of sensitivity for auroral structures located near the mid-point of 
the direct path from transmitter to receiver, and additional maxima for 
structures located well to either side of, and parallel to, the direct path. In 
contrast, the continuous-wave system will have a much more uniform sensi- 
tivity for structures located in various parts of the area of coverage. 


CONCLUDING REMARKS 

It is evident that if it is assumed that auroral radar echoes arise 

(a) from configurations of ionization that are similar to the visually 
observable auroral structures, and, 

(6) by means of volume scattering, 
then many of the observed characteristics of the echoes may be explained 
simply. The author is well aware that the experimental support for this 
approach is far from conclusive. Many of the same characteristics may be 
explained just as convincingly, though perhaps not as simply, by various com- 
_binations of aspect-sensitive scattering centers. The existence of such scattering 
structures was first proposed to explain the strong northern peak in the 
azimuthal distribution of auroral radar echoes. It seems hardly likely that this 
peak is due exclusively to the combination of effects discussed in this paper, 
but it is very unlikely that they play no part at all in its formation. If this point 
of view is accepted then it is not possible to determine the aspect sensitivity of 
the individual scattering structures by study only of the azimuthal echo 
distribution. 

Some method of separating the effects of aspect sensitivity from those of the 
volume scattering process discussed in this paper is clearly required. At first 
glance the small local maximum in the southern direction (see Fig. 2) for 
volume scattering might seem suitable. Unfortunately, this maximum tends to 
disappear with smoothing of the kind used to derive the distribution of Figs. 
2(d), 2(e), and 2(f) (see also Fig. 4). For this reason, it is not likely to be easily 
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observed. On the other hand, for the particular case of a radar with a long 
pulse duration (say, more than 1 millisecond) and a very narrow beam width 
(one or two degrees) there might occur a local maximum in the azimuthal echo 
distribution for a beam direction parallel to the orientation of the auroral 
forms. The discovery of such a maximum would confirm the assumptions of 
this paper, but it would still be difficult to arrive at a quantitative evaluation 
of the aspect sensitivity of the scattering structures. Such an evaluation 
must await some more direct measurement, possibly a radar system capable 
of resolving an individual scattering center would also provide an indication 
of its shape. 
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SPIN LATTICE RELAXATION IN NEODYMIUM ETHYLSULPHATE 
AT LIQUID HELIUM TEMPERATURES! 


J. M. DANIELS AND K. E. RIECKHOFF 


ABSTRACT 


The optical Faraday effect was used to measure instantaneous magnetization in 
neodymium ethy!sulphate. The spin populations were disturbed by pulses of 
microwave power, and by adiabatic magnetization and demagnetization, and the 
approach to equilibrium was studied. The relaxation was found to be exponential 
and spin lattice relaxation times were measured, for temperatures between 1.3° K 
and 4.2° K, and for magnetic fields between 80 and 6000 oersted. The relaxation 
time was found to decrease with increasing magnetic field, and to vary with 
temperature approximately as 1/T%. No dependence of relaxation time on pulse 
length was found. 


1. INTRODUCTION 

The Faraday effect (the rotation of the plane of polarization of light passing 
through a material in a magnetic field) is very large, in general, in paramagnetic 
salts. The rotation of the plane of polarization is proportional to the magnetic 
moment of the salt, and for a single crystal of neodymium ethylsulphate the 
maximum rotation corresponding to saturation of the magnetization is about 
114° per mm of path length for green mercury light (A = 5461 A) propagated 
along the optic axis (Becquerel ef a/. 1938). Thus, the Faraday rotation can be 
used to measure the spin temperature of a sample of one of these paramagnetic 
salts. It can also be used to detect changes in spin temperature (or relative 
populations of the spin states) under the influence of disturbing factors. For 
example, microwave radiation, of a frequency which produces paramagnetic 
resonance, alters the relative populations of the spin states, and produces a 
change in the Faraday rotation. This effect was predicted by Kastler (1951), a 
theory was given by Opechowski (1953), and the effect was verified experi- 
mentally for neodymium ethylsulphate by Daniels and Wesemeyer (1958) (see 
also Wesemeyer and Daniels (1958)). In the last publication, the authors used 
this effect to estimate the spin lattice relaxation time by the saturation 
method. The Faraday effect can, with suitable detection and presentation 
equipment, be used to measure instantaneous spin populations, and to follow 
the course of the spin populations as a function of time during the return to 
equilibrium from a disturbed state. It is thus a useful tool for a study of spin 
lattice relaxation, and this paper presents the results of such a study of neo- 
dymium ethylsulphate. Measurements were made in the temperature range 
1.3° K to 4.2° Kk, and for fields between 80 and 6000 oersted. 


2. EXPERIMENTAL PROCEDURE 
The apparatus used was originally a conventional X-band paramagnetic 
resonance spectrometer, suitable for use at liquid helium temperatures, with a 
1Manuscript received January 25, 1960. 
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transmission cavity and crystal video display. The microwave power, up to 
140 mw, was supplied by a 2K39 klystron. The following modifications were 
made to enable the optical studies to be carried out: 

(a) A magnet with a hole through the pole pieces parallel to the magnetic 
field was used. 

(6) Two holes about 2 mm in diameter were made in the broad face of the 
cavity, near the center of the top of each face, to permit light to pass through 
the cavity and the crystal. The crystal was mounted on a removable choke 
plunger at the top of the cavity, as is usual. 

(c) The cavity was immersed in liquid helium contained in a dewar vessel, 
and this was itself immersed in another dewar containing liquid nitrogen. The 
tails of these dewars were left unsilvered. A cylinder of copper foil was wrapped 
round the tail of the inner dewar and holes 4 mm X4 mm were cut in this to 
let the light pass through. Square windows of the same size were cut out of 
Scotch brand No. 33 electrical tape, and were stuck over the copper foil so 
that the inner dewar fitted tightly into the outer dewar. In this way, the path 
of the light was kept free of liquid nitrogen, and of nitrogen bubbles which 
would have introduced disturbing fluctuations into the transmitted light. 

(d) The input coaxial line was critically coupled to the cavity, to ensure 
maximum power transfer. The output coaxial line was weakly coupled, and 
the current of the silicon—-tungsten rectifier crystal was used only to monitor 
the tuning of the cavity and klystron. 

(e) The modulation coils of the magnet, which are usually energized with 
60-cvele alternating current to provide a held sweep, were not used for their 
normal function. 

The light source was a General Electric mercury arc, type H100-4A. It was 
energized from a d-c. generator via a filter which stabilized the current in the 
arc, and reduced the ripple. The light passed through a collimator, a glass 
filter which passed the 5461 A green mercury line, a Glan—Thompson polarizing 
prism, and a system of lenses to condense the light on to the sample in the 
cavity. The light which emerged from the cavity passed through another 
collimator, a Glan—-Thompson prism analyzer, and was incident on the photo- 
cathode of an RCA 5819 photomultiplier. This experimental system is depicted 
in Fig. 1. 

The photomultiplier voltages were supplied from batteries. The anode load 
was 10 kQ, shunted by a 0.01 uf condenser to cut down the noise bandwidth. A 
galvanometer was put in series with the 10-kQ resistor, and the voltage output 
from the anode circuit was connected to the Y,; d-c. amplifier of a Cossor 1049 
double beam oscilloscope. Time markers from a General-Radio 1217-A unit 
pulser were fed to the Y2 plates of the oscilloscope. The trace of the oscilloscope 
was photographed. A Tektronix model 502 oscilloscope was used in later 
experiments. Before an experiment, the electrodes of the photomultiplier were 
prefatigued by running the photomultiplier for some 10 hours with a weak 
illumination on the cathode. The linearity of the recording system was then 
checked using the cos? @ transmission law for polarized light. 

The neodymium ethylsulphate crystals were made from neodymium oxide 
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Fic. 1. Diagram of the apparatus. H, mercury arc; F, filter; C, cavity; X, crystal of 
neodymium ethylsulphate; D, dewar vessel; P, photomultiplier. For reasons of clarity, only 
one of the dewars is shown, and no attempt has been made to draw the details of the micro- 
wave cavity and its electrical connections. 


supplied by the Lindsay Chemical Company. It is described as 99% pure, the 
impurities being 0.1% to 0.4% of each of praseodymium and samarium, and 
up to 0.5% of other rare earth oxides. The oxide was first converted into the 
sulphate by dissolving it in sulphuric acid and crystallizing. Barium ethyl- 
sulphate, freshly prepared according to Beilstein (1918), was added to the 
calculated amount of neodymium sulphate both in solution. The solution was 
filtered to remove the precipitated barium sulphate, and crystals were grown 
in a vacuum desiccator. The neodymium ethylsulphate was kept in a refrigera- 
tor until required. Single crystals were grown by a variety of methods, but 
one of the best was to leave about 10 cc of a saturated solution covered and 
undisturbed in a room at a fairly constant temperature for a week or two. It 
was usually possible to cut out a portion of one of these crystals which was 
sufficiently good optically for the experiments. Two faces, perpendicular to the 
optic axis, were cut and polished on these crystals. 

When mounting a crystal in the cavity, care was taken that the hexagonal 
axis was parallel to the direction of propagation of light. This is necessary 
because the crystals are doubly refracting, and it was usually possible to adjust 
the orientation of the crystal so that the minimum intensity passed by the 
analyzer was about 1% of the maximum. This amount of ellipticity in the 
polarization of the emergent light can introduce errors of up to 2° in the 
measurement of the Faraday rotation. 

The temperature of the helium bath was maintained by pumping on the 
liquid helium through a throttle valve. The pressure above the liquid helium 
was measured using a mercury and an oil manometer, and was kept constant 
by manual adjustment of the throttle valve. The temperature of the bath was 
obtained from this measured pressure. 

A multivibrator, which delivered square pulses to the reflector of the klystron, 
was used to switch the klystron on and off. The pulse length could be varied 
between 8 and 124 mseconds. The multivibrator could either run free at about 
one cycle every 2 seconds, or could be manually triggered to give single pulses. 


The leading edge of the pulse was used to trigger the oscilloscope sweep. The 
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paramagnetic resonance line is so wide that at liquid helium temperatures 
100 mw of microwave power is sufficient to saturate the spin system at all 
fields between about 800 and 2000 oersted. 

The procedure during an experiment was more or less as follows. First the 
throttle valve on the helium pumping line was set, and the helium bath was 
allowed to come to equilibrium. Only minor adjustments of the throttle valve 
were needed subsequently to maintain this temperature. The alignment of the 
crystal was then checked. A curve of the Faraday rotation, as a function of 
magnetic field, was obtained by rotating the analyzer to a position of minimum 
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Fic. 2. Typical oscilloscope traces showing the saturation due to the microwave pulse, 
and the subsequent relaxation. The time markers are 1/30 sec apart. The irregularities on the 
“> 


trace are 60 cycle pickup from the stray light in the room. Figure 26 shows a trace in which 
the rotation is so large that the transmission through the analyzer prism passes through a 
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minimum, 
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photomultiplier current for various values of magnetic field. The rotation was 
read off from a scale of degrees on the analyzer housing. The klystron was then 
switched on, and the multivibrator adjusted so that the pulses brought the 
klystron into tune with the cavity. Next, a value of the magnetic field was 
selected, and the analyzer was adjusted to some convenient position, usually 
15° from the position of minimum transmission. Finally pulses of microwave 
power were fed into the cavity, and photographs of the oscilloscope trace were 
taken. This procedure was repeated for different fields and temperatures. The 
photographs, which were on 35-mm film, were mounted on the stage of a low 
power binocular microscope, with millimeter graph paper behind, and the 
trace was measured by visual inspection. Figure 2 shows some typical traces. 

Another series of experiments was carried out without using microwave 
power at all. It was found that, by connecting the modulation coils of the 
electromagnet to a 110-v d-c. supply, the field could be changed by some 250 
gauss, coming to equilibrium in about 10 mseconds. When the field is thus 
suddenly changed, the spin temperature is changed also, a process of adiabatic 
magnetization or demagnetization. This technique has the advantage that it 
can be used at any field, and is not confined to the region 800 to 2200 oersted 
in which microwaves are absorbed. It has the disadvantage, however, that the 
change in population, and hence rotation, which is produced is much less than 
can be produced by microwave absorption near resonance. 


3. RESULTS 
According to the phenomenological theory of relaxation, the difference in 
the populations of the two spin states of the Nd*+** ion, v, follows an equation 


© og ee ae 


(1) 

dt T 
7 is the spin lattice relaxation time, v» is the equilibrium value of v, determined 
by the bath temperature, and W is a factor which is proportional both to the 
absorption coefficient of the crystal, and to the microwave energy in the cavity. 
In the presence of some radiation in the cavity, a steady state value 1 is 
reached where 

Vy = Vo (1 + Wr), 


whence equation (1) becomes 


dv _ (v—v1)¥% 


(2) Be aie 
dt ViT 


Hence v relaxes exponentially to » with a time constant 7y,/vo. Since the 
Faraday rotation @ is proportional to v, we can expect @ to relax to a value 6;, 
somewhat less than the rotation @) measured with no microwave power in the 
cavity. This relaxation should be exponential with a time constant 76/4. 
Experimentally, it was dithcult to adjust the multivibrator so that the klystron 
was switched off completely during the off-period; there was usually a small 
radiation background in the cavity, hence this correction for background 


radiation is important. 
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The Q of the cavity was about 10,000, hence the time constant for the radia- 
tion in the cavity to come to equilibrium was about 10~ sec. This is so much 
less than 7 that radiation reaction effects are negligible, and the radiation in 
the cavity can be considered to be in equilibrium at all times. 

The Faraday rotations, 6, obtained from the photographs, were plotted as a 
function of time on semilog paper. In every case, straight lines were obtained 
covering up to 90% of the whole decay, and the slope of these lines gave 
70,/0). We estimate an accuracy of +10% in the final values of r. 

The first topic to be investigated was the effect of the length of the micro- 
wave pulse on the relaxation time. Speculations have been raised that in the 
principal mode of relaxation at low temperatures the ions give their energy to a 
relatively narrow band of lattice vibrations, which in turn gives its energy to 
the other lattice vibrations. Van Vleck (1941a, 5) suggested that the observed 
relaxation times might depend on the rate of transfer of energy between the 
lattice vibrations rather than on the rate of transfer of energy between the ions 
and the lattice vibrations, and one criticism of almost all theories of relaxation 
is the assumption of thermal equilibrium among the lattice vibrations. While 
microwave power is being absorbed by the crystal, energy is being fed into the 
lattice, and there is always the possibility that when the microwave power is 
switched off and the system begins to relax, the lattice is in a non-equilibrium 
state with a consequent change in the relaxation time. Measurements were 
made at different temperatures and fields for pulse lengths between 8 and 124 
mseconds. In no case was there found any dependence on pulse length. The 
result of a typical set of measurements is shown in Table I. Measurements 


TABLE I 


Measured relaxation times at 1.685° K and 
1540 gauss, for different pulse lengths (two 
measured relaxation times for each pulse length) 


Pulse length, Relaxation time, 


msec msec 
8 145 135 
17 136 162 
35 155 134 
63 125 166 
124 147 136 


were made of relaxation time as a function of field and temperature. These 
were possible between about 800 and 2500 gauss, outside this range the absorp- 
tion of microwaves was so small that it was not possible to disturb the popula- 
tion of the spin level sufficiently with the power available. Since an accuracy 
no better than 10% is claimed for these measurements, there is no advantage in 
trying to present them in graphical form, as the picture would be somewhat 
confused. Instead we present these results as a table (Table II) and we would 
like to point out that the values quoted are the experimental measurements. 
They are in no way smoothed, and are quoted to at least one figure more than 
is significant, in order to give an idea of the spread of the points. We see that 
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TABLE II 


Relaxation times as a function of temperature and magnetic field, measured by the micro- 
wave method. Times in mseconds. These are unsmoothed measured values of accuracy not 
more than +10% 











~—_ H, gauss 

TK —~_| 780 1030 1290 1540 1790 2040 2290 2540 
1.380 | 184 196 208 170 163 126 121 125 
1.542 | 129 134 154 144 141 121 118 116 
1.685 | 116 176 159 144 155 151 117 108 
1.875 | 57 67 75 68 67 68 64 65 
2.06 | — 69 72 73 67 65 63 61 
4.22 _ aa 11 11 13 16 aa iva 





the values of relaxation time decrease rapidly with temperature, and that they 
decrease gradually with increasing magnetic field, with a broad maximum in 
the neighborhood of 1500 gauss. 

The results of relaxation times measured by the technique of suddenly 
changing the applied magnetic field are shown in Fig. 3. These cover a much 


400 + © 1,338°K 
. 46 © 1.78 °K 
msec 
® 2.082 °K 
| aaah ie 
300 | > f° OooN\ 
im ° 
aie e ae 5 
200 | oa, 
e é an 
as 
> % 
\ 


100 |} a 


Ged eee ! Bo 
Oo | 2 3 4 5 KOe. 6 


Fic. 3. Relaxation time as a function of magnetic field. The few measurements at 1.78° K 
> 


and 2.082° K were done with crystal No. 2. 
larger range of field, but this method is inherently less accurate than the micro- 
wave method, since the relatively small change jn the field used, 250 gauss, 
produces a much smaller change in the Faraday rotation than can be obtained 
using microwaves. All the measurements made with the microwave method, 
and also the long series, using the change of applied field, at 1.338° K were 
made on the same crystal of neodymium ethylsulphate. The other results 
depicted in Fig. 3 were made from a different crystal, grown, however, from 
the same batch. 

Figure 4 is a series of log—log plots of relaxation time against temperature for 
four chosen fields, 1290 gauss, 1540 gauss, 1790 gauss, and 2040 gauss. The 
values for the last three fields for the second crystal are obtained by interpola- 





DANIELS AND RIECKHOFF: SPIN LATTICE RELAXATION 611 


msec 


100 100 





1290 Oersted 


primar lininsei Naess 
| 2 = @&" 








100 100 
10 10 
+ 790 Ocersted +. 2040 Oersted 
b b 
= t 4 1 Base : aot Aiecnpe tices 
' 2 a a I 2 3 4 °K 


Fic. 4. Relaxation time as a function of temperature (log-log plot) for four different fields. 
The solid line drawn through the points has a slope of —3. Open circles, values using pulses 
of microwave power, taken from Table II; half open circles, values using adiabatic demagnetiza- 
tion corresponding to the solid curve in Fig. 3; solid circles, values using adiabatic demag- 
netization on crystal No. 2, obtained by interpolation from Fig. 3. 


tion from Fig. 3, and should not carry much weight; however, they are signi- 
ficantly higher than the values for the first crystal. This is not surprising, 
since there is ample evidence that relaxation times are very sensitive to 
impurities of all kinds. The accuracy of the measurements is too poor to draw 
any definite conclusions about the temperature dependence of the relaxation 
time; the point at 4.22° K, which should carry much weight, is the least 
accurate of all. However, it can be said that the relaxation time varies with 
temperature in a manner somewhere between + « 7~? and +r « 7~*. Our 
results seem to favor an exponent nearer to —3 than to —2, and the lines 
drawn through the points in Fig. 4 have a slope —3. 

Figure 5 shows a curious phenomenon, observed on one crystal which was 
later accidentally destroyed, and which we have not been able to repeat. It is 
seen that, after a long pulse of microwave power, the Faraday rotation in- 
creases to a value greater than that which corresponds to the bath temperature, 
and slowly returns to its equilibrium value. It seems inconceivable that the 
spin system should cool to a temperature below that of the bath and we suggest 
the following as a likely explanation. We assume that the crystal was not a 
single crystal and that in one part the direction of propagation of light, and 
the external magnetic field, were not parallel to the hexagonal axis. We assume 
further that the relaxation time is a function of the angle between the external 
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Fic."5. The ‘‘overshoot phenomenon". Note that the overshoot occurs during the relaxa- 
tion, but not during the saturation period. 


field and the hexagonal axis. When the microwave power is applied, the 
populations of the spin levels are equalized at approximately the same rate 
in both parts of the crystal, since the approach to equilibrium is governed 
principally by the microwave field. When the microwave power is turned off, 
the two parts of the crystal return to equilibrium at different rates. The inten- 
sity of the light passed by the analyzer is a complicated function of time, since 
in at least one part of the crystal the light suffers double refraction as well as 
Faraday rotation. A qualitative analysis of this problem, which it is incon- 
venient to reproduce here, has been given by Rieckhoff (1959) and it is shown 
that such an “overshoot”’ can occur if the relaxation time is shorter in the 
part of the crystal which has the field parallel to the hexagonal axis. 


4. CONCLUSIONS 

We present these results principally as an example of a new technique for 
the measurement of relaxation times. So far, very few measurements have 
been made of relaxation time for one salt over an extensive range of fields and 
temperatures, except for the Dutch work on chromium potassium alum 
(Kramers ef a/. 1950). The methods used to measure relaxation times in the 
past have been (a) that of de Haas and du Pré (1939), in which the real and 
imaginary parts of the a-c. susceptibility are measured at audio frequencies; 
(b) the study of the steady state saturation of paramagnetic resonance 
absorption (Eschenfelder and Weidner 1953); and (c) the use of the strength 
of the paramagnetic resonance signal as an instantaneous spin temperature 
thermometer to follow the return to equilibrium after the spin popu- 
lations have been disturbed by a pulse of microwave power (e.g. Giordmane 
et al. 1958; Davis et al. 1958; Bowers and Mims 1959). The first two 
methods do not observe the relaxation directly, being steady state methods, 
and there are reasons for doubting the results of the method of de Haas and 
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du Pré in some cases. For instance, the locus of the point which represents the 
magnetic susceptibility in the x’ x” plane should be a semicircle, and the 
experimental results of Benzie and Cooke (1950) show curves far from semi- 
circular in some cases. The same criticism might be levied against the method 
of steady state saturation of the paramagnetic resonance signal. Just like the 
method of de Haas and du Pré, its interpretation rests on the assumption of 
the basic relaxation equation (1) and there is in this method no means of 
checking this equation. These results should, therefore, be interpreted with 
caution, as it is now known that the relaxation process is quite complicated. 
What is needed is a means of observing directly and instantaneously the 
populations of the various spin levels. This facility exists in the third method 
and also in our method. Unfortunately, for technical reasons, the use of para- 
magnetic resonance as a spin thermometer limits the measurements to certain 
values of magnetic field for which oscillators are available. The use of the 
Faraday effect is not limited to certain values of the magnetic field, however it 
can measure only an average of the spin populations, whereas paramagnetic 
resonance can measure the population differences between individual levels 
where the resonance absorption. lines are resolved. This is an advantage where 
more than two spin levels are involved, and, for example, where cross pumping 
is to be investigated. The Faraday rotation is, unfortunately, rather small in 
most substances, but there is good reason to believe that with a little attention 
to the instrumentation it should be possible to attain an adequate signal-to- 
noise ratio with the arrangement that we describe. Thus we may claim that 
the use of the Faraday rotation is an interesting new technique for investigating 
paramagnetic relaxation. 

The results are presented merely as measured values; it would be impossible 
to interpret them in terms of a theory of relaxation since there is no adequate 
theory at the present time. The mechanism which is believed to be responsible 
for the relaxation has been described by van Vleck (1940, 1941a, 6). However, 
the agreement between this theory and experiment has always been poor, as 
has also been the agreement between different experiments. For example, 
van Vleck predicts that for chromium potassium alum, the relaxation time 
should vary as 1/7 in the liquid helium range. de Vrijer and Gorter (1952) 
found a dependence as 7~**, while Kramers, Bijl, and Gorter (1950) found a 
dependence as 7~*, and Eschenfelder and Weidner (1953), who found a 
dependence as 1/7, are the only workers whose results agree with van Vleck’s 
theory. So far, no one has succeeded in verifying the predictions of van Vleck’s 
theory with regard to the dependence on field. 

We may review some of the more important assumptions of van Vleck’s 
theory, as these give an indication of the direction in which further efforts 
should be directed. They are: 

(i) That there is thermal equilibrium in the spin system. 

(ii) That there is thermal equilibrium in the lattice. 

(iii) That there is good heat contact between the lattice and the helium bath. 

(iv) That the frequency spectrum of the lattice is a Debye spectrum. 

(v) That spin-spin interactions play no part in the spin lattice relaxation 


process. 
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Items (ii) and (iii) have received attention from many workers, and there 
is still no really convincing experimental evidence on these two points. For 
instance, the conclusions of Eisenstein (1951) and Giordmane et al. (1958) on 
this point are in direct contradiction. Item (iv) should not be very important 
at low temperatures. The most serious assumption is that contained in (i) and 
(v). We can say generally that in concentrated paramagnetic salts the inter- 
actions are strong enough to establish thermal equilibrium in a time of the 
order of 10~* seconds, but that this is by no means true in dilute salts, so that a 
distinction must be drawn between these two cases. Where there is a strong 
interaction between the spins, it is unreasonable to expect no influence on the 
relaxation time, and it is now well known that if the separation between two 
levels is the same as that between another two levels, there is a ready transfer 
of energy between systems whose state is one of the first pair, and systems whose 
state is one of the second pair. This produces a shortening of the relaxation 
time, and is the basis for the control of relaxation times in masers by ‘“‘doping’’. 
In our experiments, we must remember that neodymium has two isotopes with 
nuclear spin different from zero: Nd1438, 12.2% abundant, and Nd145, 8.3% 
abundant, both with spin 7/2. The ions containing these isotopes have a 
hyperfine structure, and the condition for a three-spin process occurs at about 
800 oersted. This may explain the initial dip in the relaxation time as a function 
of magnetic field (Fig. 3). 

We feel that the theory of relaxation is in too poor a state, and that our 
measurements are not sufficiently accurate to attempt any further explanation 
of our results. 
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TRANSISTORIZED POWER SUPPLIES FOR A MASS 
SPECTROMETER! 


R. D. RUSSELL AND FRANK KOLLAR 


ABSTRACT 


While a mass spectrometer itself may require a few hundred watts of regulated 
power for the magnet and only a few watts of regulated power for all other 
purposes, the regulating supplies often dissipate kilowatts. The construction 
of more efficient supplies using transistors and magnetic amplifiers results in 
substantial savings in the physical size of the supplies, the cost of their com- 
ponent parts, and the effort required for maintenance. 

Circuits and descriptions are presented for a magnet current power supply 
and a filament emission control that have been used in the mass spectrometer 
built recently at this laboratory. These illustrate the simplicity possible with 
transistor-regulated supplies, as well as the savings mentioned above. The mass 
spectrometer is a 12-in. radius 90° instrument used at mass 250 with a 5000-volt 
accelerating voltage. The magnet, which supplies a field of more than 5000 
gauss over an area of 800 cm? and across a gap of 1.9 cm, and which weighs three 
quarters of a ton, is supplied from a regulated supply dissipating about 50 watts. 
The filament emission control takes about 30 watts from the mains. 


INTRODUCTION 

The electronic supplies associated with mass spectrometers have shown 
remarkably few fundamental changes in the past fifteen years, despite the 
very great advances that have been made during that time in the field of 
applied electronics. After working for a number of years with conventional 
mass spectrometer electronic circuits and then for the last year with an instru- 
ment using transistors and magnetic amplifiers, the writers are completely 
convinced of the superiority of these more modern components in almost 
every mass spectrometer circuit, and of the fact that they are not more widely 
used now simply because most mass spectrometrists are unfamiliar with 
them. The object in our making these circuits available is primarily to demon- 
strate some of the advantages of transistor-operated power supplies in mass 
spectrometers. 

In the construction of a mass spectrometer at this university, it was decided 
to use transistor regulators for both the magnet current power supply and the 
. filament emission control. The magnet current power supply was constructed 
from a combination of transistor, vacuum tube, and magnetic amplifiers; the 
filament emission control uses only transistors. Both supplies have proved 
efficient, rugged, and inexpensive, and both have required negligible main- 
tenance. 

The mass spectrometer itself is a general purpose gas-source mass spectro- 
meter, intended initially for the analysis of lead tetramethyl. It uses a copper 
tube having a 12-in. radius of curvature, and an accelerating voltage of 
5000 v. 

1Manuscript received January 4, 1960. 

Contribution from the Geophysics Laboratory, Department of Physics, University of 
British Columbia, Vancouver 8, B.C. 
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A modified grid-ratio measuring system is used to measure the ion current 
intensities. This incorporates a servo voltmeter of 0.025% linearity and 
reproducibility. The measuring system will be described in a future publication. 


THE MAGNET CURRENT SUPPLY 

The magnet is constructed from a mild steel (‘‘boiler plate’) yoke with 
armature and pole pieces from Armco iron. The area of each pole piece is 
800 cm? and the distance between them is 0.7500+.0005 inch. The reluctance 
of the magnetic circuit is substantially that of the air gap and therefore the 
magnetizing curve is linear except at very high fields. The magnet coils con- 
tain a total of 21,000 turns of No. 18} A.W.G. copper wire having a total 
resistance, when connected in series, of 720 ohms. The inductance of the 
magnet calculated from the above figures is 2100 henrys. However, the effect 
of eddy currents becomes significant at frequencies higher than about 
0.2 cycle per second. The impedance of the magnet is shown as a function of 
frequency in Fig. 1. The data for this figure was obtained at higher frequencies 
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Fic. 1. Approximate impedance characteristics of the magnet measured at its terminals. 


by determining the amplitude and phase of the current produced by a small 
sinusoidal driving voltage, and at lower frequencies by analyzing the response 
due to a voltage step. The characteristic applies at zero d-c. current and can 
be expected to be different in detail at high fields. The magnet characteristics 
are a dominant term in the over-all transfer function of the magnet current 
power supply and the knowledge of its form is necessary to obtain optimum 
stability and speed of response of the over-all circuit (see Savant (1958) for a 
general discussion of this subject). 
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The maximum magnet current used is 500 ma, at which current the magnet 
voltage is 360 v and the power dissipated is 180 w. The measured field is 
5500 gauss at this current, sufficient to analyze masses to nearly 300 with the 
maximum accelerating voltage. 

The circuit for the regulator is shown in Fig. 2. It incorporates two regu- 
lating devices, the saturable reactor which is a show device, and the transistor 
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Fic. 2. Magnet current supply. Condenser values are given in microfarads, resistors in 
ohms. Zener diodes are identified by their breakdown voltages. Transformer and reactor 
numbers in parentheses are Hammond type numbers. The condenser C is chosen for minimum 


ripple (0.2 uf). 


circuit, which has a fast response. One important function of the saturable 
reactor circuit is to keep the voltage across the control transistor at a suitable 
level. However, it also contributes significantly to the regulating action. Thus 
the regulation is seriously impaired if the saturable reactor is removed and the 
power transformer is supplied a suitable fixed input voltage by means of a 
variable voltage transformer (e.g. a Variac). 

The supply operates as follows. The voltage of the mains is fed through an 
auto transformer to the primary of the power transformer connected in series 
with the saturable reactor. The auto transformer adjusts the supply voltage 
close to the optimum value for the different current ranges employed. Since 
the saturable reactor can itself produce changes in the transformer voltage 
by approximately a factor of 10, the auto transformer could probably be 
omitted, but the regulation would be somewhat impaired, particularly at 
very high and very low currents. A full wave rectifier is employed with a 
condenser input filter. The filtering circuit used is in the same servo loop as 
the saturable reactor, which has a long time constant resulting from the large 
inductance of the d-c. winding. A conventional LC filter introduces additional 
phase lags and causes a deterioration in the transient characteristics of the 
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over-all system. Thus the tuned 120-cycle rejection circuit shown is used to 
reduce ripple. The remaining ripple, the order of 1 volt, is easily handled by 
the transistor regulator circuit. The output from the filtered power supply is 
connected to the magnet windings in series with the 2N457 control transistor 
and a series resistor, the voltage drop across which resistor serves as a 
measure of the magnet current. The sum of the voltages across the control 
resistor and transistor regulates the bias of the triode operated 6AQ5 vacuum 
tubes controlling the saturable reactor current. This serves to hold the voltage 
drop across the control transistor to a reasonable value (about 30 volts) at 
all magnet current ranges, and well below the 60-volt maximum rating of the 
transistor. The control transistor is rated at a current of 5 amp maximum and 
a power dissipation of the order of 50 watts. Therefore the maximum magnet 
current that can be delivered by the supply is here limited to 500 ma only by 
the transformer, choke, and saturable reactor ratings. 

Two 2N369 transistors are used as a difference amplifier, comparing a 
fraction of the voltage developed across the series control resistor with the 
5.36 v produced by four mercury cells which serve as the basic reference. 
This reference voltage is high-enough so that pickup and drift in the difference 
amplifier inputs can easily be made negligible. We feel that the use in this 
circuit of a very low reference voltage (50 mv) as suggested by Garwin (1958) 
and Garwin et a/. (1959) would not only exaggerate problems of noise but 
sacrifice loop gain needlessly and therefore require additional gain in the 
electronic circuits increasing the problems of stability. In this circuit the 
mercury cells are quite stable and have a long life, and therefore seem to be a 
convenient reference; if the use of batteries were considered undesirable, a 
Zener diode would provide an excellent reference. A third 2N369 isolates the 
power transistor from the difference amplifier. This transistor is essentially 
an emitter follower, matching the a-c. and d-c. impedances of the difference 
amplifier to the power transistor. The 10-v Zener diode stabilizes the collector 
supply voltage for one transistor in the difference amplifier otherwise the 
gain would be greatly reduced by negative feedback. The difference amplifier 
is designed so that the collector voltages are almost identical at all operating 
currents. The 20-v Zener diode adjusts the operating point of the 6AQ5, and 
thereby adjusts the range of voltages across the control transistor. The 22k 
resistor provides some positive feedback, thereby increasing the loop gain. 

The 54-v Zener diode protects the circuit against switching or other transients 
that might result in exceeding the voltage rating of the 2N457. The M-500 
Sarkes Tarzian diodes protect the magnet windings from a voltage transient 
when the supply is shut off. Many diodes would be equally suitable, but in 
any case best results can be obtained by measuring the reverse characteristics 
and selecting diodes of low reverse current. The characteristics of these silicon 


diodes were quite variable. 

In the construction of the unit, care must be taken to eliminate a-c. pickup 
in the inputs to the difference amplifier. In particular it was found worth 
while to use bifilar-wound control resistors. The transistor amplifier was 
built on a heat sink made from 3 in. X 2 in. X 1/8 in. aluminum channel, 
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4 inches long, which plugs into an octal socket. The performance of the remain- 
der of the circuit can be easily checked by unplugging this unit and replacing it 
by a variable power resistor of 50-300 ohms resistance. Thus the saturable 
reactor circuit can be checked independently of the transistor circuit. Similarly 
replacing the 6AQ5 tubes by a suitable resistor permits the transistor unit to 
be tested independently of the saturable reactor control loop. The a-c. power 
for the supply is obtained from a Sola regulating transformer and the d-c. from 
a 250-v Sola d-c. supply. These supplies were available for other parts of the 
mass spectrometer; otherwise the d-c. could be obtained from a simple circuit 
with silicon rectifiers. The performance of the supply does not deteriorate very 
much if the Sola regulating transformer is not used. 

A 10-v step in the input line voltage (without the a-c. Sola regulator) 
results in a 0.012% change in magnet current. The stability of the current is 
in the order of | part in 25,000 for variations of periods in the order of seconds 
or minutes, and the ripple current is less than this if the filter is reasonably 
well tuned so that the transistor amplifier is not overdriven. This is not only 
because the ripple voltage is greatly reduced by the a-c. degenerative feedback, 
but also because of the high coil impedance at 120 c.p.s. (see Fig. 1). 

The difference amplifier transistors are individually quite temperature sensi- 
tive. The slow drift resulting from this can be minimized by arranging that 
the temperatures of these transistors are always equal. A further improvement 
in long term stability might be achieved with silicon transistors in this part 
of the circuit. 

Power losses are quite small. At maximum current there are about 15 watts 
dissipated on the power transistor, 5 watts on the control resistor, 5 watts 
on the 6AQ5 filaments, and about 25 watts in the 6AQ5 plate circuit. The 
power delivered to the magnet is about 180 watts, giving an over-all efficiency 
of more than 75%. The cost of the electronic components was in the order 
of $150 not including the meter, which is a mirror scale instrument of 1% 
accuracy imported from Denmark at a cost of $40 and the magnet coils 
which were manufactured commercially at a cost of $650. The meter has, in 
addition to the current scale, a scale proportional to the square of the current, 
_ which is useful as a mass indicator. 

A number of modifications could be considered in the design. The use of 
400 cycles or higher frequency power is possible and at present we feel this 
would be worth while. The substitution of transistors for the 6AQ5 is an 
obvious change that would require only a change of resistance of the d-c. 
winding of the saturable reactor. Thyratron regulators could substitute for 
the saturable reactor, as grid-controlled rectifiers. Since the supply was built, 
silicon-controlled rectifiers have appeared on the market and these also seem 
to offer an attractive substitute for the saturable reactor. In either case the 
speed of response and gain of the ‘‘slow’’ loop could be greatly increased, and 
this in turn would make the filter design less critical. 


THE FILAMENT EMISSION CONTROL 


The circuit for this unit is shown in Fig. 3. It is a straightforward regulated 
power supply using a series transistor as the regulating element. This trans- 
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Fic. 3. Filament emission control. Component identification is the same as in Fig. 2. 


Three-phase transformer secondaries are 9-v center-tapped windings insulated to 5 kv. The 
resistance of the filament and 5A meter, including interconnecting cables, must not exceed 
1 ohm, unless the secondary voltages of the three-phase transformer are increased accordingly. 


istor can easily regulate the 5-amp current used at a voltage drop of less than 
2.5 volts. Obviously no vacuum tube will operate at these currents and voltages 
and therefore this type of circuit is only possible with transistors. 

A three-phase 120/208-v service was available at this laboratory, and 
therefore a three-phase input — six-phase output transformer was used for 
the filament supply to give a low ripple without additional filtering (approxi- 
mately 4.5% r.m.s. in the d-c. voltage). A similar result could be obtained by 
using a full wave single phase rectifier with a large condenser (several tens 
of thousands of microfarads) or an LC filter. In some cases the large ripple of 
an unfiltered supply might not be considered objectionable. 

The 2N456 control transistor is connected with a second similar transistor 
in the Darlington compound arrangement (Shea 1957) to provide a large value 
of current gain and therefore a sufficiently low base current (typically 0.7 ma). 
This pair is driven by a 2N369 voltage amplifier, which in turn is driven by a 
2N366 emitter follower. The forward voltage drops across the emitter-base 
diodes of the latter two transistors approximately cancel and therefore at 
balance the error signal, which is applied between the base of the 2N366 and 
the emitter of the 2N369, is very nearly zero. The error voltage is the differ- 
ence between the voltage developed across the 100k helipot and its series 
resistor, and the reference provided by the 27-v Zener diode. 

In addition there are two supplies of very low power. One is a glow-tube- 
stabilized supply to provide the electron accelerating voltage and the trap 
bias in the mass spectrometer source, as well as providing current for the 27-v 
Zener. This supply uses the smallest available filament transformer connected 
backwards so that 7.2 v r.m.s. from the three-phase transformer is applied 
across its 6.3-v winding, the voltage developed in the primary being rectified 
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with silicon diodes in a simple doubler circuit. This arrangement avoids the 
necessity of having a second transformer insulated for the 5000-v accelerating 
voltage. The second supply is a low voltage three-phase silicon rectifier circuit 
which supplies the emitter follower. The three-phase power was available 
from the main transformer secondary, and using it simplified the filtering. 

This supply is physically the smallest used in our mass spectrometer. The 
electronic units, including the two subsidiary power supplies and the required 
heat sinks for the power transistors and silicon rectifiers fit entirely into a 
chassis 2 in. X 8 in. X 10 in. This chassis floats to 5000 v when the mass 
spectrometer is operating and is suspended inside a somewhat larger grounded 
chassis, thus eliminating the possibility of corona discharges around small 
components and providing increased safety. The transformer is mounted on 
the subsidiary chassis and the meters are mounted separately on the front 
panel with the various controls. The supply input power is less than 30 w. 
The supply cost approximately $100 plus the cost of the meters. 

The regulation of the supply is quite adequate, the measured r.m.s. variation 
in electron current being in the order of one or two parts in 10,000. The supply 
requires no warm-up time and exhibits remarkably little drift in view of the 
fact that the base current of the 2N366 shuts the 100k control helipot. The 
electrolytic condenser in parallel with the error signal is necessary for a 
well-damped system and has not caused any of the difficulties we had antici- 
pated, perhaps because the maximum voltage across it is tens of millivolts.* 


CONCLUSIONS 

Certain of the regulated power supplies normally used with a mass spectro- 
meter lend themselves well to conversion to transistor-operated circuits. Nor- 
mally the actual power used by the mass spectrometer itself is very small, 
while the regulating supplies dissipate several kilowatts. It is to be expected 
that an increased efficiency would not only simplify problems of construction, 
but require fewer expensive components and less maintenance. Therefore 
it is well worth the effort to modernize the usual mass spectrometer power 
supplies. 

Two supplies have been described in this paper, a magnet current power 
‘supply and a filament emission control. These supplies illustrate well the 
relative simplicity of transistor-operated mass spectrometer circuits and have 
now been used long enough in our mass spectrometer so that we can be con- 
fident that they will continue to work well in their present form. Although 
both supplies are quite adequate for our needs, they would probably have to 
be modified somewhat for different mass spectrometers. Certainly the regula- 
tion and stability of either could be further improved if this were necessary. 
Circuits we are presently designing for a second mass spectrometer differ 
from these primarily in using transistorized d-c.-d-c. converters at higher 
frequencies (hundreds of cycles per second) to utilize still smaller and less 
expensive components, and obtain better servo characteristics in some cases. 


*We later replaced this condenser with a 4-uf condenser between base of the 2N366 and the 
negative side of the filament. This further improved stability and reduced ripple. 
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CORRELATION STUDIES OF RADIO-AURORA, 
MAGNETIC, AND EARTH-CURRENT DISTURBANCES! 


B. K. BHATTACHARYYA? 


ABSTRACT 


Correlation studies of the radar echo occurrence rate from aurora in half- 
hourly intervals at Ottawa, S and Sy components of the horizontal magnetic 
field H at Agincourt, and the disturbance diurnal variation of earth current 
at Crow River have been carried out. Short-time variations in auroral echo 
strength and moderate perturbations in H have also been correlated. The auroral 
echo occurrence rate seems to have a diurnal variation characteristic similar to 
that of H. It is found that auroral activity always precedes magnetic activity. 
The variation of the delay time between the two phenomena shows a local time 
dependence, being practically constant and quite small (0-15 minutes) before 
local midnight and increasing afterwards. This variation of the delay time 
appears to have a connection with reports of others regarding reversal of the 
direction of auroral ionization drift from west to east somewhere around mid- 
night with subsequent magnetic perturbations which change from positive to 
negative. 

No definite conclusion could be reached regarding the relationship of earth 
current to other factors because of a practically random variation of cross- 
correlation coefficients from month to month. 


INTRODUCTION 

The interrelationship between magnetic disturbance and auroral activity 
has been known for a long time. More recently, attempts have been made to 
establish a systematic relationship between the two phenomena with the help 
of visual, photographic, and radio data of aurora on the one hand, and mag- 
netograms recording the three components of the magnetic field on the other. 

Firstly, with the help of visual aurora data, striking relationships between 
the two phenomena have been obtained. The control of the degree of magnetic 
disturbances over the maximum elevation of aurora above the northern horizon 
at Saskatoon has been shown (Meek 1953, 1954). A relationship between the 
3-hour range K index and the auroral position in the sky at Ithaca has been 
empirically derived (Gartlein 1944). Time sequences and spatial relations in 
auroral activity during magnetic bays at College Alaska have been studied in 
great detail (Heppner 1954). Atmospheric electric currents have been shown 
to be associated spatially with aurora (Stagg and Paton 1939; Heppner 1954). 
Very recently, evidences (Bless et al. 1959) have been produced to show that 
magnetic bays and auroras are manifestations of the same phenomena. 

Secondly, photographic auroral data have been utilized to study mainly the 
motions of visible auroral forms. Attempts have been made to find the local 
time dependence of the change in the direction of motion (Meinel and Schulte 
1953; Meek 1954; and Bless et al. 1955). Most of the studies indicate a reversal 


the direction of the drift velocity from west to east near local midnight. 
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However, this reversal of the direction of drift appears to be related more to 
the character of aurora and of the magnetic disturbance than to the local 
time (Nichols 1959). 

Thirdly, from the study of radio auroral data it has been suggested that the 
drift motion of auroral ionization as observed by radar is probably the motion 
of electrons associated with the magnetic disturbance current system (Nichols 
1957; Unwin 1959). Bowles (1954) has shown that these drift motions occur 
along the magnetic lines of force while Bullough and Kaiser (1955) have 
concluded that they are either eastward or westward along parallels of geo- 
magnetic latitude. The latter workers have also found a reversal of the direction 
of motion from west to east between 2100 and 2200 U.T. Recent measurements 
by Nichols (1957), however, have disproved the observation by Bowles and 
convincingly indicated that the motions are generally horizontal and in the 
geomagnetic east-west plane. The speeds of the motions as determined by 
Nichols vary from 350 meters per second to about 3000 meters per second. 

Since the earth currents are supposed to be normally the effect of the over- 
head current systems responsible for magnetic field perturbations, correlation 
studies of aurora, magnetic, and earth-current activities have also been 
attempted (Currie and Edwards 1936; Rooney 1934). 

This paper is an outcome of studies made to probe the gross and detailed 
relationships between certain characteristics of radio aurora, variations of the 
horizontal magnetic field component HH, and earth current. 

As regards the gross effects, correlation studies are made between each pair 
of the following: 

(a) the duration of radio auroral echo in half-hour periods; 

(b) the half-hourly mean daily (S) and disturbance diurnal (Sy) variations 

of 77; and 

(c) the half-hourly mean disturbance diurnal variations of the earth current. 
A statistical study of the occurrence rate of radio aurora and magnetic dis- 
turbance indices, the 3-hour range A index, and SA index for the whole day 
is also made. 

lor the study of detailed relationship, separate individual moderate magnetic 
disturbances are chosen and the envelopes of the individual disturbance and 
corresponding auroral echo, if anv, have been correlated. 

The radio auroral data have been obtained with the help of a 50 Mes 
radar designed for back-scatter detection of auroral echoes at Ottawa (45.4° N 
lat. and 75.9° W. long.) (McNamara 1958). The time duration of radar echo 


from aurora in half-hourly intervals throughout all days in L958, hereinatter 
referred to as ‘echo duration’, has been determined from the photographic 
records taken in the form of a continuous film strip displaying range and time 


The instantaneous radar echo strength has been measured trom the Esterline 
Angus chart records displaving the integrated echo intensity with time 
Phe magnetic data have been collected from the magnetograms otf 
Dominion Observatory of Canada taken at Agimcourt (43.8° N. lat. and 
79.3° W. long.) with the help of a standard photographic variometer. Uhe halt 


hourly average values of the horizontal component #/ of the magnetic treld 


have been scaled tor all days in LOSS. 
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The earth-current records are taken at Crow River (51.5° N. lat. and 
90.1° W. long.) by the Hydro-Electric Power Commission of Ontario with the 
help of a d-c. potentiometer connected to measure the d-c. voltage across a 
portion of the power transformer neutral connection to ground. The other end 
of the transmission line, which is 161.7 miles long, is at Ear Falls (50.6° N. lat. 
and 93.2° W. long.), approximately 22.5° south of west from Crow River. 
The half-hourly range of deviation of the earth current from the mean has 
been read from the chart records for all days from July to December, 1958. 


DIURNAL VARIATION CURVES OF ECHO DURATION AND H 


With the help of the data available the monthly half-hourly mean values are 
calculated and then plotted in Fig. 1. The following results were found from 
the figures. 

(a) The diurnal minima in H and echo duration are conspicuous and occur 
around 1600 hours U.T., which is very close to the geomagnetic noon. In 
winter months the minima tend to occur at times later than those in summer 
months and in equinoxial months at an intermediate time. This is true in the 
case of both sets of data. In most of the months an auroral echo does not appear 
at all around local noon, which is close to geomagnetic noon. 

(6) There are two maxima in both sets of curves, one is very prominent and 
the other quite minor. The prominent maxima occur mostly between 2200 and 
2400 in both cases. This time seems always to be ahead of the geomagnetic 
midnight. The other maxima in case of H occur mostly between 0900 and 1100 
and in the case of echo duration between 0500-0800 (early morning). 
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Fic. 2. Linear correlation coefficient between echo duration and S component of H for 


every month in 1958. 


Krom Fig. 1 it is found that both // and echo duration have some strikingly 
similar features in their diurnal variation components. To determine an index 
of the dependence of one on the other, linear correlation coefficients have been 
calculated and plotted in Fig. 2. Except for some summer months, in particular 
June, July, and August, the correlation coefficient in quite high, above 60%. 
In the equinoxial months the correlation is highest, about 75% on the average. 
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DISTURBANCE DIURNAL VARIATION OF H AND ECHO DURATION 


In determining the disturbance diurnal variation, magnetically quiet days 
are selected with the help of the 3-hour range K indices and also the magnetic 
data available. In the selection of quiet days, attention is paid not only to the 
YK value but also to the absence of any moderate disturbances. The number 
of quiet days varies from month to month, 9 in November, 6 in April and 
June, and 5 in the other months of the year. Then the Sy component of H is 
calculated for every month of the year. These have been plotted in Fig. 1. 

The linear correlation coefficients between the Sg component of H and the 
auroral echo duration for every month of the year is then determined. The 
value of the correlation coefficient on the average is below 30% except for the 
month of August, where it is 66%. This leads to a very surprising result— 
that there is practically no correlation between the auroral echo duration and 
Sq field—a result contrary to expectation. 

The possibility of some delay time between the two phenomena, which may 
be responsible for the above result, is then explored by calculating the maximum 
cross-correlation coefficients after the introduction of optimum delays in the 
two sets of data for every month. The coefficients and the delays thus obtained 
are plotted in Fig. 3 for every month in 1958. The magnitude of the cross- 


f 





5 « S > 2 


Fic. 3. Linear cross-correlation coefficients and delay times between echo duration and 
Sa component of H for every month in 1958. 


correlation coefficient is quite high throughout most of the year, except in 
January when it is —38%. It is found that radar auroral activity at Ottawa 
always precedes the magnetic activity at Agincourt. The delay time between 
the two activities undergoes a seasonal change. In summer months this 
fluctuates very rapidly but near the equinoxes, particularly during March and 
April, it is more or less steady. The latter is also evident in the seasonal 
characteristic of the cross-correlation coefficient which is negative in mid- 
winter, during December and January, and decreases in magnitude in mid- 
summer during June and July as compared to the equinoxial months. 
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SHORT-TIME CORRELATION OF INDIVIDUAL MODERATE 
PERTURBATIONS IN H AND CORRESPONDING ECHOES 

As already stated, the auroral activity always occurs earlier than the 
magnetic activity. This suggests that the disturbances in H are mostly the 
result of the intense current flowing in the auroral zone which is proportional 
to the number of charged particles present therein. This number again deter- 
mines the radar echo strength (Moore 1952; Booker 1956). Hence a correlation 
between the instantaneous variations of echo strength and // may be expected. 

In calculating this correlation we are faced with one difficulty. Whereas the 
radar echo strength is always positive in sign, the sign of the magnetic field 
perturbations in relation to the average is sometimes positive and sometimes 
negative. This appears to be due to the change of the direction of the drift of 
auroral ionization produced by incoming solar protons (Heppner 1954). 

Whatever may be the cause, if magnetic field perturbations are caused by 
currents in the auroral zone, the instantaneous magnitudes of the negative 
perturbations will correlate well with the corresponding echo strength and the 
cross-correlation coefficient will be negative. In the case of positive perturba- 
tions the coefficients will always be positive. Hence the above coefficients 
between echo strength and magnitude of the magnetic perturbations have 
been determined. This should give information about the relation between the 
instantaneous intensities of individual perturbations in /7 and those of the 
corresponding echoes, more particularly an answer to the question regarding 
the cause-and-effect relationship between the two phenomena. 

At Ottawa a continuous record of the integrated echo strength (both with 
respect to range and time) is recorded on Esterline~Angus charts (McNamara 
1958). From these charts days are selected which have a number of separated 
echoes of moderate peak intensity and of duration about 2 hours, distributed 
practically all throughout the day. In spite of this selection, there is a period 
of about 8 hours near midday when echoes are absent because of a predominant 
tendency of aurora to occur at night. Each of these echoes is then scaled at 
intervals of 5 minutes. From the magnetograms the //-component at and 
around the time of occurrence of each radar echo is also scaled at similar 
intervals. Since for a time interval of 2 hours or so, the diurnal variation of the 
earth’s magnetic field will not have any appreciable effect, no correction is 
made to the scaled values of /7. 

In all cases it is noted that around the times of appearance of an auroral echo 
there are perturbations in the /7-component. The peak intensities of moderate 
perturbations do not seem to be in the same ratio as that of the corresponding 
peak echo intensities in most cases studied. This is not unexpected because 
factors determining the echo strength, such as echo range variation and the 
effect of the polar pattern of the radar antenna, and factors responsible for 
magnetic field perturbations, such as the direction of current flow in the 
auroral zone, are not all the same and have not identical effects on both the 
phenomena. Nevertheless, very high amplitude echoes are always found to 
accompany severe magnetic storms. 
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The cross-correlation coefficients are then determined with the help of two 
sets of scaled values for a number of days. They are practically always found 
to be quite high, above 65% most of the time and some above 80%. Only in a 
few cases are they about 50%. 

An auroral echo is always found to precede the associated magnetic pertur- 
bation. The delay times obtained have been plotted against universal time for 
a few days in Fig. 4. In this figure the delay time between the two events has 
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Fic. 4. Delay time between individual auroral echoes and magnetic perturbations against 
universal time. The thin vertical lines denote that the delay times are between 0 and 5 minutes. 


been plotted against the universal time midway between the duration of the 
radar echo. In Fig. 4 we find that starting from the afternoon hours the delay 
time remains practically constant between 0 and 10 minutes until about local 
midnight (0500 U.T.). Thereafter the delay time tends to increase. 

This variation of delay time at about the local midnight seems to be related 
to the change in the direction of auroral ionization drift. From both photo- 
graphic and radio studies of the drifts of auroral forms and ionization referred 
to in the introduction, it has been found that the direction of drift before 
midnight is usually westward causing positive bay type perturbations in the 
magnetic field, whereas after midnight the direction is reversed, causing 
negative bay type perturbations. This may be due to a change in the polarities 
of the electromotive force driving the current in the auroral zone. The delay 
time between the auroral echo occurrence and associated magnetic perturba- 
tions may be the delay between the setting-up of the electromotive force and 
the full growth of the current. 


AURORA AND BAYS 
Statistics of auroral echo occurrence with geomagnetic bays may be sum- 
marized in the following way for the year 1958. 
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Magnetic elements H D Z 
Sign of the disturbance a —- + al ee rote 

No. of days of disturbance 11 35 46 43 59 7 

No. of days aurora recorded in radar film 9 29 35 26 50 «6 


% of disturbed days on which aurora occurred 82 83 76 60 85 86 





The cross-correlation coefficients between two negative bays in H of duration 
close to 2 hours and auroral echo strengths at about the corresponding times 
are found to be above 80%, and the delay times between the two events are 
less than 5 minutes. 


AURORAL ECHO DURATION AND THE 3-HOUR RANGE K INDEX 
Correlation studies of auroral echo duration with magnetic disturbances as 


deduced from the chart of 3-hour range K indices in 1958 may be summarized 
as follows: 





Percentage of total 
disturbed periods 





Range of variations K auroral echo is 
in H, gamma indices received 
(i) >330 8,9 79 
(ii) 200-330 7 74 
(iii) 120-200 6 65 
(iv) 70-120 5 53 
(v) 40-70 4 37 
4 6.5 


(vi) 0-40 < 


The percentage of time of the total disturbed periods is shown as a histo- 
gram in Fig. 5 against the different range of variations of H. The relation 
between the two appears to be exponential. 


AURORA 


4 1 1 i 4 
5 10 20 40 70 120 200 330 500 
RANGE OF VARIATIONS IN H————® 


Fic. 5. Percentage of total disturbed periods on which auroral echo was received as a 
function of the range of variations of H. 
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If the days are divided into different categories according to the magnitude 
of 2K for the whole day, the statistics of observations may be put in the 


following way: 








Percentage of days 
auroral echo is 





=K No. of days received 

0-10 84 1.5 
10-20 149 8.36 
20-30 88 21.3 


>30 44 42.8 
The nature of variation of the right-hand side with the 2A values is also 
practically the same as shown in Fig. 5. 

From the plot of percentage of time of auroral echo duration and A,, for 
every_month in 1958 (Fig. 6), it is observed that there is practically a one-to- 
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Fic. 6 


one correspondence between the two sets of curves in winter, a moderate 
amount of correspondence in the equinoxial months, and a poor correspondence 
in the summer months. 


ECHO DURATION, S AND Sa VARIATION IN H, AND DISTURBANCE 
DIURNAL VARIATION IN EARTH CURRENT 

The earth-current records obtained from Crow River — Ear Falls stations 
are quiescent on magnetically quiet days and so the half-hourly mean ranges 
of maximum deviation from the mean for a whole month will give the dis- 
turbance diurnal variation of the earth current. This variation for the months 
of July to December in 1958 is plotted in Fig. 1. 

The atmospheric current system responsible for magnetic perturbations 
gives rise to earth currents by an induction process and so the earth current 
will appear to be related to the time derivative of the magnetic field. Since 
the data on the rate of change of the field at Agincourt is not available, and 
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it is extremely time-consuming to compute the above from the magnetograms, 
direct correlation studies between the disturbance diurnal variations of earth 
current and the magnetic field have been attempted, with expectation of a 
phase shift and consequent delay time between the two phenomena. The 
result of computation is given below: 











Delay time, Cross-correlation 
Month hours coefficient 
July —3 56.37 
August —-10 —55.95 
September +4 — 61.64 
October 0 53.88 
November 2 — 63.06 


December 0 53.37 


In the above table the negative delay time indicates a precedence of earth- 
current activity over the magnetic activity and the positive one the reverse. 
The physical significance of these results seems to be perplexing. In particular, 
in the summer months of July and August, the magnetic activity seems to lag 
behind the earth-current activity. This may be due to some sort of complex 
phase shifts introduced by the intervening medium between the atmospheric 
and earth-current systems. Moreover, the two measurements, Crow River — 
Ear Falls for the earth current and Agincourt for the magnetic field, are quite 
far apart and so a number of hypothetical reasons may be put forward. Never- 
theless, it has been found that earth-current activity is always associated with 
moderate and severe magnetic disturbances though the amplitudes of earth 
current and horizontal magnetic field are not very well correlated. In earth- 
current records micropulsations of small period, normally less than 5 minutes, 
are always present. These pulsations cannot, in most of the cases, be accounted 
for by the variation of the horizontal magnetic field component at correspond- 
ing times. 

It is always found that the auroral activity antedates the earth-current 
activity. The results of the correlation study between earth-current activity 

_and auroral echo occurrence activity are as given below: 


Delay time, Cross-correlation 


Month hours coefficient 
July 12 67.23 
August ‘< 59.59 
September 43 — 48.07 
October 5} 49.46 
November 64 84.58 
December 134 72.49 





It seems from both tables that in November, the magnetically quietest 
month in 1958, the cross-correlation coefficients are highest. This may be due 
to the least effect of extraneous disturbances on earth current in that month. 
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SUMMARY AND DISCUSSION 
Studies in this paper have shown that auroral echo activity always antedates 
the magnetic activity. This appears to be the reason why the mean disturbance 
diurnal variation of H always lags in time behind the mean auroral echo 
occurrence rate. The delay time has a diurnal variation and its monthly mean 
value betrays a marked seasonal change. 
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Fic. 7. An example showing development of a negative bay and a corresponding radar 
echo on June 2, 1958, at 06:38 U. T. 
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It seems, therefore, that the auroral phenomenon is the primary effect of 
the intrusion of solar charged particles into the upper atmosphere and this is 
then the origin of most of the magnetic disturbances. The latter is caused by 
the current flow in the auroral zone due to the drift of protons in the ionized 
gas with the surrounding wind, as postulated by Bless et al. (1959) in their 
study of relationship between bays and aurora. They have shown that only 
protons will move with the wind and they are capable of giving rise to magnetic 
bays. But it is not clear from their study how electrons, which have been 
found to move with a velocity as high as 500 m/s (Bullough and Kaiser 1955; 
Nichols 1957; Unwin 1959), are not effective in producing the bays. Their 
idealized analysis only shows that ba.» would have been much larger in 
magnitude than observed if the drift 0: .ectrons were responsible. 

The close association of bays and aurora has been lately reported by many 
workers. At Saskatoon (52.1° N. lat. and 106.6° W. long.) in Canada Meek 
and McNamara (1954) have shown that 56 and 106.5 Mc/s radar echoes from 
aurora are associated very closely in time with positive bays whereas no such 
association is detected with the negative bays. Zaborshchikov and Fediakina 
(1957) have also made similar observations at Tiksi Bay Arctic Station (71.6° N. 
lat. and 128.9° E. long.) in the U.S.S.R. This has been suggested as due to the 
occurrence of aurora at higher angles of elevation during negative bays, and 
the fact that radar echoes are not normally received from aurora at high angles. 
The above suggestion has been supported by the report of Meek (1954) that 
auroral light associated with positive bays occurs at a higher geomagnetic 
latitude than that associated with negative bays. 

In this study disappearance of auroral echoes with the onset of a negative 
bay has not been noticed. On the contrary, the auroral echo strength has been 
found to increase up to its maximum value and then decrease in a gradual and 
systematic fashion with the formation of a negative bay (Fig. 7). The high 
negative correlation coefficients between the envelopes of bays and corre- 
sponding auroral echo strengths justify our observation. The cause of this may 
be that at Ottawa, which is more southerly than the other two stations in 
geographic latitude, aurora associated with negative bays is not displaced to 
such high angles of elevation and the reflection of radio signals is still possible. 
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HIGH-ALTITUDE COSMIC RAY MEASUREMENTS AT 
FORT CHURCHILL! 


I. B. McD1armMip AND D. C. ROSE 


ABSTRACT 


Measurements with rocket-borne Geiger counters have been carried out at 
altitudes up to 250 km at Fort Churchill, Manitoba. The total primary cosmic 
ray intensity at a time near a solar maximum has been determined and com- 
pared with other measurements taken at times of high solar activity and also with 
other Geiger counter measurements obtained near a solar minimum. A low-energy 
radiation was observed whose intensity increased with altitude up to about 
25% of the primary intensity at 250 km. 


INTRODUCTION 


Cosmic ray detectors were included in the instrument loads carried in two 
Defence Research Board Aerobee-Hi rockets launched in September 1959 at 
Fort Churchill, Manitoba (geomagnetic latitude, 69°). The major part of the 
rocket instrumentation consisted of electron density experiments designed to 
study the state of the lower ionosphere in regions where ionospheric absorp- 
tion takes place. The electron density experiments were carried out by the 
Defence Research Board Telecommunications Establishment. 

The purpose of the cosmic ray experiment was the following: to measure 
particle intensity changes associated with ionospheric absorption events, to 
determine the primary intensity at a time of high solar activity, to look for 
trapped radiation at altitudes around 250 km and, if possible, to obtain an 
estimate of the intensity of low-energy albedo radiation. The detectors em- 
ployed were three Geiger counters (Nuclear Instruments D 79) surrounded by 
various amounts of shielding material. Table I lists the amount of shielding 


TABLE | 
Equivalent Minimum proton Intensity 
shielding, momentum, particles 
Counter No. g/cm? Mev/c cm”? sec”! sterad™! 

l 0.45 Al 182 0.2614 .005 
2 1.8 Cu 257 0.214+ .005 
3 7.6 Pb 350 0.212+ .005 
Telescope 1.9 Cu 260 0. 226+ 014 


(including the nose-cone skin) surrounding the counters and the minimum 
proton momentum to which the counters were sensitive. The omnidirectional 
geometric factor for the individual counters was 7.0 cm? for protons with 
momenta greater than the values listed in Table I. Counters 1 and 2 were 
also operated as a counter telescope inclined 45° to the rocket axis and with 
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a solid angle area factor of 3.1 cm? sterad. Transistorized circuitry fed both 
scaled and unscaled counting rates to subcarrier oscillators and these rates 
were monitored with the aid of telemetry equipment which was part of the 
ionosphere experiment. 

The first rocket was fired on 17 September, 1959. Throughout the flight the 
condition of the ionosphere remained undisturbed and all of the instrumenta- 
tion performed satisfactorily. The second rocket was fired on 20 September, 
1959, after the onset of an auroral absorption event. Unfortunately, the 
rocket failed shortly after firing and no useful information was obtained 
from this flight. 


RESULTS AND CONCLUSIONS 
The counting rates of the four detectors flown on 17 September are shown 
as a function of time and altitude in Fig. 1; the altitude measurements below 
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Fic. 1. Counting rates vs. time of flight and altitude. The points represent counts during 
40 seconds except the points connected by the broken line which are 10-second averages 
normalized to counts per 40 seconds. 


160 km were obtained from radar tracking data and above this the height 
was estimated by extrapolating the trajectory measured below 160 km. Above 
the Pfotzer maximum the counting rates of counters 2 and 8 and the 
telescope are essentially independent of altitude while counter | increased 
by about 25° between altitudes of 50 and 250 kim. The particles responsible 
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for the increased rate in the counter with least absorber are of quite low 
energy. For instance if they are electrons their energy must be less than 3 
Mev or they would have been detected in counter 2 with 1.8 g/cm? of absorber. 
The low energy of the particles and the fact that their intensity increases 
with altitude suggests that they are associated with the trapped radiation in 
the outer Van Allen layer (Van Allen and Frank 1959a). This would be con- 
sistent with their being electrons. 

The average intensities above 60 km calculated from the measured rates 
assuming an isotropic distribution over one hemisphere and zero over the 
other are given in Table I. Apart from the low-energy altitude-dependent 
component already mentioned the intensities as measured by all detectors 
are essentially the same. This implies that there is no measurable contributions 
from protons with momenta less than 350 Mev/c, a value considerably 
above the Stormer cutoff momentum (110 Mev/c) at Fort Churchill. 
The absence of low-energy protons in the primary spectrum at a time near a 
solar maximum would be expected (Neher 1956) but if, as has been assumed 
by Neher, albedo radiation contributes 50% of the total intensity then 
the fraction of these particles with momenta less than 350 Mev/c should 
be detectable in the present experiment. The observed absence of low-energy 
particles can be used along with a reasonable assumption about the albedo 
energy spectrum to give an approximate upper limit of 15% for albedo relative 
to total intensity. Since albedo radiation is concentrated in directions near 
the horizontal the telescope detection efficiency is considerably less than the 
single counter efficiency for this radiation. Hence a further indication that 
the albedo intensity is small is given by the near equality of the telescope 
and single counter results in Table I. This agreement also implies the absence 
of an appreciable contribution by photons. 

In Fig. 2 the intensity determined by counters 2 and 3 (with no albedo 
correction) is plotted along with other high-latitude measurements of the 
primary intensity obtained near the solar maximum; for these measurements 
no albedo correction was necessary. The Zurich sunspot number averaged 
over 12-month periods is also shown. The threshold energies of the detectors 
used in the various experiments differed considerably but all were less than 
the primary spectrum cutoff energy (at solar maximum) and hence a direct 
comparison is possible. Points labelled 2 and 3 are the sum of proton and 
a-particle intensities obtained by Meyer (1959) and McDonald (1959) with a 
combination Cerenkov-scintillation detector. Points 4 and 5 are interplanetary 
values of the primary intensity measured in Pioneer IV (Van Allen and Frank 
19596) and in the Soviet cosmic rocket (Vernov ef al. 1959). The point labelled 
1 was calculated by Neher (1959) from ionization chamber data near the top 
of the atmosphere. With the exception of the ionization chamber intensity 
there is fair agreement among the various values. A rough indication is also 
given of the expected inverse relationship between primary intensity and solar 
activity even though the fraction of a solar cycle covered in Fig. 2 is rather 
small. 

It is interesting to compare the present measurements with Geiger counter 
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Fic. 2. Measurements of primary intensity near a solar maximum. 1, Neher (1959); 2, 
Meyer (1959); 3, McDonald (1959); 4, Van Allen and Frank (1959); 5, Vernov et al. (1959); 
6, present experiment. The points due to Meyer (1959) refer to an atmospheric depth of 
13.5 g/cm? and would be slightly higher if corrected to the top of the atmosphere. 


measurements made at a latitude of 87° near a solar minimum by Meredith 
et al. (1955). It is found that the ratio of the measured intensities in 1952 and 
1959 as determined by Geiger counters is 2.3. This ratio should be very nearly 
independent of albedo effects and should therefore represent the ratio of the 
primary intensities. Neher (1959) using ionization chamber data has given a 
value of 5 for the ratio of primary intensities in 1954 and 1958. 


We are indebted to the Defence Research Telecommunications Establish- 
ment for the opportunity of carrying out the measurements described here. 
In particular, we wish to thank Dr. W. Heikkila and Mr. S. Penstone for their 
very considerable assistance. 
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SOME OXYGEN IONS FORMED AT HIGH PRESSURES IN A 
MASS SPECTROMETER! 


Wn. McGowan AND LARKIN KERWIN 


ABSTRACT 


The oxygen ions formed by electron bombardment and collision processes 
in the mass spectrometer at pressures of about 10-’ mm Hg are examined. Mass 
spectrum lines due to Ot, Of, Of, and O** as well as Aston bands due to O*, 
Of, and Of* are found. Processes leading to these species, as well as abundances 
and appearance potentials, are considered. 


1, INTRODUCTION 
When ordinary molecular oxygen at a pressure of about 10- mm Hg is 
bombarded by electrons of about 50 ev, the resulting ions form a mass spectrum 
such as that shown in Fig. 1. 







L00 OXYGEN 
+ ELECTRON ENERGY 
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Fic. 1. Low-pressure mass spectrum of Oz. 


In this spectrum, the peak at 16 may be identified as consisting of O* and 
Of*, from the reactions 
OO -- O'-F Ze (A) 


O.2+e 


and 

O; +e = Of* + Be. (B) 
The small peaks at 17 and 18 are isotopic variations of reaction (A) involving 
O" and O'8, No other stable isotopes exist to more than 1 part in 20,000 of the 
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whole (Kerwin, McElcheran, and Cottin 1957). Isotopic variations of reaction 
(B) which would appear at 16.5 and 17 are too rare to be identified in the usual 
analysis. 
The peak at 32 consists of Of from the reaction 
O2. +e = Of + 2e (C) 


while the small peaks at 33 and 34 are isotopic variations of (C) involving 
O' and O'8, Isotopic variations at 35 and 36 are normally too rare to be 
detected. 

These reactions are well known (e.g. Hagstrum and Tate 1941; Frost and 
McDowell 1958). The abundances, appearance potentials, etc. are well 
catalogued. In the remainder of this paper, we shall omit reference to isotopic 
variations, although some of these have been observed in the reactions to be 
discussed. At high-bombarding electron energies, a small peak may be seen 
at 8, as has been noted by Mattauch and Lichtblau (1939) and Hagstrum and 
Tate (1941). It is O+* from the reaction 

O22 +e = O** +O + Be. (D) 


2. HIGH-PRESSURE MASS SPECTRUM 

If the pressure of oxygen is raised appreciably above that used for measure- 
ments such as those shown in Fig. 1, the spectrometer filament is usually 
burned. However, a pure carbon filament or a thoriated iridium filament may 
be operated at pressures up to 10-? mm Hg without suffering ill effects. Under 
these conditions, and for a bombarding electron energy of 100 ev, we have 
obtained oxygen spectra such as appears in Fig. 2. In addition to peaks at 
16 and 32 we made the following observations. 
Mass 8: Spur 

The sharp spur sitting on top of the broad peak centered at mass 8 is O++ 
from reaction (D) mentioned above. 
Mass 8: Triangle 

The broad triangle at the base of peak 8 comes from the two reactions: 

O¢ +O; = Ot +0 +0; (E) 


where the O+ produced appears principally at 8, and 

Ot + O2 = Of+ + Or (F) 
where O}* appears also principally at 8. In reactions (EF) and (F), the Of is 
first produced by reaction (C). These broad peaks are known as Aston bands 
and result from ions formed by collisions which take place in or past the ion- 
accelerating region. They then appear to have an effective mass-to-charge 
ratio given to a good approximation by the expression 


bs te EE (; 7 vs) 
dt Mo J dt J 


where mp and m; are the masses of the initial and resulting ions, go and gq; are 
their charges, and V; and V are the voltages of the point within the accelerating 
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iG. 2. High-pressure mass spectrum of Op». 


region where the collision takes place and the total accelerating voltage 
respectively. When the collision takes place past the accelerating region, 
V’, = V, and the equation simplifies to 


(2) Lea (a) 20 
m 1 gd: / mo’ 


which is given by Mattauch and Lichtblau (1939). The charge-invariant form 
of equation (1) is given by Hipple, Fox, and Condon (1946). If collisions take 
place at various points inside the accelerating region, then a spread of masses 
m, will be observed, making up part of the high-mass side of the Aston band. 
The low-mass side is partially attributed to fringing field effects. The precise 
shape of this part of the band does not appear to have been satisfactorily 
explained. 
Mass 48 

A small peak at 48 may be identified with the ozone ion, Of. This may arise 
in a variety of ways, for example, 

O; + O = OF. (G) 


This peak has previously been observed (Grundland 1953; Herron and Schiff 
1958). 
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Mass 64 
This peak could result from two processes: Of has been reported by Grund- 
land (1953), who examined with a mass spectrometer the products of a dis- 
charge in oxygen. There is also the collision process, 
Of* + O2 = OF + OF, (H) 


where, as a result of being formed after the Of* has been accelerated in the 
spectrometer, the Of appears as an Aston band at mass 64. The Of* is first 
formed from reaction (B). As will be seen, our peak is not Of. 


3. DISCUSSION 
The reactions (D), (E), (F), and (H) postulated in Section 2 are those 
considered in this paper. Their assumption is based on the following con- 
siderations. 


Mass 8: Spur 

The mass spectrometer is fitted with an energy filter at the detecting end. 
When this is set to pass ions of any energy and when bombarding electrons of 
100 ev are used, the peak at mass 8 appears as in Fig. 2, and is reproduced in 
Fig. 3a. When the ion energy selector is set to repel ions not having essentially 
all of the initial ion beam energy, the triangle disappears, leaving the spur 
(Fig. 3b). This is a sharp peak with practically no energy spread. Its relative 








Fic. 38. Appearance of ion peak at mass 8: (a) low-energy barrier, high electron energy; 
(b) high-energy barrier, high electron energy; (c) low-energy barrier, low electron energy. 


intensity of 0.1% of the Of peak does not vary with pressure, indicating a 
primary process, as shown in reaction (D). Hagstrum and Tate (1941) found 
an abundance of 0.01% for this peak. The abundance does not vary with ion 
energy of from 400 ev to 2000 ev. 

An appearance potential curve made for the spur (Fig. 4) shows that the 
peak first appears at 64.9 ev, and that there are several well-defined breaks in 
the curve between this value and 80 ev. Now the total energy required to 
produce O++ according to reaction (D), if everything is in the ground state, is 
53.8 ev (e.g. Herzberg 1944). While lack of sensitivity may have failed to 
disclose the onset of the reaction at this point, it is clear that an abundant 
reaction sets in at 64.9 ev and that several others are involved at higher 
energies. 
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Fic. 4. Appearance potential curve for O*t. 


A general reaction to cover this experimental fact would be 
Oz +e = Ott* + O* + Be. (D*) 


Here the starred species are in excited states. Experimental conditions preclude 
the possibility of the O2 being in an excited state. There are several combina- 
tions of excited states for O++, and O which satisfy (D*) (see Moore (1949) for 
values). However, the experimental errors (about 0.4 ev) in our appearance 
potential curves resulting from the use of a carbon filament do not permit us 
to select a combination with certainty. The following table gives one which 


TABLE I 
Suggested excited states for reaction (D*) 












Appearance potential for (D*) 





Observed 





O*+* state O state Reference 






















2,° *P 3p *P 64.9 ev 64.8 ev 
2,7 2D 3p 3P (Not observed) 67.3 
2,7 1S 3p °P 69.4 70.2 
2,9 8S 3p °P 72.3 72.3 
2,° 3D 3p *P 79.0 79.7 









has the virtue of postulating the O atom in a single state 3p °P and the virtue 
of being relatively simple. The observed break at 74.5 ev is not accomodated 
by this scheme. A number of other schemes do fit it, and it may indicate a 
multiple process. 
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Fic. 5. Appearance potential curves for reactions (E) and (F) and for Of. 
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Mass 8: Triangle 

When the ion beam filter is adjusted to pass low-energy ions, and the electron 
energy is reduced to about 50 ev, the spur disappears, leaving the triangle. 
Thus the latter is formed by ions with a lower appearance potential, and 
having only a fraction of the initial ion beam energy. The spread of the peak 
is characteristic of Aston bands, as shown in Fig. 3c. 

Appearance potential curves were obtained from the high-mass side, center, 
and low-mass side of the triangle, all being identical. An average curve is 
shown in Fig. 5, where it is compared with an appearance potential curve 
for Of. 

At low energies, breaks were found in the Of curve at 12.2 ev, 16.1 ev, and 
18.2 ev (Fig. 4) corresponding to the ground level, 427 and 47, states of Of 
respectively (Herzberg 1950). These breaks have previously been observed by 
Herron and Schiff (1958), and Frost and McDowell (1958). These are not 
drawn in Fig. 5, but served to calibrate the energy scale. At higher energies, a 
series of breaks in the Of curve which have not been previously reported are 
shown. The correspondence of breaks in the triangle curve leaves little doubt 
that the Of is the source of the latter, as indicated by reactions (E) and (F). 
The Ot and Of* ions forming the triangle according to these reactions would 
then result from the collision of an excited Of ion with an O: molecule. 

The relative abundance of the triangle varies with pressure, as shown in 
Fig. 6. The initial linear part of the curve confirms that the process is a secon- 
dary one (e.g. reactions (E) and (F)). At higher pressures the Of beam is 
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Fic. 6. Variation of abundance of mass 8 triangle with pressure. 
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scattered. In appropriate experiments where beam intensity was measured as a 
function of pressure, the Of maximum occurred at considerably lower pressure 
than that for mass 8. In a series of experiments to be reported later, the cross 
section for the sum of reactions (E) and (F) was found to be about 2X 10-7 cm?. 

An analysis of the energy spread of the triangle ion beam was made using 
the energy filter. In Fig. 7 may be seen the successive appearances of the total 
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Fic. 7. Variation of appearance of mass 8 peak with barrier potential. 


ion beam at mass 8 as the pass-level of the filter is lowered. Note that the spur 
component and the mass 12 peak remain unchanged. The total accelerating 
potential was 1000 volts. The high-energy component of the triangle is seen to 
appear first. The main energy level of the triangle is about one-half the main 
ion beam energy, confirming the postulated reactions (E) and (F). A small 
energy loss to the reaction has been consistently noted. 


Mass 48: Ozone 

The abundance of this ion was 0.05% of the Of peak at source pressure of 
2X10-* mm Hg, and using 100-ev bombarding electrons. No quantitative 
measurements other than this could be made on this small peak. At sensibly 
lower pressures the peak was negligible. 


Mass 64 

In this charge exchange reaction (H) the resulting Of ion appears at mass 64, 
and should have twice the energy of the main ion beam (since the parent Of*+ 
was accelerated with two charges). This was confirmed with the energy filter, 
the beam (0.05% of Of at 1.5X10~ mm Hg) persisting after all others had 
been cut off until almost twice the main beam energy was applied. This 
eliminates the possibility of the peak being Of in our case. An Of abundance 
of less than 0.005% of Of would have been measurable. The cross section for 
this charge exchange reaction was established through reference to the peak 
at 16. If this is considered to be at least 90% Ot, then the cross section for 
reaction (H) is 1X10~'* cm*. It is to be expected that the contribution of Of+ 
to mass 16 is less than 10%, in which case the cross section for (H) will be 
larger. No precise measurements of the variation of this small peak with 
pressure were made. 
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An appearance potential curve measured for mass 64 (with 1800 volts 
potential barrier applied) gave an initial break at 46 + 1.5 ev and a further 
well-defined break at 58 + 1.5 ev. The reference value for the initial break 
(Herzberg 1944) is 47.2 ev. Now in Fig. 5, the various breaks in the appearance 
potential curve for the collision processes were seen to fit closely the breaks in 
the appearance potential curve of the parent Of. The same correspondence 
should attain between the curve for mass 64 and the curve for the parent Of*. 
A thorough search for such breaks at mass 16 revealed none, indicating that 
the Of+ was a very minor component of the peak at 16. Hagstrum and Tate 
(1941) on the other hand found a break in the mass 16 appearance potential 
curve at 50 + 0.15 ev. The possible reaction 

O+ + O; = 20+ + O,, ) 


whose products appear at mass 16, might be a complicating factor. Further 
data on the portion of the mass 16 peak, which may be ascribed to Of+ as 
discussed in the next section, leads to the conclusion that the cross section for 
reaction (H) is greater than 2X 10~'§ cm’. This surprisingly large value will be 
discussed in a subsequent paper dealing with cross sections. 


4. COMPLETE APPEARANCE POTENTIAL CURVES 


In Fig. 8 are shown appearance potential curves for the species O*+, Of, 
O+ + Of, and Of*. Salient characteristics of these curves are gathered in 
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TABLE II 
Appearance potential characteristics from Fig. 8 


Plot for amas, 95% width, omex relative 


Ion volts volts to OF 
OF 100 45 100 
Ot + Of* 105 85 40 
Ot+ 185 95 0.1 
O;+* 115 35 2 


From Fig. 8, the 95% widths for the diatomic ions are seen to be about half 
the widths for the monoatomic ions. The curve labelled Of* only is taken from 
data on the mass 64 peak, at which point the products of the collision-induced 
dissociation appear. A similar curve would be found at mass 16, if it were 
possible to eliminate the O* ion. The actual curve measured at mass 16 is then 
a combination of those for OF and this Of*. From the curves in Fig. 8 it may 
be estimated that the Of* contributes less than 1/20 of the mass 16 curve, 
the subtraction of greater amounts than this from the composite curve leaving 
a very irregular curve which does not correspond to those usually observed. 
It is hoped to be able to initiate a search for negative ions which might throw 
further light on some of these processes. 
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MEAN LIVES OF POSITRONS IN ALUMINUM AND THE 
ALKALI METALS'! 


R. E. BELL? AND M. H. JORGENSEN 


ABSTRACT 


The time distribution of positron annihilations in the metals aluminum, 
lithium, sodium, potassium, and cesium have been measured with a fast time-to- 
amplitude converter. The decay curves appear to be complex, with about 5% of 
the events having a mean life of approximately 5X107!° sec. The main (95%) 
components of the decay curves show the following mean lives, in units of 107!° 
sec: Al, 1.9+0.2; Li, 2.9+0.2; Na, 3.15+0.2; K, 4.0+0.2; Cs, 4.3+0.2. The 
results for the alkali metals disagree with the previously published measurements 
of De Benedetti and Richings (1952). These results are discussed, and the lifetimes 
for other metals are predicted roughly from them and from the angular correlation 
measurements of other workers. 


INTRODUCTION 

This paper reports new measurements of the lifetimes of positrons in alu- 
minum and the alkali metals (except rubidium). The measurements were 
made by the delayed coincidence method using apparatus whose time resolution 
is good enough to allow the lifetime to be read directly from the logarithmic 
slopes of the resolution curves. All previous results for positron lifetimes in 
metals have been measured by the centroid-shift technique introduced by 
Bay (1950). This method, though statistically excellent, is subject to syste- 
matic errors that can completely upset the results. As an example, the total 
transit time of the pulse through a typical end-window photomultiplier (e.g. 
RCA 6342A) is about 3X 10~-° sec; a drift of only one part in 10° in this transit 
time causes a spurious centroid shift of 3X10-" sec. Similar effects can be 
caused by changes in pulse size distribution and counting rate. On the other 
hand the logarithmic slope of the resolution curves is insensitive to small 
changes of this kind. 

The most quoted measurements of the mean life of positrons in aluminum 
are those of Bell and Graham (1953) and Gerholm (1956). Both results were 
obtained by the centroid-shift method, but with different experimental 
arrangements; they were (1.5+0.3)X10-" sec and (2.5+0.3)X10-" sec 
respectively. Other measurements have been made by Ferguson and Lewis 
(1953: 7 = (1.640.6) X10-" sec), Minton (1954: 7 = (2.9+0.3) X10-" sec), 
and Jones and Warren (1956: 7 = (1.740.3) X10-" sec). A measurement of 
the mean life of positrons in lead by Sunyar (1957) gave the value (1.6+0.15) X 
10-" sec, presumably supporting the shorter mean life values for aluminum. 
The present experiment gives a mean life lying among all these values, but 
somewhat favoring the shorter mean lives. 

Measurements of positron mean lives in the alkali metals have been reported 
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only by De Benedetti and Richings (1952). They measured the difference 
between the mean life in aluminum and the mean lives in lithium, sodium, 
and potassium, and found in every case that this difference was zero, with an 
experimental error of 0.7X10~-" sec. The results of the present measurements 
are in strong disagreement with this finding; for example the mean life in 
potassium exceeds that in aluminum by three times the assigned error of 
De Benedetti and Richings. An unpublished measurement of the mean life of 
positrons in cesium by Madansky is quoted as a footnote of a paper by Ferrell 
(1956). Madansky’s value, (3.7+0.6) X10-" sec, agrees within errors with 
the value found here. 
EXPERIMENTAL PROCEDURE 

(a) Apparatus 

The measurements were made in the customary way with sources of Na™ 
immersed in the metal under test. Two diphenylacetylene gamma counters 
(2-cm diameter X 1.5 cm high), viewed by two RCA-6342A photomultipliers, 
were placed near the source, 90° apart. Of these counters, No. 1 was biased 
to accept the upper one-third of the Compton electron distribution of the 
1.28 Mev nuclear gamma ray accompanying each positron from Na*, and 
No. 2 was biased to accept the upper one-third of the Compton electron 
distribution of the annihilation radiation. The time distribution of positron 
annihilations in a metal was measured by analyzing the distribution of delays 
of counter 2 with respect to counter 1, using a fast time-to-amplitude con- 
verter and a 100-channel pulse analyzer. The instrumental time resolution was 
determined by replacing the positron source by a source of Co™, without 
change in the counter biases. Ideally, such a pair of measurements yields the 
mean life of the positrons in the metal directly from the shift between the 
centroids of the two curves. In practice the errors already mentioned may make 
such a measurement untrustworthy, and we use the prompt Co” curve only to 
demonstrate the time resolution of the apparatus. 

The time-to-amplitude converter is similar to that described by Green and 
Bell (1958), but with a number of alterations. The actual conversion from time 
to amplitude in the apparatus of Green and Bell was made in a 6BN6 gated 
beam tube. We have found that several different devices will serve about 
equally well for this purpose; at present we are using a biased diode followed 
by an integrating condenser, as used earlier by Deutsch (unpublished) and 
Sunyar (1957). The choice was made mainly on the basis of convenience and 
simplicity. The instrumental time resolution curve for the conditions of the 
present experiment has a full width at half maximum of 5X10~-" sec, and 
logarithmic slopes of the sides of the curve corresponding to a half life of 
6X10-'! sec (mean life < 107" sec). Since the shortest mean life found in 
these measurements was 1.9X10-! sec, it was always justifiable to read the 
mean life directly from the logarithmic slope of the measured time distribution 
curve. 

The apparatus included a pile-up detector that rejected any coincidence 
event that was preceded or followed within 3 usec by another pulse in either 
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counter. The pile-up detector had a resolving time of ~30 musec, that is, it 
recognized two pulses 30 musec apart as separate. Except for the loss of a few 
genuine events, the effect is the same as if all the amplifiers and pulse height 
selectors had resolving times of 30 musec. 

The time-to-amplitude converter shows a residual effect of the pulse heights 
from the individual counters on the time analogue pulse. In the present 
apparatus this effect is cancelled to first order by adding a fraction of each of 
the individual counter pulses to the time analogue pulse; the fraction is adjust- 
able in magnitude and sign. Such a device, correctly adjusted, widens the 
range of acceptable pulse heights from the individual counters without de- 
grading the resolving time. The final effect is an increase in coincidence 
counting efficiency. 

(b) Preparation of Samples 

The samples of aluminum were prepared by depositing some drops of Na™ 
solution (as chloride) on a strip of household aluminum foil, and evaporating 
them to dryness. The strip of foil was then tightly rolled into a cylinder thick 
enough to prevent the escape of positrons. Household foil was chosen because 
it has been carefully degreased for domestic reasons. The source strength was 
usually chosen to be about 5 wc. With Na” of reasonably high specific activity, 
positrons will annihilate only in aluminum, in such a sample. 

The alkali metal samples were more difficult to prepare. In the case of 
lithium, the least reactive metal, it was sufficient simply to take two flat blocks 
of lithium metal, scrape their surfaces clean, and clamp them together over a 
Na” source deposited on 2.5 mg/cm? aluminum foil. The whole assembly was 
then dipped in molten paraffin wax, in order to exclude air. In the cases of 
sodium and potassium, the Na” activity was deposited on a 2.5 mg/cm? 
aluminum foil supported on a ring of stiff copper wire. This ring was held 
vertically in the bottom of a glass test tube just large enough to contain it, and 
the molten alkali metal was poured around it through a narrow heated funnel 
that retained the oxide layer. The alkali metal lumps that were melted for the 
purpose had been carefully scraped clean beforehand. All the dimensions were 
so chosen that no positrons could reach either the wire loop or the glass walls 
of the test tube. Molten paraffin was poured in on top of the alkali metal to 
seal it. The procedure for cesium metal was similar, but the whole process had 
to be carried out in a glove box filled with purified nitrogen. In addition, the 
melting point of cesium is lower than that of paraffin wax, and the seal could 
not be made in the usual way. Instead, a lucite test tube was used in place of 
the glass one, and a lucite stopper was cemented into it with acetone. 

All the alkali metals, seen through the transparent walls of their containers, 
were perfectly bright and metallic in appearance, and remained so for periods 
of months. In all the samples, a few positrons, perhaps as many as 5% of the 
total, will annihilate in the aluminum support foil. The effect on the results 
will be very slight, and no correction has been made. 


(c) Measuring Procedure 
The measurement of a lifetime consisted of a run of some hours on the time- 
to-amplitude converter, using the sample in question, often immediately 
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preceded or followed by a shorter run made with a source of Co® replacing the 
positron sample. The counting rates from the pulse height selectors were a few 
hundred counts per second, and the coincidence rate was around 200 counts 
per minute. The Co® comparison source was selected to give similar single 
counting rates. The scale of the time-to-amplitude converter could be altered 
by changing the over-all gain of the system. It was calibrated by delaying the 
pulses from one counter with different lengths of RG 114/U coaxial cable 
(Zo = 185 Q, v/c = 0.88). This calibration was always made by observing the 
shift of the centroid of the actual time distribution being measured, and it 
should be accurate within 2 or 3%. (The centroid shifts used for calibration 
are large, and the centroid-shift errors mentioned above are then negligible.) 
The calibration was stable within about 1% over periods of many days. The 
position of the centroid of a time distribution sometimes showed drifts of a few 
times 10~-" sec, or about one channel of the 100-channel analyzer. Whether 
these drifts were occurring in the time-to-amplitude apparatus or in the 
analyzer itself was never determined. Such drifts have almost no effect on our 
results, although they would be very serious in the centroid-shift method. 

The chance coincidence rate was observed at large values of the delay, and 
was also calculated from the measured single counting rates. The two values 
were always consistent within statistical expectation. The chance rate to be 
subtracted was usually about 3X10-‘ of the peak coincidence rate. All the 
conclusions of this experiment can be drawn from points having counting 
rates at least three times the chance rate. 

In regions where the counting rate is low, the observed data have often been 
plotted with adjacent channels of the analyzer averaged to give a single 
plotted point. This step simplifies the appearance of the curves without essential 
change in their content, and also serves in many cases to remove zero values, 
which cannot be plotted on a logarithmic scale. 

RESULTS 

One of many measured time distributions for positrons in aluminum is 
shown in Fig. 1, with a Co prompt curve plotted to the same area. The fact 
that positrons in aluminum have an easily measurable mean life is clear in 
Fig. 1, but the decaying part of the curve (right-hand side) is not fitted by a 
single exponential. All our measured curves on aluminum and other metals 
show this “‘tail” effect in comparable degree; no such effect has ever been 
observed with an ordinary radioactive source, except in those rare cases where 
an actual mixture of two lifetimes was present. 

The method of treating the data of Fig. 1 is to draw a reasonable straight 
line through the points of the tail, and subtract the component so found from 
the main curve. The line used for subtraction in Fig. 1 represents a component 
of mean life 3.85X10-" sec with an intensity of 6.5%. These values are of 
course extremely rough. The remainder curve, representing ~ 94% of all the 
annihilations, then shows a mean life of 1.93 10-" sec. The same procedure, 
applied to a number of other runs with positrons in aluminum, gave values for 
the mean life ranging from 1.85 to 1.95X10~-" sec. We therefore quote the 
value (1.9+0.2) X10-" sec. The assigned error is large enough to cover the 
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Fic. 1. One of several measured time distributions for positrons in aluminum, with a 
prompt curve (Co) for comparison. 


mean life found by omitting the subtraction of the tail, and drawing a line 
through the steepest part of the curve (about 2.0X10-" sec). 

The results for lithium and sodium appear in Figs. 2 and 3. (The lithium 
result shows the effect of a drift of centroid between the recording of the 
prompt curve and the recording of the lithium curve.) Both lithium and 
sodium appear qualitatively similar to aluminum, but the mean lives found 
are (2.9+0.2)X10~-" sec for lithium, and (3.15+0.2)X10-' sec for sodium. 
The tail effect in both lithium and sodium resembles that in aluminum. The 
curves for aluminum, lithium, and sodium were all recorded in one series of 
runs lasting 3 days, on the same scale of the time-to-amplitude converter. By 
the time we reach potassium and cesium, the mean life has lengthened enough 
to require a change of the time scale of the apparatus by a factor of about 1.5. 
Figure 4 shows the result for potassium, again with a Co™ curve for reference. 
The tail effect is now being lost under the longer main lifetime, and by the 
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time we reach cesium (Fig. 5) it is no longer distinguishable. Even for potas- 
sium, the subtraction of the tail makes no difference to the mean life of the 
main component. The mean lives of positrons in potassium and cesium are 
(4.0+0.2) X10—" sec and (4.30.2) X10-" sec respectively. 

We now return to the weak tail of longer lifetime, as seen for example in 
aluminum (Fig. 1). Owing to its lifetime, this tail cannot be due to the accidental 
presence of oil or wax in the sample. Strenuous efforts were made to try to 
eliminate the effect, and thus to demonstrate that it was instrumental, but 
without success. One of these efforts was the installation of the pile-up detector 
already mentioned. This detector should reduce pile-up effects by a factor of 
the order of 20, and in fact it did clean up the left-hand side of the curve 
appreciably, but the effect on the right-hand side of the curve was negligible. 
Another attempt consisted of a measurement on an aluminum sample like the 
first, but having only 1/10 the Na® activity in it, and hence giving 1/10 the 
relative pile-up effect. The result is shown in Fig. 6, and although the statistical 
accuracy is low, the tail effect is still present in the usual amount. This measure- 
ment was immediately followed by one using a source of Tb'® of strength 
equivalent to the Na®* in the aluminum. Tb" has a gamma spectrum which, 
when viewed with an organic phosphor, is very similar to that of a source of 
Na*. In the case of Tb'®*, however, the higher energy radiations are delayed 
with respect to the ~ 500 kev radiations by the half life of the 1507 kev excited 
state of Gd'**, 1.88 X10~'° sec. The delay is thus in the opposite sense to that 
of a Na” source. The result for the Tb"® source is also shown in Fig. 6. The 
exponential decay extending to the left is due to the half life just mentioned. 
(It was, in fact, this observation that prompted the present authors to make 
careful lifetime measurements of excited states of Gd'* (Bell and J ¢rgensen 
1959).) The right-hand side of the Tb'* curve drops rapidly to zero, and remains 
identically zero for 60 consecutive channels of the analyzer, covering the 
entire region where the curve for positrons in aluminum is in question. This 
comparison of two weak sources that are closely identical in everything but 
lifetime strongly suggests that the tail effect seen with positrons in aluminum 
is real. Other attempts to eliminate it, involving different counter voltages, 
different pulse height selection, interchange of scintillation counters, different 
angular positions of counters, and lead shielding to reduce intercounter 
scattering, were all equally unsuccessful. The old aluminum curves of Bell 
and Graham (1953, Figs. 2 and 11), made with much worse time resolution in 
a different kind of apparatus, show signs of the tail. The aluminum curve of 
Green and Bell (1958a) made with somewhat worse time resolution in 
another time-to-amplitude converter, shows an effect that is statistically poor, 
but very similar. Gerholm (1956) found a tail of somewhat longer lifetime in 
aluminum, which he has since attributed at least in part to pile-up effects 
(private communication). We are informed by Dr. Gerholm, however, that he 
has never succeeded in observing a perfectly straight decay curve for positrons 
in aluminum. There is thus a good deal of precedent for finding the extra tail 
in aluminum, and it appears to be similar in the alkali metals. All the experi- 
mental results are collected in Table I. 
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Fic. 6. Comparison of a weak source of positrons in aluminum (right-hand curve) with an 
equal source of Tb'*¢ (left-hand curve). 


TABLE I 


Results of lifetime measurements on aluminum and the alkali metals 


Main component Tail component 


Metal r (10-"° sec) r (107° sec) Intensity (%) 


Al 1.9+0.2 3.85 6.5 
Li 2.9+0.2 5.9 6.5 
Na 3.15+0.2 5.7 4.3 
K 4.0+0.2 — — 
Cs 4.3+0.2 _— — 


DISCUSSION 

The main parts of the decay curves of aluminum and the alkali metals will 
be discussed first. We plot in Fig. 7 the experimental annihilation rates 
\ (= 7"') as a function of the dimensionless parameter 7,. (This parameter is 
the radius, in units of the Bohr radius, of the sphere whose volume equals the 
volume occupied by each valence electron in the metal. It is therefore pro- 
portional to the inverse cube root of the density of valence electrons. The 
valence number is taken as 3 for aluminum and 1 for the alkali metals.) 
Within their assigned errors, the experimental points lie on a smooth curve, 
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although there is no reason to expect this smoothness to be perfect. The 
broken curve at the left in Fig. 7 shows the value of \ calculated on the basis 
of a Sommerfeld free electron gas, using the Dirac annihilation cross section 
without Coulomb effect. A second broken curve in Fig. 7 is from the work of 
Ferrell (1956), in which a rough allowance is made for the enhancement in A 
due to the Coulomb attraction between the positron and the electrons. Ferrell 
expected his curve to be unreliable at low electron densities (7, greater than 
about 4), and he suggested that in this region the annihilation rate should 
approach that of the positronium negative ion, 3.06 X 10° sec~!, shown at the 
right in Fig. 7. In fact the experimental curve approaches more nearly the 
value for the spin-averaged positronium atom, 2.0X 10° sec~, also shown at 
the right in Fig. 7. Nevertheless the experimental results support Ferrell's 
qualitative idea. 


10 
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EXPERIMENTAL 


» (10 °SEC* UNITS) 
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fs 


Fic. 7. Plot of \ (= positron annihilation probability) as a function of the dimensionless 
parameter r,. The experimental points are shown with a solid curve drawn through them. The 
broken curve marked “electron gas’’ is 15% below the experimental curve. For the significance 
of this and the other curves and points, see the text. 


The final curve shown in Fig. 7 is the dotted one labelled “electron gas’, 
and it is derived from the experimental curve in the following way. Aluminum 
and the alkali metals are members of a large class of metals in which the 
measured angular correlation of the annihilation photons near 180° consists 
mainly of an inverted parabola (Green and Stewart 1955; Stewart 1957; Lang, 
De Benedetti, and Smoluchowski 1955; Lang and De Benedetti 1957). The 
parabolas have the shape and width expected for the Fermi momentum dis- 
tribution corresponding to a simple Sommerfeld model of the metals concerned, 
together with a momentum-independent annihilation rate. In each such case 
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the parabola is superimposed on a low distribution extending to higher 
mementa. In the cases of aluminum and the alkali metals, this higher- 
momentum tail contains about 15% of the events. The tail is attributed by 
De Benedetti, Cowan, Konneker, and Primakoff (1950), and by Ferrell (1956), 
to the excluded-volume effect of the metallic ions, but Berko and Plaskett 
(1958) argue strongly that it is mainly due to annihilation of positrons with 
the outer electrons of the metallic ions (‘‘core annihilation’’). We take the 
view that it is an effect in addition to the annihilations of positrons with the 
electron gas. In that case an electron gas alone would have an annihilation 
rate about 15% smaller than that of the experimental curve in Fig. 7, as 
indicated by the dotted curve so labelled. This electron gas curve approaches 
the positronium atom value of A very closely at large 7,. 

A first goal of a theoretical treatment of positron annihilation in metals 
presumably is to produce a curve close to the electron gas curve. The details 
that distinguish a real metal could be added later. Ferrell’s enhanced curve is 
a step in this direction; he also made a less successful calculation by another 
method. The experimental lifetimes show that the enhancement of the electron 
density at the position of the positron is very large, but the degree of enhance- 
ment must be nearly independent of electron momentum in order not to spoil 
agreement with the angular results. This remark eliminates the simple 
Coulomb enhancement factor appropriate to a heavy positron and a single 


electron, 
(2rac/v)/{1—exp (—22ac/v) } 


(Heitler 1944), because it introduces a factor 1/v (v = electron velocity) into 
the annihilation rate, at low velocities. (Nevertheless, a detailed computation 
using this factor was made by Karpman and Fisher (1956).) Ferrell’s discussion 
suggests that a reasonably complete theory will not have difficulty on this 
point. In a recent paper, Kahana (1960) has attacked the problem on the 
basis of the Bethe-Goldstone equation (1957), first developed for nucleons. 
Both the Pauli principle and electron screening are included, but the computa- 
tions contain a number of approximations. Kahana’s three computed points 
are shown as open squares in Fig. 7. They display the correct sort of variation 
of X with r,, but lie about a factor of 2 too high. It is not clear whether this 
discrepancy is due to the approximations used, or to the model. Kahana also 
shows that the annihilation rate varies very little with electron momentum, as 
required in a correct theory. 

Turning to the weak component of longer lifetime, we may summarize the 
results by saying that about 5% of the positrons have a mean life of about 
5X 10-" sec in the metals where the effect is distinguishable. These properties 
are consistent with what would be observed if about 5% of the positrons 
formed positronium atoms in the metal, and these positronium atoms under- 


went rapid triplet < singlet exchange. In that case a mean life a little less 
than 5X10-" sec would be observed (Garwin 1953; Green and Bell 1958qa). 
Moreover, it has recently been shown by Bég and Stehle (1959) that slow 
positrons in an electron gas have a capture rate for electrons of the order of 























663 





BELL AND JORGENSEN: MEAN LIVES OF POSITRONS 


10 sec—!, the mechanism being the modified Auger effect first mentioned by 
Garwin (1953). Unfortunately this explanation is contradicted by the results 
on angular correlation of annihilation quanta. In amorphous materials, where 
positronium is formed, a narrow component appears in the angular distribution, 
corresponding to annihilation of positrons in slow-moving positronium atoms 
with the electrons to which they are bound (Stewart 1955; Page et al. 1955). 
This effect has been carefully sought, for a different reason, in the case of 
aluminum by Berko and Plaskett (1958), without result. Judging by their 
published curves, they would have detected as little as 1/10 as much narrow 
component as we have been discussing. If positronium atoms in a metal are 
thermalized or nearly so, the absence of the narrow component seems con- 
clusive, and the 5% tail remains unexplained. Further work is needed. 

The statement of Bell and Graham (1953) that the mean lives of positrons 
in metals appeared to be independent of the metal chosen within +25% was 
based on the work of De Benedetti and Richings (1952), so far as the alkali 
metals were concerned. The general validity of the statement is now destroyed, 
although it still presumably applies to most metals. The electron gas curve of 
Fig. 7, and the angular correlation results, enable us to predict roughly the 
experimental mean lives to be expected for other metals. We calculate r, for 
each metal from its valence, and read its electron gas annihilation rate, or 
inverse mean life 7,~!, from Fig. 7. We then turn to the angular correlation 
results, and read out roughly the fraction f of all the annihilations that lie in 
the parabolic part. This quantity can be found for some metals from Lang and 
De Benedetti (1957, their Table II), and for others by sketching and measuring 
areas on the curves of Stewart (1957). The predicted experimental mean life 
is frg. Some results of this procedure are shown in Table II. The metals in this 
table have been selected to show a wide variation in 7, (i.e. in electron density), 
so that 7, in the table varies from 2.04 to 3.62 10~" sec. On the other hand 
frg varies only from 1.43 to 2.17X10-" sec, lying in the correct range of 
absolute values, and having a reduced spread that easily fits within the 
+25°% range suggested by Bell and Graham (1953). Many metals (notably 
the transition metals and the metals around gold and platinum) have such 


TABLE II 


!Rough prediction of the experimental mean life fr, to” be expected 
in several metals, on the basis of the mean life 7, from the ‘‘elec- 
tron gas” curve of Fig. 7, and the fraction f derived from the 

angular correlation results 








Metal Te tT, (107° sec) ¥ fre (10~"* sec) 
Be 1.88 2.04 0.94 1.92 
C (graphite) 1.92 2.08 0.74 1.54 
Cu 2.67 2.70 0.53 1.43 
Zn 2.30 2.39 0.68 1.62 
Ag 3.01 3.00 0.58 1.74 
Ba 3.74 3.62 0.60 2.17 
Hg 3.40 3.33 0.47 1.57 
Tl 3.58 3.47 0.53 1.84 
Pb 2.30 2.39 0.75 1.79 
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complex angular distribution curves that it is scarcely possible to evaluate f 
for them at all. Nevertheless a theoretical calculation of an electron gas curve 
resembling that in Fig. 7, taken together with the existing discussions of the 
angular correlation data (especially Ferrell 1956; Berko and Plaskett 1958), 
would give us a good initial understanding of positron annihilation in metals. 


We wish to thank Professor N. Bohr for the excellent working conditions in 
his Institute, and Professor A. Bohr and ing. O. B. Nielsen for their helpful 
interest. One of us (R. E. B.) is grateful for the hospitality of the Institute 
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AN HYPOTHESIS CONCERNING TURBULENT DIFFUSION! 


R. BoURRET 


ABSTRACT 


It is shown that the Goldstein equation for turbulent diffusion implies diffusion 
currents dependent upon the history of the concentration gradient. An analysis of 
the stochastic model upon which the Goldstein equation is based reveals that the 
hereditary function, by which the history is weighted, is the ensemble autocor- 
relation function of velocity. Heuristic arguments and an appeal to the theory of 
irreversible thermodynamic processes lead to the postulation of an integro- 
differential equation for turbulent diffusion involving the velocity autocorrelations 
of the diffusate particles. 


I. THE GOLDSTEIN EQUATION 


In order to consider a proposed extension of the theory of diffusion as 
developed by Taylor (1921) and Goldstein (1951) and extended by Michelson 
(1954) and Davies et al. (1954), a brief review of Goldstein’s analysis is pre- 
sented first. His notation has been altered slightly. 

Consider a one-dimensional random flight approximated, as in Fig. 1, with 
discrete intervals of space n and of time r+. If the probabilities of occupation 


j+l 


Fic. 1. Diagram of diffusive movements described by Goldstein’s analysis. 


of all the sites are known at an initial time, what are these probabilities at a 
later time? This is the problem to be solved by specific assumptions about the 
site-to-site transition probabilities. In these transitions, it is assumed that the 
diffusing particle moves with a constant velocity, v = 7/7, either to the right 
or to the left. The probability that the 7th site will be occupied at time k by a 
leftward- or a rightward-moving particle is designated by Ly, and Ry, respec- 
tively. The total probability for the occupation of the site is then 


(1) Px = sit Ru. 


1Manuscript received December 29, 1959. , ; ; ; - 
Contribution from the Research Laboratories, Hughes Aircraft Company, Culver City, 
California. An abstract of this paper was given at the International Oceanographic Congress 


held in New York, 31 August-11 September 1959. 
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If p, and p_are the respective probabilities for unaltered and reversed motion 
in unit time, then it can be inferred from Fig. 1 that 


(2a) Rit, 1 = P+RatpLeu 
and similarly that 
(26) Lin, en = PxLutp_Ru. 


These difference equations become differential equations when we explicitly 
introduce the time interval 7 and the space interval n, expand to first order, 
and take the limit 7 — 0. Thus, from (2) it is seen that 


R(x+n, t+r) = p+(7)R(x, t) +p_(r)L(q, 2) 
L(x—n, t+r) = p_(r)R(x, th +p4(r)L (x, t). 


Performing the expansion, dividing by 7, and taking the limit yields 


(3) 


O0R, OR 
(4a) ate = RL 
and 

aL 06 
(46) es = r_R+riL. 


It has been stated that v = lim (/r), and 74 and r_ have been defined as 
follows: 750 


(5) r, = lim [ #2@)—1 | ; r = lim Eg , 


730 730 


By the use of (1), the following equation can be obtained for P alone: 


aP 10P 10°P 
(6) a Da vor” 


where the “diffusivity” is defined by 


7) Din es 


the latter equality holds by virtue of the definitions (5). This modified dif- 
fusion equation, originally stated by Goldstein, has a very interesting feature 
to which attention seems not to have been called. Using the equation of 
continuity with J representing the diffusive flux, i.e. 


oP 
(8) att? = 0, 


the three-dimensional form* of (6) may be expressed as 


_ *This is the three-dimensional form given by naive induction from (6). The strict generaliza- 
tion to more than one dimension is somewhat involved. The case of two and three dimensions 
has been developed by the author and will be described elsewhere. 
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1 1 aJ 
(9) v (vps I+3 2») = 
and hence it can be inferred that 
a3 ,v 2 
(10) a tDs = —v VP+H(t)+V XG(r, ¢). 


Since the functions G(r, ¢) and H(t) have no physical significance in the 
present inquiry, they are neglected here. It is clear that in the steady state 
(10) reduces to Fick’s law. The general solution, however, may be written in 
the form 


(11) I(r, t) = —DVPrer(r, t) +5 (4, Oe", 
where 
D 
T= 3 
Uv 
and 
ot 1 - 
(12) Pre(r,t) = | —e7*!" P(r, t—t')dt’. 


This result is a consequence of the non-Markovian character (in the co- 
ordinates) of the Goldstein process in which persistence effects are introduced. 
In the difference equation from which (6) is obtained (for this equation see 
Goldstein (1951)), values of P at time & depend upon the previous values of P 
at time k—2, and if the dependence extends beyond k—1, the process is no 
longer Markovian. The stochastic process for the velocity, however, is Mar- 
kovian, since the L and the R in (2) depend only upon the values at the 
immediately preceding time. By a general theorem of Doob (1942), it is 
known that the autocorrelation function for a Markov process must be of the 
form e-*!"!|, This suggests that the hereditary function in (12) may have 
something to do with the velocity autocorrelation, and in fact, the exact 
relation will be demonstrated in the next section. 


Il, ENSEMBLE VELOCITY AUTOCORRELATION 
Consider the velocity autocorrelation 
(13) R(te, th) = (v(ty) v(te)) 


in which the braces indicate an average over the ensemble. This average may 


be expressed in terms of the second probability distribution, W(%, fi|v2, t2): 





(14) R(te, ty) = J fos, ti|v2, te)dv,dv2. 


This distribution is the probability (density) of a measurement of velocity 
v at time ¢,; and of a measurement v2 at the later time tf. The velocity is 
associated with the same diffusing particle. In the present model, that of 
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Goldstein, we consider velocities of v and —v only. In this case, for example, 
W(v, ti\v, t2) means the residual rightward-moving fraction at time f of the 
original rightward-moving fraction at time ¢,. For simplicity, the total number 
is assumed to be unity so that fractional number and probability may be used 
more or less interchangeably. Employing the notation N,(t) and N,(t) for 
rightward- and leftward-moving fractions at time ¢, respectively, the second 
probability distribution in terms of delta functions may be written as follows: 


(15) Wa, tile», t2) = 6(a1+2)8(v2+0) Ni (th) Ni (te) +6(21 —2)6 (22-2) Nati) Vr (lo) 
+6(11+2)6(v2—v) Ni (th) Nr(te) +6(21-0) 6 (02 +7) Ne (h) Ni (bh). 















Substituting (15) into (14) gives the result 
R(to, tt) = ULNR(h) Nr(te) + Ni (tr) Ni(t) — Neri) Vi (te) — Ni (th) Nr(t)I- 





(16) 









This expression can be explicitly evaluated from an examination of (4). The 
Laplace—Fourier transforms of the dependent variables RK and L are introduced 







by 
ae (FFG. 2») i [ae fal B@) a 
) nas oo hee 
with the inversion formula 
R(x, )) et [ . (4G »») = tertot 
(18) coe ws On J. dk Ss dp L¥(k, p) € : 





As may be seen from (18), the rightward- and leftward-moving fractions are 
readily expressed in terms of the transforms as shown below: 





+a woe a+ lex 
(19) Na(t) = | R(x, t)dx = — | REO, p)e'dp. 


27 J- 





Ny (t) may be similarly expressed. The transform version of (4) is 





(p—ikv—rs)Re(k, p)—r_L¥(k, p) = Re(k, 0) 


9 
(20) (ptikv—r,)LE(k, p) —r_RE(R, p) = Le(k, 0), 





where the right-hand terms are just the Fourier transforms ef the initial dis- 
tributions. Solving these equations algebraically and using the short cut given 
by (19), it is found that 






Na(t) = e'* (NR (0) cosh r_t+.V1(0) sinh r_t] 


y) ¥ 
(21) Ny (t) = e’*[NR(O) sinh r_t+N z(0) cosh r_t]. 





If at 4; = O there were, for example, rightward-moving particles only, then 
Nx (0) = land Ny, (0) = 0.* From (21), then, the velocity autocorrelation is 






simply 












(22) R@® =ve"™”, 


*In terms of the Lagrangian analysis of the motion of an individual particle we would say: 
at time ¢,; =0 its motion was certainly to the right (or, alternatively, certainly to the left). 











e 
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and neglecting the transient term, (11) may be written in the one-dimensional 


form 
+ 

(23) J (x,t) = -2f Ri(t')P (x, t-t')dt’. 
Ox 0 


Before considering generalizations of (23), it is instructive to inquire how 
the thermodynamics of irreversible processes leads to a similar result. 


III. DIFFUSION AS AN IRREVERSIBLE THERMODYNAMIC PROCESS 

For an ideal gas of like particles, a result equivalent to (23) can be derived 
under very reasonable assumptions. 

The first law of thermodynamics leads to an expression for the entropy 
change in a reversible process which may be written 
’ ee 
(24) és = TtP 
The first term here is a measure of the change of entropy, caused by a diffusion 
of molecules, for example, in momentum space. This term corresponds to a 
process in which there is diminishing localizability of the molecules in momen- 
tum space. The second term corresponds to the spatial spreading of the 
molecules, and may represent, when positive, a loss of information as to the 
whereabouts of the system in configuration space. Now, the total mean 
entropy for a system may be expressed in the form (to some extent, the view 


dV. 


of entropy proposed by Jaynes (1957) is adopted here) 


(25) S(t) = - | [ro r) In P(p, r)d’ pd’x, 


where we note that, in general, the entropy is time dependent. Utilizing the 
additive property of the entropy function, 


(26) S(t)=S,+S,=—- | Py(r) In Pi(r)d*x— | P:(p) In P2(p)d’p. 


In what is called an isothermal process, the second of these terms is constant. 
In the slow diffusion of a thermally isolated quantity of gas, there is a process 
of this kind. The slowness of the process ensures that the momentum dis- 
tribution will at all times be in equilibrium. Therefore, it may be said that all 
entropy changes are due to changes in the first term of (26) that represent 
configurational entropy. Taking the time derivative of (26) and observing the 
constant normalization of the position-dependent concentration P(r), it is 
seen that 

In Pi(r)d°x. 


97 we i Sen 
(27) ry 


dS = fe 
dt e 


Henceforth, the subscript will be dropped, since it is no longer needed. Employ- 
ing the equation of continuity (8) and vector identities gives 


(28) S -{ J-Vin Pact (In P)J-d. 


dt vol surf 
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As is usual in the thermodynamic theory of irreversible processes, we examine 
the rate of entropy production per unit volume. The contributions of the 
surface integral in (28) may therefore be neglected. Writing So for entropy 
production per unit volume gives 


(29) P Sy = J.X, 
where the abbreviation 
(30) X = -vVP 


is used, and X is called the ‘“‘force’’ or ‘‘affinity’’ associated with the ‘‘flux’’ J. 
Ordinarily one finds, or assumes, that the fluxes and the affinities are related 
by equations of the form 


(31) J = DagXg (summation on repeated indices). 


The Onsager theory demonstrates the symmetry of the Dag on a basis which 
appears inapplicable to the present problem. This demonstration includes 
treating the J, as being due to fluctuations in the number of particles per unit 
volume. This may be illustrated simply by integrating the equation of con- 
tinuity over a very small spherical volume to obtain the number of particles, 
N, contained therein. This is written 


(32) J-x = D:X-%, 

7 
where = represents the surface normal and the double dot represents the inner 
product of the diad D with the affinity vector X corresponding to (31). This 
rate of change, however, must be expressed more rigorously by 


(33) *(V(t+1)—N())e = D:X-x, 


where £ means ensemble average and rt does not approach zero but must be 
small compared with the time required to reach configurational equilibrium 
and yet large compared with the time required to reach momentum, or flux, 
equilibrium. If the latter has a relaxation time of the order 79, then + > 7» is 
required. For convenience, the more definite stipulation that t= 27» is made. 
It can now be seen that the fluxes have time to attain their prescribed 
values (31) in relation to the affinities. Expanding the expression in braces 
to second order in 7 and allowing 7 to assume its limiting value gives 


= ZaN d’N\, : 
(35) Xa tr al’ = —J-3>. 


Using the first equality in (32), the average values of the fluxes can be written 


as 


(36) 








ne 
he 


ay 
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This has the solution (a generalization of (31)) 


(37) J(t) = J - e *!D:X (t—t')dt’, 
070 





which is of the same form as the expression given in equation (23) for the 
Goldstein equation. The assumption that the system is always in a state of 
momentum equilibrium and that the entropy production rate, due to a 
changing flux distribution, is negligible is justified precisely by the retention 
of rate-of-change differentials of higher order, either as in equation (36) or as 
stipulated in equation (33). This point is discussed from a similar standpoint 
by De Groot (1952). 

If the system under consideration is an ideal gas, it is reasonable to consider 
the velocity changes to be governed by a Markov process. This is because the 
dynamics of a particle may be described by a Langevin equation in which the 
fluctuating (driving) force is caused by collisions with other particles. For an 
ideal gas there are no “distant encounters,” and so the collisions will produce 
intense, short impulses whose power spectrum is essentially white. The power 
spectrum of the resulting particle motions corresponds (by the Wiener- 
Khinchine theorem) to an exponential velocity autocorrelation function with 
a correlation time equal to the mean time between collisions. It is well known 
that the time between collisions is approximately equal to the time of relaxa- 
tion of the velocities to the Maxwell distribution; hence, the appropriate 
velocity autocorrelation function for the ideal gaseous system is exactly the 























hereditary function appearing in (37) above. 












IV. PROPOSED GENERALIZATION 






An examination of (23) and (37) leads to a formally straightforward generali- 
zation involving the velocity cross-correlation matrix defined by 


Rix(t') = (v,(t)o,(t+t’) ), 







(38) 









where, as before, the average is to be taken over an ensemble of similar systems. 
The generalization of Fick’s law now reads 






vt 
(39) J (r,t) = 2 | Rut’) P(r, t-t')dt’. 


k 
OX e 0 






Taking the divergence of this relation and employing the equation of continuity 
yields the postulated equation for turbulent diffusion: 






aP_ 
ot ax dx". 





ot 
| Rut’) P(r, t-t’)dt’. 


0 


(40) 






Since the velocity correlations need not be of the simple exponential form, this 
constitutes a significant generalization of the Goldstein equation. It is interest- 
ing to note that Batchelor and Townsend (1956) remark that a satisfactory 
theory of turbulent diffusion might be expected to take the form of an integral 
equation because of the persistence of motions. 
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The postulated equation leads to a well-known kinematical result, deducible 
independently of any diffusion equation, and due originally to Kampé de 
Fériet (1939). This relation expresses the mean square deviation of diffusing 
particles in terms of the velocity autocorrelation. This result is obtained by 
multiplying equation (40) by the square of the radius vector and integrating 
it over all space. After interchanging the order of the differentiation and inte- 
gration, it is seen that 


2 


d - — 3. {  - A Dp ae , 
(41) at (r*) = fa xT a oxt Rat) P(r, t-t')dt’. 


By a number of partial integrations coupled with the assumption that the 
concentration and its gradiant diminish at large distances according to 


(42) rP > 0 and r° = — 0 
or 


the desired relation is obtained: 


(43) = a *\ = 2 > |. ‘Relt’) dt’. 
It is required here, too, that the turbulence be statistically homogeneous, i.e., 
that the R;,(t') be independent of position. It should be pointed out that this 
result is strictly correct only if the source functions are sufficiently well be- 
haved. For example, an impulse point-source solution to the Goldstein equation 
has, in addition to the well-behaved wake solution, a wave-front solution with 
delta-function character. This latter contributes an additional transient term 
to the quantity given by (43). 

It is interesting to compare (40) with another diffusion equation, also 
involving velocity autocorrelation functions, which has been given by Pai 
(1957). He derives 


aP a°P 
(44) — Dix(t) tats 


where the diffusion tensor is given by 


(45) Di(t) = ae (v;(t" vg (t) +0, (t)0,(t')) dt’. 


Pai’s equation may be written with fo = 0. 


2 
(46) raf Radt = a" 


Pai’s analysis assumes a normal distribution of aie oniiiy per unit time. 
This is essentially an assumption regarding the form of the ‘propagation 
function”’ for the diffusing particle and completely determines the result. It is 
also an assumption which has the same validity as the central limit theorem 
and, like the latter, applies only to the results of a large number of independent 
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events. Specifically, Pai’s assumption is valid in the present context only for 
times that are long compared with a characteristic time of the diffusion process. 
In gaseous diffusion, this would mean the average time between collisions; 
here it is of the order of the relaxation time, or time scale of the turbulence, 
and is given by 

1 


= 
9) 
a 


(47) r= | Ri(t')dt’. 
0 
V. DIFFERENTIAL FORM OF PROPOSED EQUATION 
The proposed description (40) can be cast in purely differential form when- 
ever the autocorrelation functions have Laplace transforms expressible as a 
ratio of polynomials. The transform is defined by 


(48) ® .(s) = { Ru (tye “dt. 


( 
Writing formally, 
‘ *+ia a ttica 
( 19) Rix(t) (2 wl) j ei, (s)ds = (2 ri) Py aI | e ds 
. Ol} Jia 


—io 


aeaae =) 5(t) 
= ik at 0 ). 


Substituting this operational result into (40) gives 


(50) a = oe | oul : ) a | P(r, t—t’)dt’+6(r)5(2), 


ot ox' ax" Ji) ot ; 


where an impulse point source at the origin has been added to the original 
equation. Finally, from the properties of the delta function, 


4P y° ‘) 
(51) oF -—— P x = P+6(r)é(t). 
ot 


_ 7..A 
ot Ox' Ox 
When #,, can be expressed as a ratio of polynomials, we may write 
(52) Pi, — Vir Xs 


where x is the largest common denominator of the ®,. Operating on both 


sides of the preceding equation with x yields 


a 0 -— a (2) ‘ (2) 
(53) 2) at ~ axtax® Vita) PTO) x\ a7) 8). 


This is the differential form of the equation for turbulent diffusion. The point- 
source term may be omitted in applications but has been included here because 
of the unusual form it takes in the differential equation (though not in the 
integro-differential equation). 

The solution of the above equation, with the singular source terms, corre- 
sponds physically to the diffusion of unit amount of diffusate released instan- 
taneously at the origin at time zero. The Green’s function solution, however, 
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must be distinguished from this one. The Green’s function is the solution of 
the above equation with the source function replaced by one in which the 
operator x does not appear. The physical solution corresponding to a unit- 
impulse point source then can be obtained from the Green’s function solution 
by multiplying by the x(0/dt) operator. 

Some examples of the turbulent autocorrelation functions, their transforms, 
and the polynomial operators Y;, and x will be shown below. 










A. The Classical Case of Fickian Diffusion 
(54) Ra(t) = D bx 646) 
(55) @,(s) = Dbx 
Vi = Dbu 


ee ) x=1 















(56) 










where 6,(¢) is so defined that 


at 
(58) | 5,(t)dt = 1 t>0 
ei 





and 6, is the Kronecker delta. 












B. The Case Reducing to the Goldstein— Michelson Equation 


(59) Ralt) = (vt) be!" 










(60) ® .(s) = rove )0 ix /(1 +705) 
(61) Vn = 5 u(vk)T0 
0 
- a 4 
(62) x = 1+70 ar” 





C. The ‘ Phenomenological” Case (given by Frenkiel (1958, 1953) and applicable 
with suitable parameter values to a variety of turbulent regimes) 








(63) Rat) = dx (v7re7!*"”"* cos wnt 
(64) Pi.(s) = 715 ix (We) (Ste +L) /| (rex) + (Ste +1)"] 
ss 2 0 
(65) Vx = esletia( 147 #) 
2p 0 . 2 
(66) X= (+n = + (wT, ) . 
at 





D,. The Case Corresponding to the Turbulence Described by von Karman- 
Howarth (1938) 







(67) R(t) =rive 
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(68) Rut) = Raa(t) = Wide *!(1— 4ele)) 
with cross-correlations vanishing. 
(69) @11 = (vj1)/(s+8) 
(70) ee = ss = (01)(s+38)/(s +8) 


= wi) (248) 
2 = vas = (v?) C4) 


2 
(73) « (2+8) 


VI. DISCUSSION 


The diffusion equation given here was originally conceived, and has been 
presented, in a purely inductive fashion. However, since the first draft of this 
material was prepared, the author has learned that J. M. Richardson* has 
obtained, on very general deductive grounds, a quasi-differential diffusion 
equation, for the case of negligible spatial velocity correlations, that is equiva- 
lent to (40). The absence of any spatial correlation term in our equation is 
conspicuous; consequently, use of the new formulation is probably justified 
only in application to regimes in which the spatial coherence is negligible in 
comparison with the time coherence. In particular, application to plasmas with 
collective excitations may be precluded for this reason. From a consideration 
of the form of the equations of Richardson and from results obtained in a 
rigorous extension of the Goldstein model to two and three dimensions (to be 
presented elsewhere), it appears that the equation proposed here is valid 
providing at least one of the following two conditions are met: 

(1) the spatial velocity correlations are negligible; 

(2) space derivatives higher than the second are negligible. 

Applications of the diffusion equation proposed here to systems satisfying 
one of these criteria will be presented in subsequent reports. The effects of 
spatial velocity correlations will also be discussed. 


REFERENCES 

BATCHELOR, G. K. and Townsenp, A. A. 1956. In Surveys in mechanics, edited by G. K. 
Batchelor and R. M. Davies (Cambridge University Press, New York). 

Davies, R. W., DramMonp, R. J., and SmitnH, T. B. 1954. American Institute of Aerological 
Research, OSR-TN-54-62. 

De Groot, S. R. 1952. Thermodynamics of irreversible processes (Interscience Publishers, 
Inc., New York), pp. 222-223. 

Doon, J. L. 1942. Ann. Math. 43(2), 351. 

FRENKIEL, F. N. 1958. NACA Technical Memorandum. 1436. 

1953. In Advances in applied mechanics III, edited by R. von Mies and T. von Karman 
(Academic Press, New York), pp. 62-106. 
GOLDSTEIN, S. 1951. Quart. J. Mech. and Appl. Math. 4(2), 129. 


*Private communication. 





676 CANADIAN JOURNAL OF PHYSICS. VOL. 38, 1960 


Jaynes, E. T. 1957. Phys. Rev. 106(4), 620. 

KAMPE DE FEriET, J. 1957. Ann. soc. sci. Bruxelles, 59, Ser. I, 145. 

von KArMAN, T. and Howartu, L. 1938. Proc. Roy. Soc. A, 164, 192. 

MIcHELSON, I. 1954. American Institute of Aerological Research, AF 18(600)-476. 

Pat, SutH-[. 1957. Viscous flow theory, II, turbulent flow (D. Van Nostrand Co., Inc., New 
York), pp. 181-182. 

RICHARDSON, J. M. Hughes Research Laboratories. Private communication. 

TayLor, G. I. 1921. Proc. London Math. Soc. 20, 196. 





THE THIRD SPECTRUM OF COBALT (Co III)! 
A. G. SHENSTONE 


ABSTRACT 


A nearly complete analysis of the third spectrum of Co (Co III) is given, based 
chiefly on observations in the vacuum ultraviolet region. The chief structures 
found are 3d’, 3d° 4s, 3d° 4p, and 3d° 4d. The last of these is not commonly found 
in third spectra and it is here very incomplete. Some cases of unusual intensities 
in intersystem combinations are pointed out, especially ones in which the intensity 
in the combinations with the central of three / values is zero or very small. The 
complete line list is given. 


During the nineteen-thirties when there were many laboratories engaged in 
the analysis of atomic spectra, the third and fourth spectra of the metals of 
the first long period were neglected in favor of the lower and much higher 
ionizations. After the war we started at Princeton to fill this lacuna by extend- 
ing the analyses of third spectra which had been partially analyzed and by 
making analyses as complete as possible of the ones which had been passed 
over. In 1950 we started the observation of Ni III and in 1951 Co III. The 
preliminary term tables derived from those analyses appeared in Atomic 
Energy Levels, Vol. II, in 1952. A more complete description of Ni III was 
published in 1954 (Shenstone 1954), and additions to it were made in Pro- 
fessor Catalan’s laboratory in Madrid (Garcia-Riquelme 1955). Recently | 
have returned to the spectrum Co III and have been able to make such con- 
siderable additions to it that I believe it is now in as complete a state as is 
possible with the present observations. A description of the sources and 
instruments used in this investigation is given in some detail in the paper on 
Ni III (Shenstone 1954) to which reference should be made. 

The spectrum of Co III is somewhat more complex than that of Ni III and 
its description is not quite as complete in spite of the fact that the isoelectronic 
spectrum Fe II (A.E.L. 1952) has been very thoroughly analyzed, whereas 
Co II, which is isoelectronic with Ni III, has not. Mn I is also isoelectronic 
with Co III but many of its important electron structures are of a different 
character involving two excited electrons. 

In Fe II and Co III, the chief structures are 3d’, 3d° 4s, 3d°4p, 3d®* 5s, 
3d° 4d, all of which are of considerable complexity. The main multiplets are 
those arising from the transition 3d°(°D)4s—3d°(°D)4p. In Fe II they fall in the 
region 2700 A to 2300 A where observations of all kinds, including Zeeman 
effect, can be made with great accuracy, giving a consequent certainty in the 
analysis even though many other multiplets from the same electron transition 
are mixed with them. In Co III all of these multiplets are at shorter wave- 
lengths, 2000 A to 1800 A, where resolving power and accuracy of measurement 
are considerably poorer. There is a consequent overlapping of lines, with the 
result that several have to be assigned to more than one transition. On the 
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TABLE II 
Odd levels of Co III 
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TABLE II (Concluded) 
Odd levels of Co III 
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other hand, the structure 3d’, which in Fe II is mixed with 3d® 4s, is in Co III 
lower by over 40,000 cm~'. It therefore produces a set of strong lines reaching 
down to 650 A, completely separated from the main transitions. This is of 
considerable aid in the analysis, because the short wavelength lines must be 
transitions to one or another of the levels of 3d’. 

The details of the analysis are contained in the tables. Tables I and II give 
the even and odd levels arranged in terms in order of decreasing J. Table III 
gives a resume of all the identified terms and Table IV is the list of lines from 
which the analysis was made. The latter includes all of the observed lines 
except for a few of doubtful origin or very low intensity and it will be seen that 
most of the lines have been assigned. There still remain unassigned, however, 
a’number of moderately strong lines in the region 1300 A to 1600 A and a few 
even stronger lines between 1800 A and 2000 A. These groups require 
explanation. 

The shorter of the two regions contains mainly lines from the transition 
4p—5s and 4p-4d. The higher structures involved, especially the 4d, are very 
complex. Even 3d°(°D)4d alone yields 37 levels, of which 24 have been found. 
It seems, therefore, that the undiscovered levels of the complete structure can 
very easily account for all of the remaining unidentified lines of this group. 

The longer wavelength group is in the region occupied by the 4s—4/ transition 
and the lines are in all probability due to the unidentified higher levels of those 
two structures. The chief even levels not yet found are the high d’ 2D and those 





& 
| 
~ 
foe] 
° 
oO 
te 
° 
= 
ee) 
~ 
e 
oO 
Q 
A. 
un 
Q 
m 
je) 
e 


SHENSTONE 


HeP 2D YzD a? dz? 





An 


QQ 
~ 
- 


Ay? 
As? ds 


nm 2 


(b<udu 


AyD dy? 


suid} paasssqa 


Ill O-) JO sudo} paasasqo 
Il ATAVL 


(peu) su 


dy f'0 
dof‘? 


(de 9) 


(d1) 
(D1) 


(ds) 
(I) 
(De) 


(ds) 

(He) 

(ds) 
XU(Cs)9 PE 


PE 


99 FoF ZeSZe5 





CANADIAN JOURNAL OF PHYSICS. VOL. 38, 


TABLE IV 
Observed lines of Co III 
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TABLE IV (Continued) 
Observ ed lianas of Co III 
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844.866 118362 
845.793 118232. 
848.088 3 117912. 
849.210 117756. 
849.485 é 117718. 
850.067 ‘ 117637. 
850.424 : 117588. 
850.965 117513, 
858 . 663 
858.975 
865.567 


a*Py;-x *P2y 
a 4Py-x {Py 
a *G34-x 2 Foy 
a °G4y-w {Dsy 


a *G34-w °G33 


CcUWNK.o 
as 
coc 


a °G34-w 4 Foy 


CMHCNReENLS 


2G44 x 2s 
Geax "Has 


ay 


> i 00 00 Hm He 
Dim OMT On 


or 


iD www 
Cowen 


Omer 


a 4 Fy 3 {Day 
a *Foy-24F3y 
a 4F;, 3 4Fyy 
a 4*Foy-2 *Dyy 
a*P,4-2 *Day 


enc 
— OI. 


= 


co 


a*H,,—x *Asy 

a*Gyy-y 43 

865.898 a *Hyy-x *Ay 

870.007 4 a 2Hsy =¥ °G44 

874.294 4378. 2H44-x *Gsy 

876.594 é 2G,43-2 *sy , 
1030. 


877.881 

878.080 2Piy-y *Diy 1037. 96394.5 a ?F3y-y *Day 

878.543 *G43-v *Gay 1043 . 2 95854.9 a4Piy-24Py 

878.963 °G34-2 2H, 1044.2 2 95759.7 a4Poy-2 *Piy 

879.480 1046.76 é 95532 .9 a 4Pi4-24Piy 

880.950 °G34-x 4G34 1047. ‘ 95467 .9 a4Py-2'Py? 

882.025 4Pyy-x “Diy | 1048.87! ’ 95339 .§ 

883.154 2Doy-y 2F 34 1050.762 95169.0 a 4P24-2 *P24 

883.345 Gaya °F 5y 1050.§ 95149.6 a4Py-2'Py 

883.703 2Py-y "Diy 1053. : 94943.6 a4Pyy-s ‘Poy 

883.960 4Py-x 4D, 1064.56 93934. 

884.192 ‘Piy-x 4Day 1077 .5: : 92804. 

885.011 ‘Py-x Diy 1078.1: p 92749.: 

886.378 4Poy-x 4D3y 1079.06 ‘ 92670. 

887.777 2G34-2 * Foy 1088 . 486 91870.7 a*Pyy-2 Diy 

887.988 4Py-y "Diy 1089. 89° 91751.6 a *Poy-z *Doy 

891.902 a *Py4-2 2S, 1090. 55 91696.5 a4Py-2 ‘4D, 

893.045 a *Hs;-y *Hs4 1090 . 66 91686. 

893.095 a *Hyy-y *Hy 1092.: 91526.4 a4P14-2 4D; 

893.713 a*Dyy-y 2 Foy 1093 .066 5 91485.8 a4P,-24Dy 

897 . 686 a?Py-2 Sy 1095. 44: f 91287.3 a4P24-2 4Ds; 

899.025 1114.405 89734.0 a 2F3;-2 Gay 
a*Do4-y *Doy 1270. 646 : 78700. 


903 .730 
908.491 a *G44-2 Gay 1272.416 ‘ 78590. 


SON RE ORO ROIS 


Jie bobo bo 
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TABLE IV (Continued) 
Olneve ed lines of Co III 





A (vac.) I v | X (vac.) I v 

1286 . 626 Qo  Ga22.7 1473 .303 2h =67874.7 2 °Foy-f °Fey 
1308 . 850 L 76402 .9 || 1473 .528 5 67864.3 

1371.419 2 72917 .2 2 ®Day-f "Diy | 1474.278 3 67829 .8 

1371.779 3 2 °Diy-f Dy 1474 .953 3 67798 .8 

1372.536 0 || 1475.266 1 67784.4 2 °F5y-e "Gsy 
1376.576 0 7: 2 °Doy-f ®Day || 1475.635 20 67767 .4 2 °F5y-e "Gay 
1378. 208 2 a z*Di-f "Diy 1476 .366 5 67733 .9 

1378. 665 10 a 2 °D34-f *Dsy || 1476.895 5 67709 .6 

1380.775 5 7 2 °Di4-f *Dey | 1477.260 2 67692.9 2° Fyy-f ®D3; 
1383 .971 () 7225! 2 ®Doy-f ®Dsy | 1478.367 50 67642.2 2 °Fy3-e Gs, 
1384. 187 50 7224: 3 *Diy-f "Day || 1478.858 2 67619 .7 

1389 .079 10 71990.1 2 °D3,-f [Da || 1479 .393 3 67595 .3 

1396. 210 2 71622.5 2°Dy-e Py || 1479. 880 2 67573 .0 

1402.415 10 71305 .6 2 D4, € ®Dyy | 1481.905 100 67480.7 2 °F54-e Gey 
1403 .755 5 71237 .5 2*Day-e "Day || 1483 .363 10 67414.4 

1404.009 10 71224.6 2 ®D3;-e ®Day || 1483.800 20 67394.5 

1405 .032 5 71172.8 2°Dsj-e Diy || 1484.565 1 67359.8 b4Fyy-x 2 F 54 
1405. 451 ] 71151.5 2°D,y-e*D, || 1484.780 I 67350. 1 

1408. 808 5 70982 .0 2 ®D3y-e *Day || 1485.170 15 67332 .4 

1409 340 30 70955 . 2 || 1485.510 2 67316.9 

1412.012 5 70820 .9 a *D3;-y 4*D3y || 1487.962 2h ~=67 206.0 b4Fyy-w * Fy 
1412.974 10 (0772.7 3 ®Doy-e ®D3y || 1488.174 5 67196.4 

1413 .884 10 70727 .2 2 ®Dyy-e ®Doy | 1488. 465 1 67183 .3 

1425 .329 I 70159.2 3 °Dyy-e Fay || 1488.753 2 67170.3 

1425 .650 3 70143 .4 2 ®D3y-e ° Fay || 1488.975 1 67160.3 

1425 .984 | 70127 .0 || 1489. 243 0 67148.2 2 °Fiy-f *Day 
1427 .308 2 70062.0 2 Doe °Fyy || 1489.550 1 67134.4 

1428 .331 2 TOOLL.8 || 1489.899 20 67118.6 

1429 236 i) 69967 .5 2 *Diy-e Fy || 1490.511 | 67091. 1 

1429 .324 0 69963 .1 a 4D34-2 4H34? | 1491 .273 10 67056.8 

1430.177 30 69921.4 2 ®Dyy-e °F yy 1492 .528 5 67000. 4 

1430.530 30 69904 .2 2 °D34-e ®F3, | 1493.152 20 66972.4 b4Piy-x *Py 
1431 .323 30 69865.4 2 ®Doy-e °F, || 1493 .595 3H 66952.6 b4F;y-w 4*Fiay? 
1431 .800 10 69842 .2 s *Diy-e SF, || 1497.700 10 667690 

1431 065 10 69829 .2 2*Dy-e °F, 1498. 144 2 66749 .3 a 4De,-2 4Sy 
1432 .556 10 69805 .3 1498.830 20H 66718.7 

1434.259 200 69722.4 s *Dyy-e °F; | 1499.124 50 66705 .6 2 4D3y-f 4Dey 
1434.642 I 69703 .8 z°D;-e 6Fyy 1499 .877 5 66672.1 6 ‘Psy w *F, 24 
1435.426 100 69665 .7 3 ®D3-e eth 1500. 269 20 66654.7 2 4D24- f*'Diy 
1435 .840 5 69645 .6 2°Diy 5 Foy 1502.531 5 66554.4 4D\y-f 4D, 
1436.236 20 69626.4 2 °Do)-e o*Fy 1504. 362 3 66473.4 b4F., wsFy 
1445 .427 2 69183 .7 1505. 249 10 66434.2 b4P\4-x 4P24? 
1447 .216 3 69098 . 2 1508 .683 l 66283 .0 2 °F5y-e Day 
1448. 890 1 69018 .4 1510.489 2 66203 .7 b 4 F54-w 4Daoy 
1463 037 3 68350.9 and 

1463 .936 2 68309 .0 2 *Fy3- e *P 34? 
1464. 733 4 68271.8 1511.404 5 66163.6 

1464.932 2 68262.5 1511.730 3 661594 

1465. 101 2 68254 .7 1511 .920 5 6614101 2 °Fyy-e 8Dyy 
1465 .558 3 68233 .4 1512.997 5 660940 5 8 Fyy-e 8Dsy 
1466 .579 5 68185.9 and 

1466. 919 3 68170.1 2 ®Fy-e SGiy z*Dsy-f Ds, 
1467 394 10h =68148.0 2° Foy-f Diy 1513 .706 5 66063 .0 

1467 .663 10 68135.5 2 °Fy-e Gi, 1515.027 50 66005 . 4 

1468010 2h ~=68119.4 1515.638 2 65978 .8 

1468 . 992 10 68073 .9 2 °Fiy-e °Gay 1515.776 2 65972.8 

1470.282C? 1 68014.2 1515. 982 3 65963 .8 

1470. 700 I 67994 .8 1517 .070 2 65916.5 

1471 .140 2 67974.5 1517 .153 2 65912.9 

1471 .328 2 67965.8 2 °F34-f ®Doy 1517. 482 15 65898 . 6 

1471.575 1 67954 .4 a 4H3y-w 4F 34? 1517 .724 10 65888 .1 2 *Dsy-e 4 Fs, 
1472.019 20 67933 .9 2 °Foy-e °Gay 1518. 636 1 65848 .6 s4D2,-¢ 4 Fey 
1472.756 3 67899 .9 1519.326 15 65818.7 24Fyy-f *Dsy 





SHENSTONE: THIRD SPECTRUM OF COBALT 


TABLE IV (Continued) 
Observed lines of Co IIT 


d (vac.) D v r (vac.) 


1608 . 950 

54 Fyy-f *Doy 1618.491 
1620. 051 

SFyy é 4 Fy 1620. 254 
1622. 923 

§$Hyy-x Gay? 1629. 157 
1629. 666 

Fy f ‘Diy 1635 .521 
1645. 986 

4 Fy y—e 4 F yy? 1647. 100 
'Fiy-f 4D, 1649. 153 
1649. 265 

SF -e * Fi, 1651.783 
1652 .024 

et Fay 1652.182 
1652.791 

1659. 757 


boy | 


v 


62152.: 

61785. § 

61726.5 b?Hs54-x 7s, 
61718.7 6b 2F3y-x ? Fay 
61617.: 

61381.- 

61362.3 b?Fyy-x ? Fi; 
61142.6 a 4Gs5y-x G43 
60753 .9 a *Hey—-y *Hey 
60712 4H5,-y 4H; 
60637 .2 a ®Do;-2 4 F: 
60633. 

60540 . 6 

60531 

60526. 

60503 . 

60249. 


- 


1522 .049 3 65700. § 
1522 .885 20 65664. 
1523. 153 10d? 65653.¢ 
1524.113 15 65611. 
1524 .472 3 65596 
1524 .942 2h = 65576. 
1525. 163 2h 65566 
1525. 407 15 65556 
1526.734 30 65499 .: 
1527 .780 20 65454 
1527 .948 3 65447. 
1530.385 3 65343 
1530. 906 10h 65320. 
1531 .068 15k 65313. 
1532 .032 10 65272 
1534.268 200 65177. 
1536. 554 5 65080 
1536 .955 5 65063 
1539. 255 2 64966 
1539. 458 20 64957 . 
O94 100 64846 
.616 10 64824 
981 30 64767 
129 3 64761 
$11 30 64749 


02 00 Co C1 SD OS CO 
a 
Cie eho en 


co 


orcs 
nt 
ty : 
w 
ores 


oS 


se 


@ 
a 


Scun 


Qoyoanw ah 


b 27153-x "G43 
b*Ay4y x °G34 
a ®Do; 3 'Fyy 

6D4y-2 4*D5, 
a ®Dsy z ‘Do 
a *Hyy-x 'G} 
a ®D33 z 1F yy 
a *H6y-x *Gs; 
a Diy 24D, 

§Piy-x Go) 

6Do4-s 4Doj 

‘Hs x ‘Gay 

6D3 s 4D3; 


1661 .422 60189. - 
1663 . 500 ‘ 60114.: 
1665. 269 60050. 
1668 .032 ¢ 59950. ¢ 
1672. 133 59803 
1673 .325 59761 .: 
1673 .679 59748. 
’ 1677 .901 59598 .: 
S80 40 64730.0 24F\y-e 1F 1679. 481 5 59542 
309 5 64628 .: 1679578 20 59538. § 
5.413 20 64582: 1680. 734 1 59497. 
940 20 64518 1681 .074 59485. 
.036 5 64452 .3 1683. 180 1 59411 
3.015 5 64390.‘ 1687 .910 5 59244 
328 64212 1689. 283 3 59196 
746 64030.‘ : 1689.858 100 59176. 
1563 .158 63973 'Gyy-w 4 F 1690. 399 5 59157 
1564.090 63934 .9 b2H;4-w ? 1693 085 30 59063 
1567 .173 63809 |: 1695 .090 1 58993 . § 
1567. 917 63778 .9 b 2H; 1696.008 1000 58962. 6D33-s ®P 24 
1568 .308 637638 s*Pyy-¢ 1696.875 5 58931. ¢ 
1569. 264 63724.1 b? w? 1696. 980 5 58928 .2 a *Hs;-y Gay 
1572. 266 63602 .5 z®&Py4-e ® 1697.988 800 58893 .: 6D,\-2 *° Pry 
1573 .921 10h 63535 .6 1699. 880 58827. 2G4y-w 7s 
1576 .025 2 634508 1701.770 58762 .: 4Fyy-x 4 Foy 
1578 .0938 63367 .6 33-2 *Poy 1702 .506 : 58736.9 6 4Fsy-x 'Gs, 
1580.619 63266 1702.790 58727 *D,-2 * Py 
1580. 735 63261 s1Piy 1704.817 58657 .3. a 4Go;-y 2 Fay 
and 1705. 306 58640. é 
a 'G4y-w 7G3y 1707 .348 58570. ¢ 6Dyy-2 * Psy 
1581. 502 63231 .0 1707.951 6 58549. 6Do 4-2 * Poy 
1582.044 : 63209 .4 1709 .395 58500. : 
1584. 663 63104.9 a 4Gyy-w 4F sy 1711.531 58427 .2 6°Gs3y-w Hay 
1586. 180 63044 .6 1712.489 583945 
1587. 186 63004.6 b4Fay-y 2 Fay 1715.768 58282.9 a *Gyy-y 7H, 
1588 . 642 62946.8 b 2H,y-w Ga, and 
1593. 446 62757. 1 b°Gyy-x 2 Fay 
1504. 182 62728.1 a4Goy-w t Fy 1716.251 58266 ®Di4-s ° Poy 
1594.953 62697 .8 a‘4Gyy-w *Dsy 1718. 887 58177 1G5y-y 2Hs4 
1601. 903 ‘ 624258 1719.383 ‘ 58160 
1603 . 563 c 62361 .1 @ 4G3;-w *D2, 1720.068 58137. 
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bobo 
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CHKnmowa 


pt feet fos 
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oS 


bobs 


4H33-3' 


1608 .777 t 62159.0 1721.678 : 58082 
1608 . 864 ‘ 62155.7 a Gey-w'Dy 1722 .911 58041 


5 

2 

4 

2 F2; 
1607 .410 62211.9 1721.151 58100.7 a 4Hyy-s 7 Fs, 

9 

3 


'F sy x ‘Fy, 
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TABLE IV (Continued) 


Observed lines of Co III 


v d (vac.) 


58027 .: 1774.577 56351.5 a 2154-y 7164 
58020 .3 a 4Gsi-y "Ay 1776.630 56286.3 b 4Dsy-w *F34 
58005 .6 a ®§D54-2 ®P 34 1777 .145 56270.0 b 4Fyy-y 4Gsy 
57932 .9 a *He6y—v 4Gs4 1778.091 56240.1 a *D.4-2 °F, 
57913.1 a 4H34-2 ?F sy || 1779.536 56194.4 b4Pi4-y 4D, 
57867 .3 b °G34—-x 2 Foy 1779.896 56183.1 b4Pi4-x *D2, 
57812.7 b4Fyy—-x * Fay 1780 .046 56178.3 a ®D24-s ° Fo, 
57781.1 b?*Hyy-y 7H 1780. 686 56158. 1 
57724.8 b *G34—-x °F sy 1782.195 56110.6 
57718.6 a *Hs4-y *Gay || 1782.966 ;¢ 56086 .3 a °D24-2 °F 
57682.3 a +Hs4-y 'Gsy 1784.055 56052.1 b 4F3y-y *G4, 
57623 .6 a 4Hyy-y *Gsy || 1784.790 56029. 
57593 .3 a *D24-2 ® P34 1784.944 : 56022 .6 
a4Gy,-y 7s 1785. 116 : 56018 .8 b 4D2y-w *Ds; 
b 4 Fyy-x 4G || 1785.705 : 56000 .3 6 ?G34-y 7H 
b 4 Fyy—y *Gay | 1785.965 55992.1 a °D,y-2 °F, 
a *H33-y 'Goy || 1786.342 55980.3 b 4F.y-y 'Gsy 
a *H34-y *Gs4 |, 1786.679 gf 55969.8 b4Di4-w 4Doy 
C °*G4y-w ?F 34 | 1787 .082 55957 .1 a ®Dyy-2 ° Fy 
and 
b 2 Fyy-y 2 Fy b4*Dy-w 4D, 
a *H34-y *G44 | 1787.502 55944 b 4Dey-w *Doy 
6 *H;,-y *Hs4 | 1787.575 ; 55941 b4Dy-w'D, 
| 1788.627 f 55908 . 
b 4 F4-y "Gay || 1789.070 55895 .0 a *Dyjy-s ° Fay 
1} 1789.373 55885 .5 6 4Fiy-y 4Goy 
b *P24-y *Day || 1789.549 55880 .0 b 4D34-w 4D. 
1790. 258 55857 .9 a 4Gsy-y *Hey 
b 4P.4-x 4Dsy 1790 . 389 55853 .5 b 4Diy-w 4Diy 
56987 .8 b 4D24-x * Foy | 1790.892 55838.1 b4Dy-w'Diy 
56948 .3 a ®Dyy-2 ° Fay i} 1791.153 55830.0 a 2J54-w 2H, 
56920.0 b?2Pi4-y *7Dey | 1791.277 55826.1 a °*D,-2 °F, 
56898 .0 b 4Fyy-2 2 F5y | 1792.144 55799 .1 a 4Hs4-2 7165 
56868 .1 b 4F34-2 2 Fo, | 1792.410 55790.8 a ®*Dy-2 *Fiy 
é 56845.8 b4Doy—x ?Fyy || 1793 .924 55743 .7 b4Piy-y 4Dey 
5000 56806 .8 a °D,4~-2 °F sy || 1794. 804 55716.4 a 21ey-w ?Hsy 
1 56783 | 1795 .426 55697 .1 b?Piy-x 4 Foy 
10 56774.1 b*G4y-x #Hsy || 1796. 200 55673 .1 a 2154-w *Hs, 
100 56706.5 b 4D34-x ?F yy | 1796.664 55658 
15 56704.4 b?Fyy-y ?Foy 1797 .095 55645 .4 b?Py-y*Diy 
0 56643 .9 b 4F34-2 27F 54 || 1798.064 55615.4 b 4*Pay-y *Dsy 
1 56639. || 1798.876 55590.5 b4Py-x *Diy 
56590 .4 a °D34-2 ® Foy || 1800.469 55541.1 b4Diy-x 4 Py 
50 56583 .2 b 4D24-w *G3; || 1800.974 25 55525.5 b4Dy-x *Py 
5H 56561 || 1801.030 25 55523 .8 6b 2Pi4y-2 72S, 
200 56553 .5 6b 4D3,-w 4 Fyy | 1801.940 0 55495. 
100 56546 .0 b 4 Fyy-2 * Fay | 1802.917 1 55465.7 a 4Gyy-y 4H, 
10 56513.: || 1805.535 500 55385.2 a 4Gy,-y 454 
500 56498 .5 a "Ds4-2 °Fsy 1806 .096 10 55368 .0 a 4Gsy-x 'Gay 
10 56481 |: 1806 . 350 10 55360. 
100 56457 .0 6b 4 F.4-2 Gay | 1806 .492 5 55355 «§ 
200 564381 b4Dy,y-w * Fay || 1808.384 300 0 DG gy-x Gay 
300 56426.0 b4D,y-w Fy, 1809 . 233 40 
and 1809 . 904 5 5 
b 4D2y-w Fs, 1810.464 15 234.5 a tHyay-s Gay 
56419.5 b ADs, -w Gay 1811.317 100 55: 4 b4Py-y'Diy 
56415.‘ 1811.466 400 55203 .9 a 4Gsy-y 4Hay 
56412.0 b4Dy-w Fy 1811. 694 3 5514 b4Dy-x *Piy 
and 1812.336 5 f 4A a tsy-s*Gay 
b4Dzy-w 4 Fy 1812. 550 50 bf 9 b4Dy-x'*Pyy 
56394 276y—-y Tey 1813 .044 50 a *Gyy-X 1G5y 
56383.5 a °Dsy 2° Fay 1813.186 300 a ‘Gay-y 4Hs4 
56356. ! 2754-y *Ley 1814. 084 10 a *Hyy-y *F yy 


57515 
57476 
57436 
57393. 
57365. 
57362 
57322 
57284 
57269. 
57261 
57245 
57151 
57109 
57082 


PON INK POKMWOMNOW 


& 00 
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ed 
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Qn 
ro 














1814. 
1814. 
1814. 
1814. 
1815. 
1815. 
1815. 
1816. 
1816. 
1816. 
1817. 
1817. 
1817. 
1818. 
1819. 
1819. 
1819. 
1819. 
1819. 
1820. 
1821 
1821 
1821 
1821. 
1822 








1822. 
1823. 







1823 
1823. 
1824. 
1825. 
1825. 
1825. 
1826. 
1827. 
1829 
1829 
1830 
1830. 
1830 
1830 
1831 
1831 
1832. 
1832. 
1833 
1834 
1835 














1835 
1835 
1836 







1837. 
1837 
1838 
1839 
1839 






1835. : 


d (vac.) 


219 
573 
683 
865 
063 
596 
686 
084 
250 
617 
518 
626 
750 
684 
A) 
070 
261 
330 
733 
064 


.004 


. 262 


688 
766 


046 


215 
079 


414 
622 
874 
365 
464 
947 
716 
094 
392 
674 
093 
581 
780 
870 


.439 
.916 


201 
784 
876 
840 
000 
255 
617 
687 
200 


630 
840 
659 
535 
636 


I 


100 
2 
100 
100 
20 
200 
200 
5 
10 
20 
100 
100 
20 
200 
2 

1 
30 
200 
100 
100 
1 
400 
400 
400 
100 


1? 
1000 


100 
150 
100 
400 
300 
750 
1 
300 
5 
200 
2000 
300 
20? 
150 
2000 
750 
400 
1 
50 
75 
5000 
100 
150 
20 
300 


1000 
100 
9 


50 
20 





v 
55120. 
55109. 
55106 
55100. 
55094. 
55078 .3 
55075 .6 
55063 .5 
55058 .5 
55047 .4 
55020. 1 
55016.8 
55013. 1 
54984 .8 
54975 .0 
54973. 1 
54967 . 4 
54965 .3 
54953 . 1 
54943 .1 
54914.8 
54907 .0 
54894. 1 
54891 .8 
54883 . 4 


54878 .: 
54852: 


54842. : 


54835 
54798 
54783 
54780 
54766 
54743 
54731 
54663 
54654 
54642 
54627 
54621 
54618 
54601 


54587 . 


54579 
54561 


54529. 
54500. 
54495. § 
54488 .: 
54477. 
54475: 
54460.¢ 


54417 
54411 
54387 
54361 
54358 


6 


aac he 


TABLE IV (Continued) 


SHENSTONE: THIRD SPECTRUM OF COBALT 


Observed lines of Co III 


a 4H33-2 G34 
C°G3y-w *Ayy 
a *G24-x 4 Foy 


b 4F3y-x *Day 
b4*Dey x Poy 
b 4 Fyy-2 754 

a *Hyy-2 Gay 
b *Hsy ¥ 1H; 


b SFiy x ‘Diy 
a 4‘Hyy - {Fy 
C*Gyy-w 7s, 
a §G34-x {G34 
b *G34-x *G3j 
b*Piy-y *P4 


and 


b*Hy-y ‘Hs, 


a 4H 64-2 *1s4 
a *Hey-2 ‘164 
and 
a ‘Gs “x 4Gsy 
a 215y—w *Gay? 
b 4D3y-x *P2y 
b 4 Fyy-x *Day 
b 2 F3y-y *Doy 
a *H3y-2 Gay 
a *Gyy-x Gay 


x ‘Ds, 
x ‘Dey 
y*Dy 
2 ‘Poy 
z {Tsy 
y ‘Fy 


b4Fy 
b4Fy 
b 2 Foy 
a ‘Dy 
a ‘Hs, 
a ‘Hs, 


a 4Hsy Zz ‘Tey 
a {Doy-s *P iy 
a 'Gs;y-x 1 Fyy 
b ‘Fy y ‘Dy 
b*Piy-y 4 Fay 
a *G34-x Gay 
a ‘Hey Zz ‘Thy 
{‘Dyy-3 *Py 
4*Fyy-x *Dsy 
2Hsy-x “Gay 
a *Gyy-x * Fay 

and 

( G34 x Foy 
a *Hyy-2 4154 
a “Gay-X Gay 


= 


co 


€ Gay x °F yy 
b *Fay-y *Day 


1851 
1852 
1852 


1853. 
1853. 


1854 
1854 


1854. 


1855 
1855 


1857. 
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1859 
1860 
1860 


1861. 
1862. 
1863. 
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1863 
1863 


1864. 


1864 


1865. 


1865 
L866 


1866. 


1867 


1867. 


L868 
1870 
1870 


5050 
5157 
3.514 


7.300 
.825 
.030 
9.299 
9.464 
9.932 
.503 
.780 


.509 


937 
645 
919 
266 
841 
194 
393 


763 
331 

954 
657 
487 
510 
428 
866 
775 
660 
134 
467 
615 
$26 
187 
519 
424 
456 
497 
615 
490 
930 
796 
O12 
634 


125 
125 

25 
150 
300 
200 
300 

10 
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500 
10 
10 
125 
400 


400 
25 
50 
50 
50 
50 
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400 
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100 
LOO 
20 
100 
50 
30 
30 
30 
io 
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54301. 
54291. 
54285. 
54249. 
54246. 
54243. 
54231. 
54209. 
54208. 
54198. 
54181. 
54169. 
54166. 
54156. 


54133. 


54117 


54111. 
54074. 


54069 
54056 


54039. 
54031 .< 


54010. 


53997 


53976. § 


53968 
53958 
53942 
53931 
53926 


53915 
53898 
53880 


53831 .< 


53807 
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53738 
53712 
53686 
53673 
53663 
53659 
53653 
53642 
53633 
53607 
53606 
53576 
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53535 
53510 
53475 
53457 
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b {Foy oe {Dy 
b *F2y-y *Day 
a ‘Gyy-x {G54 
a‘*D,-2 ‘Py 

b4Py-2*Diy 


€ *Gay-w ‘Fy 
a *G3y-x ‘Fy 


a {Diy z ‘Py 

a 4Hyy-2 ‘Tay 
and 

a *G3;-2*Ayy 


b*Py 2*Piy 


C *G4y-w 7G 4y 

b 4Psy-y ‘Pay 

C *G34-w 7G34 

b*Piy-y *Piy 

b 4 Fyy-2 *Gay 
and 
€ 2G34-w °G4y 

b*Py-y *Py 
and 
b*P14-2 *Py? 

a *D.4-2 *P2y 

a *Gyy-x t Fay 

a *Hy4-2 *T44 

b *Hs4-x *Gsy 

C *Gyy-w *F sy 

a *Gay-y 74 

a‘*D,-2'Piy 
and 
b*Gyy-y ?F sy 

b*Hyy-2 *7Ay 

c °G34 w ‘Fy, 

b*Py-2 4S, 


a ‘Gay x ‘Gs, 
bt Fyy-2 Gy 


b Fy 2 ‘Fy 
b*Hsy-x * Fay 
a *Hey-2 *Hs, 
a *Igy-x *H sy 
b*Fyy-y {Hy 
a *H4-2 ‘Hey 
b *Hs4-2 *Hsy 


a ‘Hy 2 *Hs,; 
b 4 Fyy-y *F sy 
a*D,y 2*P2y 


a *Hsy-2 ‘Hy 
€*Gay-w *Dsy 
b 4 F34-2 "Gay 

b *Hyy-x *Gsy 
b ‘F;, y ‘Fy 





CANADIAN JOURNAL OF PHYSICS. VOL. 38, 1960 


TABLE 1V (Continued) 


d (vac.) v A (vac. ) v 


1871.870 53422.5 a *Hsy-2 4Hs3 1912.465 52288. 
1871 .952 53420.2 a 4H33-2 4H, 1912 .952 ‘ 52275. 
1872 .532 53403 .6 a 2I54—x *Agy 1916.112 ‘ 52189. 
1872.575 53402.4 a 4H54~-2 *He6y 1917 .321 : 52156. 
1873 .014 53389 .9 b?Hsi~y *Gay 1919.120 5 52107. 
1874.355 53351.7 b4Fyy-y 4 Fy 1919. 980 52083 . § 
1874 .822 53338 .4 a 4H 44-2 4H 44 1922 .000 52029. 
1875 .094 53330.7 6b 4F3y-y 4 Foy 1923 .801 ; 51980. 
1875.609 53316 1924 .529 ¢ 51960. 
1877 .464 53263.3 b 4Fiy-y 4Dsy 1924.72: 51955. 
1877 .544 53261.1 b4Fyy-y 4Fiy 1925.26 51941. 
1879.170 53215.9 b4Fyy-2 4T 54 1925.51: 5 51934 
1879. 244 53212.9 a 4Hyi-2 Hs 1925. 56: t 51932. § 
1879 .385 53208 .9 b ?G3y-y 2 Foy 1925 .§ 51923. 
1879 .753 53198.5 b ?Fo-y 4H34 1926 .365 51911. 
1880 .321 50 53182. 1926 . 656 ‘ 519038 
1880.449 30 53178 2Hy4-y Gay 1927.16 f 51888 
1880.912 50 53165.7 b4Foy-y Fy 1927.7: s 51874. 
1881.080 150 53161 \G44—2 °F 34? 1928. 46 51854 
1881.427 150 53151.1 a4H3;-24Hy 1928. 51851 
1881.702 1000 53143.4 a 4Hei-2 4Gsy || 1928.80: 51845 
and 1929 .756 51820 
a 4H yy Z 'G34 1930. : 51800 
and 1932.44: 51748 
b4Foy-y4Foy || 1932.766 51739.: 
ISS1 300 53138 1F154-2 *Gay | 1933.25 51726. - 
1882 150 53125 ‘Fiy-y 4 Fiy 1933.74: 51713 
1883 200 53098 1H34-2 'Goy 1934. 27- 51699 
1884 2 53063 .6 6 2 Foy-x 4 Foy 1934. 73+- § 51686 
1885 15h = 53037 IGs4-y 4Gsy 1935.02: 51679 
1886. 50 53009 'Fs,—y 4D34? 1936.39: 51642.: 
and 1936. 78: 51632 
b 4 Fiy-y 4Foy 1936.93: 51628 
1886.742 200 53001 1937.66 51608. 
1888.34: 10. 52956 1H 5-2 'Gsy 1940 51542 
and 1941.46 : 51507 
b2Fyy-x 'Goy 1941.7: 51500! 
[S88 .451 10 52954. 4 Fyy—c 4144? 1942.36! 51483. 
1889 .090 10 52935! ; aT 
1889 .321 20 52929 
1890 .642 2 52892 1942. 46 51480. 
1892.011 150  52853.8 ¢°Gyy-x 2H 1942.796 51472.: 
1893.095 20 52823. 6 b4Foy-y 4Ds; 1945. 23: 51407 
1893911 5 52800.8 a 'G34-2 2 F sy 1946.79: 51366 
1895.368 500 52760. 27, 2hy 1947 626 51344 
1895 .495 30 52756 1948655 51317. 
1896 .031 ] 52741.8 33-2'G 1949 ) : 51308 
1897 .560 5 52699 .3 6 4F34~s 41; 1949 53: ‘ 51204 
1899183 5 52654.2 b?F,4—x 'F; 1949. 805 51287 
1899 .795 50 52637 .3 ¢ ?G4y-x 7, 1950. 668 
1900. 763 50 52610.5 b 2Fy;-2 2H, 1950. § 
1901.357 300 52594 _ 1950.96 
1903. 199 5 52543 1952. 158 
1905.35 300 5248: °G34-y *# ay 1953 .§ 
1905.674 20 2474.9 b2*Py-22Py 1954 
1905.810 5 5s Gey IGay 1954 
1906. 265 30 52458 6 b ? Foy-x 4Gay 1955. 
1907 .458 30 52425 .8 1955.7 
1909. 666 52365.2 b?P y-2?P, 1956 
1910. 151 1 52351.9 b2Fy)-x 1Gyy 1956 
1910. 840 300 52333 *Gai—y 2s 1958 .£ 1 
1911.174 20 §2323.9 a4D3,-2 4 Foy 1959. 41. 51035.7 
7 
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b 2Fay-x "Fay 
a Dy, 2 ®D3y 
b*Hsj-y Gay 
a*Do; 2{*Fy 


CweOnwc 


> OO be 


b °Day—-w 2F 44? 


b 4Fyy-2 4G 


we 


2H yyy "Gay 


SCwOWNwo 


~ a a —_ 
oo onNs 


oon 


Pm StS OS No 


host 


and 1961. 450 


a'Gsy-y Gs} 1963 243 3 50936 





SHENSTONE: THIRD SPECTRUM OF COBALT 


TABLE IV (Continued) 
Observed lines of Co III 





d (vac. ) A (vac. ) I v 
1963 .743 50923. : 41Do4-2 2052 .770 10 48714.7 6 
1964.765 < 50896. 1G4y-x 4D3 2053 .769 200 48691 
1966 .029 : 50863 . ¢ 2053 . 977 2 48686 
1966. 680 : 50847 24-2 7154 2056.869 100 48617 
1969.312 2 50779.2 b ?Fay-y 4Gs4 2062. 825 10 48477 
1970.054  ¢ 50760 4D34-2 4D34 2072.67 48246 |! 
1971 .889 50712. 4D\3-2 ‘Dy 2073.91 48218 .: 
1973 . 667 50667 . 2P 13-2! 2074.18 48211 
1973 .767 50664.5 6 ?F34-y 4G 2075.78 48174 
1974.003 50658. £ 2078.01 48125 
1974.883 2 50635.9 b 2Hyy-s 2754 2080.40 48067 
1977 .031 Z 50580. ¢ ‘Dos *De 2081.67 48038 
1978 .948 5 50531. ¢ Gsy 2°G. 2081.90 48032.‘ 
1980.113 2 50502.‘ 1D4-2 2082.91 48009 | § 
1981 .000 50479.6 a'Gsy-y ' ey 2084.89 A794 
1981 .345 50470.8 a4*Dy-s 4D, 2085.61 ; AT9AT 
1982 .049 50452 

1982 .361 50444. § 
1984 .324 3 503905 2Hyy—x 4D; 2088 .06 47891 .° 
1986. 704 3 50334 .6 2089. 24 t 17864.: 
1987 .197 30 50322 2 Foy 3Goy 2090 09 47844 
1987 .952 50303 2091.17 47820.° 
1988 . 806 50281 2 Fyy-y 4G 2092.88 47781. 
1989.598 — - 50261 2H4-2 716 2097 .62 3 47673. 
1989 .645 50260. ’ 2098 .30 47657 
1990. 841 50230 2105.84 a 47486. § 
1991 .929 50202 2111.82 17352 
1992. 158 50196 2120.49 é 47158. § 
1993 .625 50159 2121 8) 47127.§ 
1995 .397 5 50115.: 2124.74 17064 
1999 .913 3 50002 2134.82 46842: 
2000. 447 , 49988 |! 2135.55 46826: 
2000 . 702 49982 2136.81 46798 
2001 .735 49956 2137 .98 46773 .° 
2002. 559 19936 2172.94 16020. : 
2003 . 437 49914. * 94 45580. 
2004 .555 , 49886 
2004 939 ’ 49876. 
2006. 194 r 49845 
2011. 252 19720 .: 
2012.266 : 19695. : 
2013 .379 : 19667 
2014.53 ; 19639: 
2016 .473 ; 195901! 
2016.685 49586.: 
2020. 064 F 19503 
2020 .775 i 19486 
2025.118 < 19379 
2029 . 229 19279 
2031 .33 : 19228 
2032 464 $9201. 
2033 .728 : 19170 
2034 64 { 19157 
2034 939 49141 
2039. 820 5 19023 
2040. 498 5 49007 
2040. 884 48998 | « 2978.01 
20438. 700 2 18030. ¢ 2986.13 
2044. 882 § 48902. 6 154-2 2990.04 
2046. 817 48856.: — 2 2991.89 
2048 .019 : $8827 3010.92 
2048 . 409 : {S818.- wt 3018.41 : 
2049 .093 48802 AF, 3024.92 r 33049 
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TABLE IV (Concluded) 
Observed lines of Co III 


d (air) v A (air) I v 


3075.89 5 32501 .5 
3 
2 
) 
) 


3305.38 15 30245.1 a 4Gsy-2 4 Fy 
3097.74 2 32972. 3451.25 10 28966. 
3116.68 32076. 1 3526.24 2 28350. 

3151.40 2  31722.8 | 3634.21 27508.5 a *Diy-2*Py 
3180.64 2  31431.2 3636.31 27492.6 a 4Dy-2*P, 
3232.11 2 30930.6 a ‘Gs *Fy 3667.52 27258. 

( 27186.7 b4Dy-2*Piy 
27160.8 b4Do4y-z*Piy 
26570.6 b ‘Dey-2 *Poy 
26431.7 b4D3y-2 *Pay 


3249.24 2 30767. 3677 . 23 
3259.68 20 30669.1 @ 4G34-2 4 Foy 3680.74 
3269.23 2 30579.5 a *Goy-s 4 Foy 3762.50 
3287 .68 10 30407.8 a4Gyy-2 4 Fs, 3782.27 


Cre White et 


based on the 3d° terms *D(part), 'S, 'D(part), 'F, b°P, b°F, b'G. It might 
seem that some of these levels could be found from the remaining lines, but a 
thorough search has been unsuccessful. 


INTERSYSTEM LINES 


The terms of different multiplicity are connected by the numerous combina- 
tions which are expected in a spectrum which shows evidence of j—j coupling. 
Some of the intensity of such lines is due to the sharing of properties by 
neighboring levels of the same J value but different origin as in the case of 
wiDs;, w'Fy3, wGs3, and x*F33; but there are many for which such an 
explanation is not tenable. The intensities of the lines are also frequently very 
anomalous as in the following case (Table V). 

TABLE V 


3h 2h 

2G 44 2 30 
34 t ‘ 50 

2 Fy, f 0 100 

. : 00 100 


On the evidence of these intensities alone 2?Gy might be ‘/4 or ‘Gy or 
‘Fy and 2°G3, could be 4H or *F33. Likewise, 2?F could be given a number of 
different names. Especially difficult to understand is the not uncommon habit 
of intersystem combinations of missing the central of three levels of successive J. 

Another curious example of intersystem combinations is given in (Table VI). 


TABLE VI 


Hs, 4h Fs, 


YH 200 
64 20 200 
44 
34 
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HIGH EVEN TERMS 


In Fe II the high terms due to 5s and 4d are separated by over 4000 wave- 
numbers and they are therefore easily distinguished. In Co III these electrons 
have nearly the same energy so that terms due to them overlap completely. 
Each of the structures produces a §D term and both have now been found in 
Co III. These terms in Fe II make combinations of very different intensity, 
the 4s term giving strong lines, while the 4d *D is represented by a few levels 
giving very weak combinations. There are two °D terms in Co III having these 
characteristics and they were assigned to 5s and 4d by comparison with Fe II. 
It is regretted that the only high level reported to Mrs. Sitterly for inclusion 
in A.E.L. was wrongly identified as 3d°(°D)5s *D4, whereas it is now with some 
certainty identified as the 5d *Dy level. The term e °G is based on single lines 
except for one satellite, but there seems little doubt of its reality. A comparison 
of the 4d and 5s terms of Fe II and Co III is given in Fig. 1. 





. 44 SF 
83} OOO Fell Coll S 


se ee 
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Jo$ tp ab3y 4s sy 


Fic. 1. 4d Levels of Fe Il and Co IIL. 


TWO-ELECTRON TERMS 
In Ni IIL the observations of the short wavelength end of the spectrum 
made it possible to find several of the terms of 3d° 4s(®D)4p and subsequently 
Garcia-Riquelme identified two even levels as belonging to 3d* 4s?°D 
(Garcia-Riquelme 1958). Although none of the levels of 3d°4s 4p of Co III 


have been found, it is possible by a comparison such as is given in Catalan, 
Rohrlich, and Shenstone (1954) to place 3d° 4s? °S at about 182,000. Attempts 


to find such a term have been unsuccesstul. 
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Observations of the pulsed cobalt hollow cathode have been made at Lund 
by Garcia-Riquelme and Velasco and they have kindly put the line list at my 
disposal. The wavelength range is 2308 A to 3998 A and the list includes most 
of the lines of Table III where the lists overlap. In the longer wavelength 
part, one might expect the lines due to the transition 4d—4f, but I have been 
unsuccessful in a search for such combinations. The long wavelength observa- 
tions probably must be made more complete before such a complicated 


structure as 3d° 4f can be identified. 

In A. E. L. the value of the series limit 3d°°D is given as 270,200 cm™. 
The additions to the analysis given in this paper include no evidence on which 
to change that value and it is therefore to be accepted as the best available 
at present. The corresponding ionization potential is 33.49 volts. 


REFERENCES 
Atomic ENERGY LEVELS, Vol. II, page 55. 
CatTaLan, M. A., Ronevicn, F., and SHENSTONE, A. G. 1954. Proc. Roy. Soc. A, 221, 421. 
GarcIA-RIQUELME, O. 1955. Ann. reale soc. Esp. (A), 51, 59. 
—— 1958. J. Opt. Soc. Am. 48, 183. 
SHENSTONE, A. G. 1954. J. Opt. Soc. of Am. 44, 749. 





NOTES 


ASYMMETRIES IN THE BETA DECAY OF POLARIZED NEUTRONS* 


M. A. CLARK AND J. M. ROBSON 


It has been shown by Jackson et al. (1957) how various asymmetries in the 
decay of polarized neutrons give important information about the nature of 
the interaction responsible for the decay. Allowing for failure both of parity 
conservation and of time reversal invariance, these workers showed that the 
angular distribution of the decay products includes the following terms: 


45 -B432 85 pi. (P: x Ps) 
eB I 6; FE NES OE: 

using their notation. The terms A and B are non-zero only if parity is not 
conserved and D is non-zero only if time reversal invariance is violated. The 
coefficients depend in various ways on the coupling constants and thus measure- 
ments of .1 and B provide information about the nature of the interaction. 

We report here a final value for our measurement of D, and also the results 
of a measurement of B. 

We have continued our measurement of the coefficient D as already described 
elsewhere (Clark et a/. 1958a and 6) and we have obtained a final value of 
D = —0.14+0.20. For full violation of time reversal invariance one expects 
D = 0.45, using a mixture of vector and axial vector interactions, where 
Cy Cy) = 1.19. Our result is thus consistent with time reversal invariance 
in the neutron decay. (Burgy eft al. (1958b) have reported D = —.04+.07.) 

The principle of the measurement of the coefficient B in the neutron decay 
can be seen from Fig. 1 in which a neutron with spin ‘‘up”’ is assumed to decay 
with the proton going to the left, as indicated by the arrow ‘“‘p’’. If the electron 
is detected in the direction £1, as shown in the left-hand diagram, then, by 
conservation of linear momentum, the antineutrino must go upwards. If, on 
the other hand, the electron is detected in the direction #2, as shown in the 
right-hand diagram, the antineutrino must go downwards and be approxi- 
mately antiparallel to the neutron spin. Thus a comparison of the (p, £1) 
and (~, £2) coincidence rates will, after suitable corrections for geometry and 
electron energy, yield a value for B, the coefficient of the correlation between 
the neutron spin and the neutrino direction. We have used the apparatus 
already described (Clark et a/. 1958a and 6) for this measurement with only 
one modification. The two electron counters were remounted so that the plane 
defined by their axes was now perpendicular to the neutron beam direction, 
instead of containing it as in the time reversal experiment. With this arrange- 
ment, the axes of the proton and electron counters satisfy the conditions of 
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Fig. 1. The same methods of beam purification and polarization were used for 
this experiment as for the time reversal experiment. As before, any effects of 
instrumental drifts were minimized by counting with hourly cycles of alter- 
nately polarized and unpolarized neutrons. These hourly cycles also served 


~ 


ee 
= 


Fic. 1. A neutron with its spin ‘“up’’ decays into an electron, EZ, a proton, p, and an 
antineutrino, ¥. With the proton and electron detected as in the left diagram the antineutrino 
must, by momentum conservation, have travelled approximately parallel to the neutron spin 
direction. With the proton and electron detected as in the right diagram the antineutrino 
must have travelled approximately opposite to the neutron spin direction. 


the useful purpose of providing a convenient monitor of the degree of neutron 
polarization. Daily checks were made of the range of energy sensitivity of the 
electron counting systems, and of the coincidence resolving time and of all 
other aspects of the electronics which could affect the observed number of 
coincidences. The neutron spins were turned over every few days by reversing 
the magnetic fields of the neutron polarizing and guide systems. Data were 
accumulated over several such cycles thus eliminating any asymmetries due 
to these fields or due to differences in the geometries of the two electron 
counters. 

The measured coincidence asymmetry, after correction for detection 
geometry, electron energy, and polarization of the neutron beam gives 
B = +0.96+0.40. The corrections include the effect of the .1(J/J).(p,. FE.) 
term in this measurement. For this latter correction which was small (—0.03) 
we have used the value of 4 determined experimentally by Burgy et a/. 

Since B = 2x(x—1)/(1+3x?) if « = C,4/Cy, our result implies a value of 
C,/Cy either less than —0.3 or greater than +6. The latest value for CZ. Ci 
derived by Kistner and Rustad (1959) from the half-life measurement of 
Sosnovsky ef al. (1959) is 1.42+0.06. These data are shown in Fig. 2 and 
clearly establish opposite phase for the V and A interactions in the neutron 
decay, a conclusion which has been reached independently by Burgy ef al. 
(1958a). An interaction of the form V—A is in agreement with the electron 
neutrino angular correlation results (Allen et al. 1959; Robson 1958; 
Trebukhovskii et al. 1959). 


We are grateful to Dr. R. Nathans for assistance in the production of the 
polarized beam, to Dr. B. N. Brockhouse, who suggested the use of a crystalline 
quartz filter for beam purification, and to Mr. G. C. Dixon for his assistance in 
constructing, testing, and using the experimental apparatus. 
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Fic. 2. Comparison of theory and experiment for the B coefficient. B is plotted as ordinate 
and x = C,4/Cy as abscissa. The curve gives the theoretical dependence of B on x. The vertical 
lines show the value of |C4/Cy| as determined from the neutron decay (Sosnovsky et al. 1959). 
The sectioned region encloses our experimental value for B allowing for the error quoted in 
the text. 
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THE SPECIFIC HEAT OF CRYSTALLINE QUARTZ BETWEEN 2°K AND 4°K 
G. H. S. Jones anp A. C, HoLitis HALLETT 


In 1956 a theoretical analysis of the lattice vibrations of natural quartz 
was made by Betts et al. and a value of 555.8° K was obtained for its Debye 
temperature. This value is considerably greater than the values of 166° K and 
204° IX deduced respectively by Anderson (1936) and by Saksena (1941) from 
experimental results, and it was of interest to measure the specific heat of 
natural quartz at low temperatures where there is reason to hope that the 
specific heat of quartz would obey the Debye 7* law thus enabling the Debye 
temperature to be determined from the measurements. 

Previous measurements at low temperatures have been made by Nernst 
(1911), Anderson (1936), and others and a compilation of these has been given 
by Sosman (1927). The lowest temperature at which measurements had been 
made was 25° K (by Nernst) and even at this temperature it was necessary 
to express the temperature variation of the specific heat of quartz as the sum 
of a Debye function plus a variety of Einstein functions; in Saksena’s (1941) 
analysis there are given 24 Einstein functions corresponding to a variety of 
frequencies discovered by spectroscopic methods. In all cases the Debye con- 
tribution to the specific heat was considered small compared with the Einstein 
functions, and there is probably considerable uncertainty in the values 166° K 
(Anderson) and 204° K (Saksena) given for the Debye temperature. 

Accordingly, measurements have been made at Toronto on a natural quartz 
crystal in the range of temperature from 2° K to 4° K and it has been found 
that the Debye 7* law is well followed by the results and a value of 470° K 
with an uncertainty of less than 10° K is deduced for the Debye temperature 
of quartz. This value is much closer to the theoretical estimate of 555.8° K 
(Betts et a/. 1956) than were the previous estimates from experimental results, 
but it would appear that some discrepancy still exists. 

The specimen used in this investigation was a piece of natural quartz 
crystal chosen for its freedom from internal defects as judged by visual 
inspection of its clarity. Its surface was not, however, perfect since the ends 
were somewhat chipped and one face was marked by a small roughened 
area. Another face had a complex step structure, the size of the steps being of 
the order of a few thousandths of an inch. It was approximately 4+cm long 
and was roughly Lem by lcm in section and had a mass of 10.926 ¢. To 
the specimen were bonded the ‘‘addenda” comprising a heater and carbon 
thermometer. The heater consisted of about 9 feet of Formel-coated Manganin 
wire of diameter 0.002 in. and had a resistance of about 600 ohms. This was 
wrapped around the crystal and cemented to it with a small quantity of shellac 
which was subsequently baked. The thermometer was a nominal 25-ohm, 
}-watt commercial radio resistor (Ohmite) whose outer casing had been ground 
away to reveal the carbon. This was attached to the specimen using a small 
quantity of cold-setting araldite. The total mass of the addenda, i.e. the 
heater, thermometer, and bonding agents, was 0.567 g. 
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The specimen was suspended within a brass chamber by means of fine nylon 
threads. These threads were secured to metal springs so that they were in 
tension, and the specimen hung freely in the center of the chamber without 
touching the walls or supports. The electrical leads from the thermometer 
and heater, coiled near the specimen, passed through a lead-in seal at the 
top of the chamber and through the liquid helium bath in which the chamber 
was immersed. This lead-in seal was a commercial seal of the type which had 
a number of fine tubes sealed in ceramic. It was soft-soldered into a hole 
bored in the top of the chamber, and the leads were passed through the 
tubes without join. The tubes and the socket for the seal were then filled with 
cold-setting araldite. The seal was vacuum tight even when immersed in 
helium IT. 

The chamber was supported at the bottom of the cryostat by a thin-walled 
monel tube which served as a pumping tube for the chamber and was con- 
nected to a pressure gauge, a liquid air trap, and diffusion and backing pumps, 
all outside the cryostat. A baffle inserted into the pumping tube near its 
lower end served to reduce the stray heat input into the chamber due to gas 
currents and to radiation. The resulting temperature drift of the specimen 
was somewhat variable depending primarily upon the efficiency of removal 
of the exchange gas within the chamber and upon the difference of tempera- 
ture between the sample and the outside helium bath. For all the results 
reported here, however, the drift rate was sufficiently small compared with 
the deliberate heating that corrections for it could be determined readily and 
reliably. 

The experiment consisted essentially of supplying a measured amount of 
power for a measured time to the heater, and observing the resulting change 
of temperature. The powers used varied between 3.7 and 0.22 microwatts 
and the heating times varied between 15 and 60 seconds. The power was 
determined by measuring the current flowing through the heater with a 
sensitive meter and the potential drop across the heater with a potentiometer. 
The change of temperature was observed as a change of resistance of the 
carbon thermometer using another potentiometer connected across it. These 
changes were measured either by observing the off-balance deflections of the 
potentiometer’s galvanometer or by balancing the potentiometer manually, 
depending upon the rate at which the temperature changed. The thermometer 
was calibrated against the vapor pressure of the helium bath using the 1948 
scale of temperature. Recalculation of the results in terms of the more recent 
1955L scale did not produce changes which were significant when compared 
with the random scatter of the results. 

The calculation of the heat capacity from the observed variation of tem- 
perature with time betore, during, and after the heating period tollowed 
standard methods. For those measurements at 4° K and above, good thermal 
isolation of the specimen was achieved and the temperature drift rates before 
and after heating were negligible. Thus in this region the heat capacity was 
calculated from the measured heat input and the difference between the 


initial and final temperatures. For measurements below about 2.5° K, good 
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thermal isolation could not be achieved and the temperature drift rates before 
and after heating were appreciable. However, the rates of change of tem- 
perature with time before and after heating were constant and equal, and 
the temperature change due to the heating could be determined readily. 
Between about 2.5° K and 4° K the temperature drift before heating was 
made very small by maintaining the outside helium bath temperature very 
close to the initial temperature of the specimen. The rate of cooling of the 
specimen after the heating period was not constant in this temperature range, 
and the method used to calculate the heat capacity was that described by 
Kouvel (1956). In all cases the temperature of measurement was taken to 
be the mean of the initial and final temperatures. It is felt that in general 
the accuracy achieved in measuring the heat capacity was about 2%, and the 
mean temperatures are known within about 2 millidegrees. 

The heat capacity measured by this procedure was the sum of the desired 
heat capacity of the quartz and the heat capacity of the addenda, com- 
prising the heater, carbon thermometer, and bonding material. Although 
attempts to measure the heat capacity of the addenda directly were not 
successful, an indirect determination was made by measuring the heat 
capacity of another specimen consisting of a rod of pure copper fitted with 
addenda as similar as possible in appearance to that on the quartz. Since the 
heat capacity of the copper alone was well known from the work of Corak 
et al. (1955), the heat capacity of the addenda could be determined by sub- 
traction. This heat capacity was surprisingly large, of the order of 40% of 
the heat capacity of the quartz alone, and was proportional to 7* within the 
errors of measurement. This last fact was used to calculate the heat capacity 
of the addenda at each of the temperatures at which measurements were 
made on the quartz specimen, and values of the heat capacity of the quartz 
alone were obtained by subtraction. 

The results obtained, expressed as the heat capacity of the quartz per 
gram, are plotted against the cube of the temperature in Fig. 1. It can be 
seen that the results lie reasonably well on a straight line and that therefore 
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Fic. 1. The heat capacity of the natural quartz specimen plotted against the cube of 
the absolute temperature. The slope of the straight line corresponds to a Debye temperature, 


8, of 469° K. 
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the specific heat of quartz is proportional to 7* in accordance with the Debye 
theory. Thus, assuming that this quartz can be represented by the chemical 
formula SiOz, it was possible to determine an effective Debye temperature, 
@ = 469° K, corresponding to this temperature range, and the straight line 
drawn in Fig. 1 has a slope corresponding to this value of 6. The scatter of 
the individual points around this straight line might be considered as a 
variation of @ with temperature, but since the departures of the individual 
results from the drawn line are more random than systematic, it is felt that 
these results give no definite information regarding a possible variation of @ 
with temperature. In view of the large contribution to the total heat capacity 
made by the addenda, it is difficult to estimate a figure for the uncertainty 
in the value 469° K quoted for 6. However, with all effects considered, the 
uncertainty in @ should not be greater than 10° Kk. 

Comparison of the present results with those obtained by Nernst (1911) 
and Anderson (1936) is shown in Fig. 2. Here the specific heat of natural 
crystalline quartz in calories per gram formula weight per degree has been 
plotted, using a logarithmic scale, against 7/0 also on a logarithmic scale 
assuming @ = 469° K from the present measurements. The full curve drawn in 
Fig. 2 is the Debye theoretical curve for an atomic solid equivalent in atomic 
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Fic. 2. Comparison of the results of this research with those of Nernst (1911) and Anderson 
(1936) and with the Debye theory (full curve). 
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mass and total number of atoms to the quartz. The full curve fits the present 
results since the value of 6 was derived from them. An exact fit with the 
results of Nernst and Anderson is clearly not present, but an exact fit could 
hardly be expected in view of the nature and structure of quartz. Figure 2 
does, however, show that the present results are in line with the earlier 
measurements, and also indicate to some extent the deviations of the behavior 


of quartz from the Debye theory. 


ACKNOWLEDGMENTS 
The authors are indebted to the National Research Council of Canada and 
to the Ontario Research Foundation for their grants-in-aid of research which 
made this work possible. One of us, G. H.S. J., is indebted to the California 
Research Corporation for financial assistance. 


ANDERSON, C. T. 1936. J. Am. Chem. Soc. 58, 568. 

Betts, D. D., Buatia, A. B., and Horton, G. K. 1956. Phys. Rev. 104, 43. 

Corak, W. S., GARFUNKEL, M. P., SATTERTHWAITE, C. B., and WEXLER, A. 1955. Phys. 
Rev. 98, 1699. 

KouveEL, J. S. 1956. J. Appl. Phys. 27, 639. 

Nernst, W. 1911. Ann. Phys. 36, 395. 

SaksENA, B. D. 1941. Proc. Indian Acad. Sci. A, 12, 93. 

SosMAN, R. B. 1927. Properties of silica (The Chemical Catalog Co. Inc., New York). 
See also International Critical Tables, 5, 105; 4, 21. 


RECEIVED JANUARY 30, 1960. 
DEPARTMENT OF PuysIcs, 
UNIVERSITY OF TORONTO, 
TORONTO, ONTARIO. 


ON THE AXIAL SYMMETRY OF THE NILSSON MODEL NUCLEUS* 
T. D. NEwTON 


The strong coupling approximation for Bohr’s collective model (1952) of 
the nucleus, particularly when the single particle states are realized as Nilsson 
(1955) functions, has had many successes (e.g. Mottelson and Nilsson 1955; 
Gove and Litherland 1958). It is usual to determine the best shape of the 
potential well of the single-particle Hamiltonian by seeking the minimum 
energy of the assigned number of particles as the well shape is modified. In 
the references quoted above it is always assumed that the nucleus possesses 
an axis of symmetry. It is therefore of some interest to inquire whether or 
not the spheroidal shape of the Nilsson (1955) single-particle model nucleus 
is intrinsically stable; that is, does the minimum energy of a given number 
of fermions in an anisotropic oscillator potential with Nilsson spin-orbit 
coupling always correspond to a potential with axial symmetry? 

A modification of the Bohr strong coupling model has recently been pro- 


posed by Davydov and Filippov (1958). In this model the collective motions 
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are represented entirely by rotations but the nuclear shape is permitted to be 
aspheroidal. This Davydov—Filippov model provides a simpler description of 
the states of even-even nuclei than does Bohr’s original form of the theory, 
although its application to other types of nuclei has not yet been fully worked 
out (cf. Davydov 1959). It is, therefore, also of interest for applications of 
this model to know what shape of oscillator potential will provide the minimum 
energy of a system of single particles. 

The energy levels and eigenfunctions of a completely anisotropic oscillator 
have been computed (Newton 1960). The single-particle Hamiltonian is 


h: 


Qu 


(1) H = — V+ 3ut(w, x +o, y +o, 2°)+Cl,.5+D1/. 
Here the /, are the infinitesimal pseudorotation operators; that is, if we 
define 


(2) ¢ Z wwe 2 n oo Bay 2 :* = u*we 2 
. feat jeer ; = 


then the /, are the infinitesimal rotations in (é,7,¢) space. C and D are con- 
stants. The mass parameter u* is the single-particle effective mass (Newton 
1959) of a nucleon in the nucleus. 

Since we wish to discuss (1) in cases in which w,, w,, and w, are all different, 
the single distortion parameter used by Nilsson is not sufficient and we have 
chosen to use parameters more nearly equivalent to those of Bohr (1952). We 


write 

(3a) = apace, wy = woA,, @; = weds, 
(3d) hw,ww, = I, 

(4a) A, = 1-\f5/42 6 cos(y+ 2/3), 

(40) A, = 1—-\f/5/4r B cos(y— 7/3), 

(4c) A, = 14)/5/42 B cos ¥. 


The frequencies*w,, ®,, #, are thus inversely proportional to the axial 
intercepts of the distorted nuclear surface defined by Bohr (1952). Equation 
(36) ensures constant nuclear volume and defines the energy scale. When 
y = Oor y = 7/3, equations (3) and (4) reduce to the Nilsson (1955) case, 
and the parameter 8 can then be related to Nilsson’s parameters 6 and e, 


B = 1.0569 (26+3y/[1+ (2/3)8}[1 — (4/3)5]—3) / (1-26) 
= 1.05696 (1 — (3 /2)6— 36°— (13 /24)6°) / (1 — 26) 


and 


B = 1.0569e/(1—e/3) 


for spheroidal systems. 
Having determined the single-particle energies, the total energy of a model 
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nucleus can be obtained. It is the sum over all occupied levels of the quantity 
(5) e& = 4(1+3u*/M)E,—(iu*/M)(Cl,.5+D1?), 


in which the subscript Rk labels a single-particle state of energy E, and (f), 
represents the expectation value of the operator f in the state k. M is the 
true nucleon mass. In the model, neutrons and protons are not differentiated; 
two particles of each kind are allowed in each state k. 

The Hamiltonian (1) is so constructed (Nilsson 1955, Appendix A) that 
the set of states of a given principal oscillator number N = n,+n,+n, is 
closed under /7. The eigenvalues E, have been computed for NV < 4, 35 states 
in all, for each of 46 values of the parameters 6 = —0.5 (0.1) to 0.5 and 
y = 0 (1/24) to 7/6. At B = 0, y is irrelevant. This provides positions for 
70 particles of one kind, but since at large distortions some of the NV > 5 
levels may be brought very low in energy we have not attempted to discuss 
systems containing more than 50 particles of one kind. 

The terms Cl,.8+D1? in (1) are primarily postulated for the purpose of 
adjusting the order of levels and making major shells of the model close at 
the true particle number. Nilsson considers the values 


(6a) C = —2khws, 
(60) D = —Kkpnyhw". 
On our energy scale fw» = 1. The parameter wy varies with the principal 


quantum number. Nilsson (1955) suggests wo = wi = we = 0, ws = 0.035, 
uy = 0.045, and x = 0.05. We have used the Nilsson values of D, but have 
varied C independently of D. 

We use the term ‘‘preferred shape” of a model nucleus containing Z protons 
and N neutrons to designate the potential of the form (1) for which the 
ground state energy is minimal under variation of 6 and y. The preterred 
shape can be specified by the values By and yo at the minimum. 

It is clear that the values of C and D affect the preferred shape. Consider 
a system of 40 nucleons of one kind. If C = 0, D = 0 we have a pure oscillator 
potential and 40 particles just fill the N = 3 shell. The preferred shape is 
spherical. If still with C = 0, D = 0 we consider 50 particles, we find the 
preferred shape is aspherical because of the partly filled N = 4 shell. How- 
ever, an acceptable shell model must include spin-orbit coupling with C < 0 
in order to make the shell closure occur at 50 particles. This term makes the 
greatest decrease in energy when the nucleus is spherical, and so as C decreases 
from 0 the preferred shape for 50 particles tends towards the spherical form. 


The preferred shape is not spherical at C = —0.1, but it has become spherical 
at C = —0.3. 

In Figs. 1 and 2 we show the approximate preferred values of (Soyo) for 
even sets of one kind of nucleon with C = —0.1 (Fig. 1) and C = —0.8 


(Fig. 2). These shapes possess the minimum energy, « = Le, from among 
the 46 calculated shapes. In Fig. 3 we show the way the energy of the svstems 
of 28 and of 50 nucleons varies with distortion when C = —0.8. It should 
be noted that for these calculations C was changed without changing D and 
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Fic. 1. The preferred distortion parameters with spin-orbit coupling C = —0.1 of systems 
containing an even number of one kind of fermion. The solid circles indicate the value of 8 
shown on the left ordinate scale; the crosses indicate the non-zero values of yo shown on the 
scale on the right side. If no cross appears for a given even number of particles the corre- 
sponding value yo is zero, i.e. the system prefers axial symmetry. 


the large value is not necessarily one which provides a good shell model 
level order. It is almost certainly too large to use in the shells of higher V. The 
purpose here is to indicate the form of the energy minimum near the spherical 
shape for these ‘‘magic numbers’? when the spin-orbit coupling is consider- 
ably increased. u*/ = 0.609 (Newton 1959) was used in computing «. This 
parameter has no very large effect on the properties now under discussion. 

The remaining figures illustrate our main point. They are ‘‘energy maps” 
of the ground state of five model nuclei. The systems illustrated are not 
particularly significant. The same general behavior is observable in greater 
or less degree in every system we have studied. 

The maps show the variation of the energy ¢« of the ground state as a 
function of 6 and y; the minimum of ¢« marks the preferred shape. It is more 
illuminating to restrict 8 to positive values and to represent the deformation 
(—8,y) by the equivalent parameters (+8, 7/3—y). Curves are drawn for 
each of the computed values of y, and are labelled p (for y = 0, prolate) 


1, 2,..., 7, o (for y = 7/3, oblate). At some points the variation of ¢ 
with y at fixed 6 is plotted as an insert. Between calculated points the curves 
are smoothly interpolated by eye. The contours are actually composed from 


ares with discontinuities of slope at the points where the lowest state changes 
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Fic. 2. The preferred distortion parameters with spin-orbit coupling C = —0.3. The 
meaning of the figure is otherwise the same as Fig. 1. 


from one configuration to another. Consequently the interpolations may be 
somewhat inaccurate. The general features of the energy surfaces are, how- 
ever, well established by the calculated values. 

The salient feature of the energy surfaces is the marked changes which 
occur with changing spin-orbit coupling. With C = 0 the preferred shape is 
usually fairly well defined but as C is decreased from zero the energy minima 
are broader. The preferred value 8) usually remains relatively well marked 
but the value of yo is often ill defined. If one requires that the closed shells 
of 28 and 50 particles should be spherical, then the spin-orbit coupling in 
the .V 3 and N t states should be negative and greater in magnitude 
than the value —0.1 hws. 

Our few examples suffice to show also that the preferred shape for some 
sets of particles is not axially symmetric. One can see from ligs. ta-4e that 
the energy of the system of independent particles may be decreased by 10 
to 15° ¢ of hwy by the change from axial symmetry to a non-axial shape. This 
energy change is comparable to the energy of a collective mode of excitation. 
Consequently the effects of departure from axial symmetry in Bohr’s (1952) 
strong coupling theory cannot properly be treated by perturbation methods. 

kinally we remark that the preferred shape is not always well defined, i.e. 
the ground state energy sometimes varies but little over a wide range ol 


distortion parameters. Thus it is difficult to justify the concept of an asym- 


metrical top with unvarying principal moments of inertia. 
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Fic. 3. The ground state energy of systems of 28 and of 50 particles of one kind when 
C = —0.3 as a function of distortion. The abscissa is the parameter 8 and the function E(3) 


is drawn for each available value of y from » = 0 (prolate shape, curve p) to y = r/3 (oblate 


shape, curve a) in steps y = w/24. 


In fact in the region of the light nuclei both the Bohr strong coupling 
approximation (e.g. Litherland, Paul, Bartholomew, and Gove 1956) and 
the Davydov—Filippov model (e.g. Litherland 1960) describe well some of 
the features of the excited states but both approximations have some dis- 
agreements with experiment, usually in predicting transition probabilities. The 
present calculation indicates some of the difficulties in the justification of 
these models, which may be the reason for their incomplete agreement with 


experiment. 
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225.5 





Fic. 4. The ground state energy of various systems of two kinds of particle as a function 
of the distortion. Each energy map is shown for three values of the spin-orbit coupling con- 
stant C. The parametrization is the same as in Fig. 3. 
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ON THE CONDENSATION COEFFICIENT IN THE GROWTH 
OF SILVER FROM THE VAPOR PHASE 


R. A. Rapp,* J. P. Hirtu,t AND G. M. Pounpt 


Chandra and Scott (1958) have recently reported condensation coefficients 
which were appreciably less than unity for the growth of gold, silver, and 
copper from the vapor at high beam fluxes. For silver at fluxes of about 10" 
atoms/cm? sec, the condensation coefficients reported by Chandra and Scott 
ranged between 0.982 at — 28° C and 0.879 at 171° C. Frauenfelder (1950) and 
Devienne (1953) had previously reported condensation coefficients for silver to 
be 0.3 to 0.8 and 0.9 to 1.0, respectively. Yang, Simnad, and Pound (1954) 
have determined the condensation coefficient for silver to be 1.0 at low fluxes 
of 10'! atoms/cm? sec at 192° C. Although different silver fluxes were used in 
the separate investigations, these data are still contradictory. 

In view of this contradiction, the present authors determined condensation 
coefficients for silver at relatively high fluxes of 1.3-4.5X 10" atoms cm? sec 
In these experiments, the 


. 


and temperatures between 45°C and 440°C. 
residual gas pressure in the vacuum system was varied between 4.4 X 10~° and 
6X10-" mm Hg in order to determine whether the adsorption of gaseous 
impurities on the growing silver substrate affects the value of the condensation 
cceificient. The apparatus, shown in Fig. 1, is similar to that used by Chandra 
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Fic. 1. Silver condensation chamber. 


and Scott. A glass cover slide, 25 mm square and 0.15 mm thick, was held 

snugly in a copper jig which was heated by tungsten filaments. The tempera- 

tures of the copper block holding the “‘hot plate’”’ (referred to as the “‘retlector 
*Max-Planck Institut fiir Physikalische Chemie, Gottingen. 
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plate” by Chandra and Scott) were measured with a chromel—alumel thermo- : 
couple and are reported in Table I. The actual “front face” temperatures of 
the glass hot plates may have been somewhat lower than these listed tempera- 
tures. For a copper supporting block temperature of 200°C, calculations 
showed that the heat from radiation sources alone will maintain the silver 
deposit on the hot plate at 140°C. Since the amount of heat received by 
conduction from the copper block cannot be estimated (exact area of contact 
is unknown), one must conclude that the temperature of the silver deposit 
was between 140° C and 200° C. However, the indeterminacy of the substrate 
temperature does not affect the calculated condensation coefficients or the 
conclusions from this investigation. A silver bead which was held in a tungsten 
foil boat and heated to 900 to 950° C served as the source of silver vapor for 
deposition on the hot plate. Facing the hot plate, but not the silver beam, was | 
another identical glass plate called the ‘‘cold plate’ (referred to as the ‘‘collector 
plate’’ by Chandra and Scott). A ‘‘cold plate” temperature of approximately | 
—78° C was maintained by contact between the plate and a copper block 
held at —78° C. The function of the cold plate was to collect even the smallest 
flux of silver atoms which might reflect or re-evaporate from the growing 
silver substrate (hot plate). In some experiments, the initial glass hot plate 
was gold-plated to insure rapid initial nucleation of the silver on the hot plate. 
Maintenance of all the surfaces of the system, except the hot plate, at —78° C 
precluded re-evaporation of any silver atoms to the cold plate from these | 
extraneous surfaces. This was a notable departure from the radiation chamber 


Urs i Sol then ean anita 


apparatus of Chandra and Scott. 

For the ultra-high vacua experiments, cleanliness of the substrates from 
condensable impurities was insured by a 3-cycle bake out of the system at 
450° C, double-trapping of the Octoil pump with liquid nitrogen, and the cold 
bath (—78° C) surrounding the condensation chamber. For the two higher 
pressure experiments, a single liquid nitrogen trap and the surrounding cold 
bath were used. 

Prior to an experiment, a flux shutter operated by a magnet outside of the 
vacuum system shielded the hot plate from metal vapor flux while the silver 
source was outgassing. When the residual gas pressure and the temperatures 
of the hot plate, cold plate, and silver source had become steady, the shutter 
was rotated, and silver was allowed to deposit on the hot plate. Only about 
10 seconds was required to deposit a monolayer of silver on the hot plate for 
a flux of 1.43 10'* atoms/cm? sec. Hence the condensation on the hot plate 
was essentially silver on silver for an experiment of 1-hour duration. 

The value of the condensation coefficient was determined by quantitative 
chemical analyses of the silver deposits on both the hot and cold plates. A 
calculated geometric factor related these data to the amounts of silver which 
would have been found on infinitely large parallel plates. The calculated 
geometric factor for our arrangement of plates was 0.29. The condensation 


coefficient was calculated from the formula 


(1) a@ = wt. Ag hot plate/ wt. Ag hot plate + (wt. Ag cold plate/0.29). 
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A description of the analytical methods is given by Sandell (1944). The 
small amounts of silver which were found on the cold plates were determined 
by using a Beckman D. U. spectrophotometer. 

The results of this investigation are presented in Table I. In this table, the 
error in the hot plate weight determination is +2X10~7 g. In a set of blank 
runs it was found that on the average both plates picked up 1.6X10~ g of 
silver during bake out. This weight was subtracted from both hot and cold 
plates for the baked out tests at 45°, 202°, and 440° C and a@ values determined 
from the weights thus corrected. The estimated error includes only the un- 
certainties in the original weights of silver. The maximum error includes the 
error in original weights and the factor 1.610-* g for the cases where this 
latter factor was subtracted from the raw data. Because of the necessity for 
the bake out each test was 1 week in duration. Hence insufficient tests could 
be performed to determine the reproducibility of the data. 

In the experiments at 200° C the condensation coefficient was unity within 
experimental error. In all cases, the small amount of silver found on the cold 
plate is felt to represent silver re-evaporated from the hot plate in the initial 
portion of a run prior to or during nucleation of the first silver layers on the 
hot plate. Thus the apparent values of @ for the tests at 45°C and 440° C 
are lower than the values at 200°C only because of the smaller weights 
deposited on the hot plate. The adsorption of residual gases at 199° C and 
205° C did not decrease the value of the condensation coefficient from unity. 
Although the experimental conditions of these experiments were greatly 
different from those of Yang, Simnad, and Pound, our data are in agreement 
with the results of their investigation. Our data are in disagreement with the 
results of Chandra and Scott, and Frauenfelder. Explanations for the low 
condensation coefficients reported by these authors can only be speculative, 
but, as a result of our work, we would suggest that the seemingly re-evaporated 
atoms in the other investigations were from some source other than the 


growing silver substrate. 


This work was sponsored by the U.S. Office of Naval Research, and is 
based on a portion of a thesis submitted by R. A. Rapp in partial fulfillment 
of the requirements for the degree of Doctor of Philosophy to the Graduate 
Faculty, Carnegie Institute of Technology. 
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Under this heading brief reports of important discoveries in physics may be published. These 
reports should not exceed 600 words and, for any issue, should be submitted not later than six weeks 
previous to the first day of the month of issue. No proof will be sent to the authors. 


Frequency Measurement of Standard Frequency Transmissions! :? 
Measurements are made at Ottawa, Canada, using N.R.C. caesium-beam frequency resonator 
as reference standard (with an assumed frequency of 9 192 631 770 c.p.s.). Frequency devia- 
tions from nominal are quoted in parts per 10!°. A negative sign indicates that the frequency 


is below nominal. 


GBR, 16 ke/s 


Date, —_— 

February MSF, 7}-hour 24-hour KK2XEI, 
1960 60 ke/s average* average 60 ke/s 

1 —177 —173 —177 N.M. 

2 —179 —177 —175 —8s9 

3 —174 —175 —176 —S88 

+ —171 —174 —175 —s8y 

5 N.M. —173 —174 —88 

6 173 —179 —176 N.M. 

7 —184 —179 — 180 N.M. 

8 — 198 —178 —179 —87 

9 N.M. N.M. N.M. —96 

10 — 187 N.M. N.M. —97 

11 — 190 N.M. N.M, N.M. 

12 -171 N.M. N.M. —105 

13 N.M. N.M. N.M. N.M. 

14 —182 N.M. N.M. N.M. 

15 173 N.M. N.M. — 94 

16 N.M 178 N.M. —101 

7 sy —175 —178 —92 

18 185 — 182 -177 —93 

19 170 — 182 - 180 —o+ 

20 183 —178 - 183 N.M 

21 186 —183 176 N.M. 

$3 186 183 - 183 N.M 

23 sy 179 183 96 

24 175 184 LSS —4 

25 N.M. ISS 1s4 —96 

26 N.M. 1St 186 36 

27 195 194 188 N.M 

28 197 193 193 N.M 

20 N.M. igo 192 —45 

Midmonthly mean 1838 180 180 8S 

Midmonthly mean 
o WWYV Lid 


Note: N.M 


*Time of observations 
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